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PHOTONEUTRON REACTIONS IN C® AND O'! 


J. P. Roatsvic,? IsHwar C. Gupta, AND R. N. H. HAsLAm 


ABSTRACT 


Absolute yields of the reactions C'%(y,n)C" and O'(y,2)O have been deter- 
mined at, 22-Mev maximum bremsstrahlung energy using the University of 
Saskatchewan 24-Mev betatron. For the reaction C'%(y,n)C" a yield curve 
from threshold to 24 Mev has been obtained and the cross-section curve for the 
reaction computed. A thorough comparison with other results has been made. 


INTRODUCTION 

The reactions C!(y,n)C" and O'*(y,2)O" have been investigated by several 
authors: Price and Kerst (1950), Haslam et al. (1951), Montalbetti e¢ al. 
(1953), Nathans and Halpern (1954), Katz et al. (1954), Barber et al. (1955), 
Cook (1957), and Thorson (1957). 

However, in view of the discrepancies of the results obtained (see Tables I 
and II) a redetermination of the activation and the cross-section curve for 
C(y,2)C" was thought desirable. It was of particular interest to find out 
whether most of the reaction cross section took place by photon absorption 
into discrete nuclear levels or into the giant dipole resonance. The results in 
the case of the O'*(y,2)O" reaction also cover a wide range. Due to this dis- 
crepancy and since the technique of absolute 8 counting is now a little better 
understood (Roalsvig and Haslam 1959) it was considered desirable to re- 
measure the absolute yield of the O'%(y,z)O reaction. The values obtained 
in the present work are thought to be of greater accuracy than those published 
from this department before, due to greatly improved knowledge of the 
betatron energy scale at the present time. 


EXPERIMENTAL TECHNIQUE 


(a) The carbon samples consisted of cylinders of polystyrene, 7.5 cm long 
with an outside diameter of 2.6 cm and a wall thickness of 2 mm. The samples 
were irradiated in the center of the X-ray beam from the University of 
Saskatchewan 24-Mev betatron. They were placed in a fixed position 60 cm 
from the betatron target, with their cylindrical axes perpendicular to the 
beam direction. 


'1Manuscript received January 17, 1961. 

Contribution from the Physics Department, University of Saskatchewan, Saskatoon, 
Saskatchewan. 5 s 

2Present address: Physics Department, St. John’s University, Jamaica, New York, U.S.A. 
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The secondary activity from C" with a half life of 20.5 minutes was counted 
with a Victoreen IB85 cylindrical Geiger—Miiller counter placed inside a thick 
lead castle. The 8+ activity was counted in the time interval between 10 and 
20 minutes after the irradiation was stopped. Any activity due to O'”, which 
decays with a half life of 2.1 minutes, was then found to be negligible. 

(6) For the determination of absolute yields of O'%(y,2)O and C#(y,2)C", 
samples of mylar polyester (CisH12O¢) in the form of circular disks of 2.5 cm 
diameter and various thicknesses were irradiated in the center of the X-ray 
beam. The samples were placed at a fixed position of 18.3 cm from the target. 
The irradiation time was 2 minutes for studying the O'%(7,7)O" reaction and 
15 minutes for the study of the reaction C?(y,2)C". 

The samples were then mounted on different backscatterers (C, Al, Cu, Cd, 
and Pb) and the total secondary activity produced (i.e. O with a half life 
of 2.1 minutes and C" with a half life of 20.5 minutes) was measured. The 
counting was done using a proportional flow counter and the activities were 
followed for several half lives. 

From the plots of counts/minute versus time, the C" activity (which forms 
the background for the O" activity) was obtained and subtracted from the 
total activity, yielding the O" activity produced in the sample. 

The counting rates were corrected for background and for counting losses 
due to the dead time of the tube. In the case of the proportional flow counter, 
the counting loss is less than 0.1% for a counting rate of 200,000 counts/ 
minute. No correction for counting loss was necessary, since the counting rate 
was usually in the neighborhood of 10,000 counts/minute. In case of counting 
by the Geiger tube, the counting loss is relatively higher and therefore the 
following formula was used for this correction. 


x, = t- I) (n-n) 
fn Ap(71—e M70) 


where Ky is the corrected number of counts, 
Ky is the measured number of counts, 
d is the disintegration constant (0.0338 min-! for C"), 
p is the dead time of the tube, 
T, = el, ro = e~", where t, and ¢, are 10 and 20 minutes respectively. 
The energy stability of the betatron was checked regularly by irradiating 
samples of boric acid at a fixed integrator setting (called monitor point) and 
measuring the induced 2.l-minute 8+ activity from O', During the course 
of the experiment the total energy drift was about 50 kev, but the integrator 
setting was correspondingly adjusted. The maximum unaccounted for drift 
was about 9 kev. 
Two slightly different procedures were used to obtain relative activation 


curves: 

(a) The samples were given a known relative dose, measured by a monitor 
in connection with a condenser (see Johns et al. 1950). Values for roentgens 
per unit relative dose for different maximum photon energies were obtained 
using a Baldwin—Farmer type r-meter. Thus relative numbers of counts per 
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roentgen could be obtained. These were corrected to relative saturated activity 
by multiplication by the factor tg{1—exp(—A¢g)}~!, where ¢p is the irradiation 
time and ) is the disintegration constant. 

(6) All samples were irradiated for the same length of time. The relative 
dose was measured by a monitor in connection with an R-C circuit, whose 
time constant equaled the mean life of the activity under consideration. 
Relative values of roentgens per unit dose were obtained by using a Baldwin- 
Farmer type r-meter. 

The two relative activation curves for C!*(y,n)C" reaction when normalized 
at 22 Mev showed good agreement, deviations being within statistical un- 
certainty. 


NORMALIZATION AND RESULTS 


(a) The relative yield curve for the reaction C(y,n)C" was normalized 
by using the reaction Cu®(y,z)Cu® as a standard. Circular disks of copper 
and polystyrene were used, with diameter 2.5cm and thicknesses ranging 
from 10 mg/cm? to several hundred mg/cm?*. All disks were irradiated at 
22-Mev maximum photon energy, the relative dose being measured with a 
monitor as under (a) above. The secondary activities were counted in a pro- 
portional flow counter and followed for several half lives. The number of 
counts/minute at ¢ = 0 were corrected for finite irradiation time and divided 
by the number of g-atoms of C” and Cu® respectively in the samples, and 
also by the number of relative dose units given to the disks. The correction 
factor for the self-absorption in the sample, backscattering, and absorption 
of the backscattered 8+-particles was calculated using the empirical formulae 
given by Roalsvig and Haslam (1959). 

The ratio of the saturated activity of the reaction Cu®(y,z)Cu® in neutrons 
per g-atom per relative dose unit at 22 Mev to that of the reaction C(y,m)C" 
was found to be 196+2. The absolute yield value for the former reaction 
has recently been determined in this laboratory (Roalsvig, Haslam, and 
McKenzie 1959). Uncorrected for absorption of X rays in the lucite block 
surrounding the r-meter thimble chamber, this value was found to be 2.42 x 108 
n/g-atom per 100 r at 22 Mev. The uncorrected yield value for the reaction 
C!2(y,2)C" thus becomes 1.23 X 10° n/g-atom per 100 r. The correction factor 
for absorption in the lucite block was calculated to be 1.070 for the latter 
reaction, thus giving an absolute yield value of 1.15X10° n/g-atom per 100 r 
for the reaction C!?(y,n)C"™ at 22 Mev. The absolute yield values quoted in this 
paper are those that would be obtained if the sample were placed in the 
center of an 8-cm lucite cube in the same position as the r-meter. These 
values are approximately 7% lower than if the sample were placed in air at 
the same distance from the betatron target. 

(6) The absolute yield of C(y,2)C"™ at 22 Mev was also determined without 
using the empirical relations for the correction factors for self-absorption in 
the sample, backscattering, and absorption of backscattered 8+-particles. The 
experimental technique followed was similar to that of Roalsvig and Haslam 
(1959). The dose in this case was calibrated against a Victoreen thimble in 
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an 8-cm lucite block. Dose calibration was checked by determining the 
Cu®(7,2)Cu® yield also. 

For the determination of the absolute yield, the counting rates after sub- 
tracting background and any other undesirable activity were plotted against 
time. From these plots the counting rate at ¢ = 0 after irradiation was 
obtained. This quantity was divided by the number of g-atoms of C” and the 
relative dose (which was already corrected for the decay of the observed 
activity in the sample during the irradiation). The values of this new quantity, 
called N’(x,z), where x refers to the thickness of the sample and z to the 
atomic number of the backscatterer, were plotted against 2 for each value of 
x and against x for each value of z (see Figs. 1, 2, 3, and 4). From the extra- 
polation of these curves the values of N’(0,0) were obtained. The absolute 
yield was obtained by using the average value of N’(0,0), multiplying it by 
the geometry factor (for the system used, the geometry factor was 2), and 
dividing by the number of 100 roentgens per unit registered dose (for further 
details, see Roalsvig and Haslam 1959). 

The absolute vield curve, normalized to the value 1.15 10° n/g-atom per 
100 r for the reaction C!(y,7)C" at 22 Mev, is given in Fig. 5. The cross-section 
curve was calculated by using the ‘‘modified total spectrum method”’ of Pen- 
fold and Leiss (1958) and is shown in Fig. 6. The cross-section curve is seen 
to have its maximum at 23.0 Mev with omax = 7.9 millibarns. The half width 
of the curve, Ty, is 3.2 Mev, and the integrated cross section from threshold 
to 24 Mev is 0.022 Mev barn. 

(c) The absolute yield of the reaction O'%(y,7)O was determined in the 
same way as outlined in Part (4) of this section. Three independent runs of 
the experiment were made. The average value for the absolute yield for 
O'%(y,2)O" reaction at 22 Mev (maximum energy) X rays is (2.837+0.2) x 10° 
n/g-atom per 100 r. 


DISCUSSION 


(a) The Yield Values 
The measured yield values at 22 Mev are given in Table I. The absolute 











TABLE I 
Absolute yields of O'%(y,7)O% and C%(7,n)C" at 22 Mev 
Absolute yield of Absolute yield of 
O'(,)O"8 C#(y,n)CH 
References (n/g-atom 100 r) (n/g-atom 100 r) 
Price and Kerst (1950) 0.67 X 108 0.67 X 10° 


at 


Haslam, Johns, and Horsley (1951) 0.75 X 108 


Johns, Horsley, Haslam, and Quinton (1951) f 1.68 X 10° — 
.2X108 


Montalbetti e¢ al. (1953) 2.7 X10° 3 

Nathans and Halpern (1954) — 2.8106 
Barber et al. (1955) — 2.3X105 
Cook (1957) — 1.5X10% 
Thorson (1957) — 0.76 X 10° 
Present work 2.37 X108 1.15X10* 





“According to the note added in proof the author mentions that this value should be reduced by 10%. 
A redetermination of the absolute yield gave results in close agreement with this value. 
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Fic. 1. Relative counting rate N’(x,z) versus sample thickness, x, for a fixed value of the 
atomic number of backscattering material (for C™ positrons). 

Fic. 2. Relative counting rate N’(x,z) versus z, the atomic number of the backscatterer, for 
a fixed value of x, the sample thickness (for C" positrons). e and d designate two typical 
sample thicknesses. 
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Fic. 3. Relative counting rate N’(x,z) versus sample thickness, x, for a fixed value of the 
atomic number of backscattering material (for O" positrons), 

Fic. 4, Relative counting rate N’(x,z) versus z, the atomic number of the backscatterer, for 
a fixed value of x, the sample thickness (for O" positrons). e and d designate two typical 


sample thicknesses. 
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Fic. 5. Absolute activation curve for the reaction C!%(y,2)C". 
Fic. 6. Cross-section curve for the reaction C%(y,n)C™. 
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yield values of the C(y,z)C" reaction show a large discrepancy, varying from 
0.67 X 10° n/g-atom per 100 r to 3.2 X 10° n/g-atom per 100 r, that is by a factor 
greater than 4. To investigate this great discrepancy, the activation curves 
from the investigations in which the photon source was a betatron were drawn 
on the same figure, Fig. 7. From this figure it is clearly seen that the shapes 
of the activation curves are generally quite similar, the curves being, however, 
displaced from each other along the energy axis. Thus, for example, a dis- 
placement of the curve of Montalbetti et al. of about 1.35 Mev to the right 
would make it coincide with the curve of the present work. This may thus 
indicate that the main source of error in obtaining the absolute yield at the 
standard point 22 Mev is in the energy calibration of the betatron. Due to 
the rapid variation of the yield with energy, a rather small displacement in 
energy will give a high difference in the absolute yield. 
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Fic. 7. Activation curves for the reaction C'%(y,")C", as given by different authors. 
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Our value of the absolute yield of O'8(y,7)O" at 22 Mev is 2.37 X 10° n/g-atom 
per 100 r. This is twice as great as our value for the yield of the C?(y,”)C" 
reaction. Price and Kerst (1950) and Montalbetti et a/. (1953) report nearly 
the same yields for the two reactions. However, if as suggested earlier, in 
the curves of Montalbetti et al., the energy scale is shifted by 1.35 Mev, the 
ratio of the yields of O'%(y,z)O to C®(y,n)C" at 22 Mev becomes ~ 1.7. 
The ratio of the yields obtained in our experiment is in very good agreement 
with those of Haslam et al. (1951) and Johns e¢ al. (1951). Carver and Lokan 
(1957) report the yields of C#(y,v)C" and O'*(y,2)O" in arbitrary units as 
0.8 and 1.5 respectively. If the arbitrary units are the same in both the 
reactions then this result is also in agreement with ours. Further checks on 
the C(y,2)C" yield were made by determining the ratio of the yields of 
O'%(y,2)O" to C?(y,2)C", which was found to be 2.28. Using our value of 
oxygen yield-we get for the C"(y,2)C" reaction at 22 Mev a yield of 1.04108 
n/g-atom per 100 r. 

The similarity in the shape of the activation curves in C?(y,z)C" suggests 
that the cross-section curves would also exhibit similar shapes, though dis- 
placed with respect to the scale of photon energy. Thus the variation in 
Emax, the energy at which the cross-section curve has its maximum value, 
should provide a measure of the variation in the energy scale. The greater 
the apparent energy Emax, the lower should be the yield at some fixed energy 
below Emax, for example, 22 Mev. That this is the case is shown quite clearly 
in Fig. 8, in which the yield values reported for a bremsstrahlung energy of 
22 Mev are plotted against Emax. 

If the shapes of the activation curves had been exactly similar, and if the 
same method had been used by all authors to calculate the cross-section 
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Fic. 8. Relation between the yield at 22 Mev and Emax, the energy at which the cross- 
section curve has its maximum value. 
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curves, the values of the yields at E,,, should have been the same. These 
vield values are given in Table II, column 38. Although there is some variation, 
from 1.95X1C®n/g-atom per 100r to 3.45X10%n/g-atom per 100r, the 
spread is much lower than for the yield values at 22 Mev. It would, therefore, 
seem more valuable, while energy scales are in doubt, to use the energy Emax 
as a comparison point for the absolute yield, rather than an arbitrary energy 
value such as 22 Mev. 

From Table II, it is seen that the average value of E,,x for the C!(y,2)C"™ 
reaction when taking all results into consideration is 22.4 Mev with a standard 
deviation of 0.6 Mev. Correspondingly from Fig. 8, one might deduce a ‘“‘mean 
value” of the yield at 22 Mev of 2.4X10°n/g-atom per 100 r with an un- 
certainty of +0.5X10° n/g-atom per 100 r and —1.4X10* n/g-atom per 100 r 
for this reaction. However, the actual mean value of the yields for C!?(y,2)C™ 
at 22 Mev is (1.6+1.0)10°n/g-atom per 100 r. The relatively low value 
of the yield of C(y,n)C" at this energy obtained in the present work, 1.15 x 10° 
n/g-atom per 100 r, corresponds to a relatively high value of Emax. 


(b) The Cross Sections 

The integrated cross sections from threshold to an energy E; as reported 
by the different authors are given in Table II, column 6. The different values 
of the energy E, are also given in parentheses in Mev. In order to get com- 
parable values of | ¥1.0(y,2)dE, a constant value of E, has to be used, the 
value 24 Mev has been chosen here. The integrated cross sections from thresh- 
old up to this energy have been tabulated from the reported cross-section 
curves, the values being given in Table II, column 7. The agreement between 
these values is rather good, the average value is found to be 0.026 Mev barn 
with a standard deviation of 0.003 Mev barn. 

The difference in the values is mainly due to inaccuracy in computing the 
cross section from the measured activation curves. If this procedure had 
been more accurate, decreasing values of S24 Mevo(y,n)dE should be expected 
to appear with increasing values of Emax. No such relationship is, however, 
indicated by the reported values as seen from Fig. 9a. Also the values of 
Cmax, the maximum cross section, should be expected to be constant. As seen 
from Table II, column 4, these values have, however, a rather high standard 
deviation of about 25%. 

If the error introduced by computing cross-section curves was independent 
of the energy scale, the ‘“‘normalized integrated cross section’, defined as 
S24 Mevo(y,)dE/omax and given in Table II, column 8, should show decreasing 
values with increasing Eas. This is generally not the case, as seen from 
Fig. 96. Thus a distortion with energy does also take place, which might also 
be seen from the divergent values of the half widths Ty as given in Table II, 
column 5. A more accurate method of computing cross section from the 
activation values is therefore highly desirable. 

However, one might say that to a first approximation, at least up to Enaz, 
the shape of the cross-section curve is fairly well established. This is also 


indicated by Fig. 10, which shows the decrease of [24 Mevo(y,n)dE with decreasing 


values of omax: 
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Fic. 9. (a) Relation between the integrated cross section from threshold to 24 Mev, and 
Emax, the energy at which the cross-section curve has its maximum value. 

(6) Relation between the ‘normalized integrated cross section” from threshold to 24 Mev, 
and Emax, the energy at which the cross-section curve has its maximum value. 


(c) Cross Section under the Breaks 

The integrated cross section under the breaks in the fine structure of the 
activation curve Lt MevGfodE has been investigated by Katz et al. (1954) 
and Thorson (1957), the values being given in Table II, column 9. These 
two values are, however, not directly comparable, since the value of Katz 
et al. is obtained by using the yield value of Montalbetti et al. of 3.2X10° 
n/g-atom per 100 r at 22 Mev, while Thorson determined this value to be 
0.76 X 10° n/g-atom per 100 r. 

When both these integrated cross-section values are normalized according 
to the yield value of 1.15 10° n/g-atom per 100 r as obtained in the present 
work, they become 0.019 Mev barn and 0.013 Mev barn respectively. The 
average value 0.016-+-0.003 Mev barn should then be compared with the value 
0.022 Mev barn for the continuous integrated cross section. This should 
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Fic. 10. Relation between the integrated cross section from threshold to 24 Mev, and the 
maximum cross section. 


thus indicate that a substantial fraction, namely 75+20% of the C?(y,”)C"™ 
reaction cross section is due to photon absorption into discrete nuclear levels 
located at the energies of the observed breaks in the activation curve. 


SUMMARY 


The absolute yield for O'%(y,72)O” reaction has been determined and is 
found to be (2.37+0.2) X10° n/g-atom per 100 r. The absolute yield and the 
activation curve for the reaction C!(y,n)C" have been measured and the 
corresponding cross-section curve calculated. These have been compared with 
previous investigations, and reasons for the discrepancy in absolute yields 
and cross sections have been discussed. By taking all works on the reaction 
C(y,n)C" up to the present, average values of various parameters are: 


(a) yield at 22 Mev: (1.6+1.0) 10° n/g-atom per 100 r, 


(6) Beness 22.4+0.6 Mev, 
(c) yield at Enaz: (2.6+0.3) X 10° n/g-atom per 100 r, 
(2Y ima? 10.8+2.5 millibarns, 
(ce) Pa: 3.1+0.8 Mev, 

224 Mev 
(f) J o(y,n)dE: 0.026+0.003 Mev barn, 

thresh 

24 Mev 2 

(g) XG J o dE 
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The values obtained in the present work are: 
(a) yield at 22 Mev: (1.15+0.02) X 10° n/g-atom per 100 r, 


(b) “Bias? 23.0 Mev, 

(c) ‘yield at-Eaax: 2.77 X 10° n/g-atom per 100 r, 
(2) Does? 7.9 millibarns, 

i) ‘Fe: 3.2 Mev, 


24 Mev 
(f) f a(y,n)dE : 0.022 Mev barn. 
t 


hresh 
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THE ELECTRICAL CONDUCTIVITY OF GELATIN FILM 
HUMIDIFIED WITH HEAVY WATER VAPOR! 


C. D. NIVEN 


ABSTRACT 


When protons were replaced by deuterium nuclei in gelatin films, by exposure 
of dry films in 75% relative humidities, the resistance increased very markedly. 


At the present time proton migration is believed to account for the electrical 
conduction of hydrophilic substances. This conduction is much smaller than 
that of a true semiconductor and Riehl (1956) has suggested that a mechanism 
depending on the formation of chains of oxygen atoms, each of which has a 
proton suitably located to move to the next oxygen atom in the chain, explains 
the migration of the protons in a field. 

If protons are responsible for the conduction in these hydrophilic substances, 
and, as Dr. Ta-You Wu* pointed out to the writer, if every proton taking 
part in the conduction process is replaced by a deuterium nucleus, the resistance 
should increase markedly. The work described below was undertaken to deter- 
mine whether this was so, because if there was no increase in resistance the 
proton-conduction concept was probably incorrect: on the other hand if there 
was an increase, this concept was confirmed. 

The problem demanded the provision of two comparable atmospheres, one 
humidified with ordinary water and the other humidified with heavy water. No 
relative humidity tables were available for use with a wet and dry bulb 
psychrometer when heavy water vapor instead of ordinary water vapor was 
in the air. Accordingly some method had to be found to determine the relative 
humidity from wet and dry readings using heavy water. Since the results of 
the conductivity measurements are highly dependent on the reliability of this 
method and on the figures actually used in applying it, some details concerning 
this calculation may not be out of place. 


CALCULATION OF RELATIVE HUMIDITY FROM WET AND DRY 
BULB READINGS IN AN ATMOSPHERE HUMIDIFIED WITH 
HEAVY WATER VAPOR 

The method was that described by Hougen and Watson in their book 
Industrial Chemical Calculations (1936). The basic assumption is that the heat 
lost by the humid air must equal that gained by the water in vaporization 
and superheating: hence an expression is obtained which can be written in 
the form 

‘Manuscript received January 11, 1961. 


Contribution from the Division of Applied Physics, National Research Council, Ottawa, 


Canada. 
Issued as N.R.C. No. 6255. 
*Chief, Theoretical Physics Section, Division of Pure Physics, National Research Council 


of Canada. 
Can. J. Phys. Vol. 39 (1961) 


Oe 


657 








CANADIAN JOURNAL OF PHYSICS. VOL. 39, 1961 


_ LyXHw—6.95D 


pe 
L,+8.4D 
where H is the molal humidity of the air at dry bulb temperature ¢ and 
wet bulb temperature t,, 
D = t-ty, 


L, = molal heat of vaporization at ty, 

H, = molal humidity of saturated air at f,, 
and the two constants 6.95 and 8.4 are the heat capacities of air and of water 
vapor respectively. These, incidentally, can be taken as constant, but the 
assumption is made that the heat capacity of heavy water vapor is approxi- 
mately equal to that of ordinary water vapor. As a matter of fact in the 
present calculation the whole contribution of 8.4D to the denominator is less 
than 3 parts in 1000: accordingly this assumption does not affect the results. 

The particular relative humidity at which the heavy water readings were 
made corresponded to dry and wet bulb readings of 87° and 81°F. (The 
Fahrenheit degrees enter in because the psychrometer was equipped with 
Fahrenheit thermometers.) Converting to centigrade t = 30.5°, ty = 27.2°, 
and D = 3.3°. From the table drawn up by Whalley (1957) on the thermo- 
dynamic properties of heavy water, L, at temperature 27.2° is 45.278 kilo- 
joules per mole (amended value) or 10816.5 calories per mole, and the vapor 
pressure at 27.2° is 0.03125 bar. Therefore 
0.03125 


gi a ES 9 
Hw 1.01325 — 0.03125 Onerrae. 





Hence 
LM = Stat 


and since 6.95D = 23.164 the numerator in the expression for H is 321.049. 
The denominator is 10816.5+8.43.3° or 10844.5; hence H = 0.0296048. 
Converting molal humidity to relative humidity by the formula 


a ia entihnes —_1.01825K100— 
1+H ‘saturated vapor pressure at 30.5° C in bars 





we reach a percentage relative humidity figure of 75.4% from the above 
quoted psychrometer readings. 

For ordinary water vapor a relative humidity of 75% at dry bulb tempera- 
ture 87° F corresponds to a depression of 63° F: so this determined the condi- 
tions at which readings had to be taken when ordinary water was used. 


EXPERIMENTAL METHODS AND RESULTS 


Apart from the matter of providing atmospheres which were reliably com- 
parable, there was the problem of providing suitable samples, because previous 
experience with these gelatin films had indicated that samples did not all act 
alike. If the investigation was to yield results of any value, several samples 
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had to be used. Accordingly 13 small pieces of gelatin film about ? in. X} in. in 
size were painted with conducting rubber paint on either side so as to give an 
area of about in. X% in. through which current could pass. They were mounted 
on a piece of polystyrene carrying 13 pairs of electrodes suitably spaced to 
take the films: a switching arrangement could connect any particular pair to 
the megohmmeter. The gelatin resistors were alternately dried out and humidi- 
fied. The drying was done by putting calcium chloride in the humidifier 
(Niven 1958) with a sensing element for relative humidity 19%. 

Approximate relative humidities were of no use whatever for this investi- 
gation; accuracy in reading the psychrometer had to be pushed to the limit 
but since the sensing elements permitted remarkably good repetition of con- 
ditions, psychrometer readings would repeat to one quarter of a degree 
Fahrenheit in depression values. Before taking the resistance readings the 
megohmmeter was allowed to heat up, but as soon as the leads to the switching 
mechanism were connected, readings were taken just as quickly as they could 
be written down: leaving voltage on the film completely upsets the resistance 
reading possibly for days (Niven 1959). It seemed preferable to run a series 
of tests with heavy water and then change to ordinary water rather than to 
alternate heavy and ordinary tests, because there was then less chance of a 
residue of the wrong kind of hydrogen nucleus being in the gelatin. 

The results are shown in Table I for the 13 different samples: each set of 
results was taken after the film had been dried out and allowed to come 
back to its value in the humid atmosphere. 

The blanks in Table I for No. 13 arose because the voltage was accidentally 
left on, thus vitiating the usefulness of these readings. The ratios calculated 
from the averages show that the replacement of the protons with deuterium 
nuclei does increase the resistance quite markedly. Just why this is not 
approximately the same in every case is not clear: certainly the conditions 
in the humidity cabinet cannot account for this because all of the resistors 
were exposed to the same ambient humidities. 

The average of the ratios in Table I is 1.62; Dr. Wu expected +/2 or 1.41, 
according to the following reasoning. The drift (v) of an ion of mass m and 
charge e in an electric field E is the acceleration 


(2 mean free path — eE d 


m 


thermal velocity om x V (8kT m) 

where A is the mean free path of the ion, 7 is the temperature, and k is Boltz- 
mann’s constant. If there are 7 ions per cc, resistance = E/nev and is therefore 
proportional to +/m. 

Admittedly the disagreement between theory and experiment may be due 
to experimental errors. But it is also conceivably possible that it arises from 
much more fundamental reasons: the idea of collisions due to thermal agitation 
holds good in the theory proposed by Riehl because these collisions are respon- 
sible for making and destroying the chains of water molecules or ‘‘canals”’ 
down which a proton is passed. These ‘‘canals”’ consist of momentarily correctly 
oriented water molecules, and the proton is visualized as moving only from 
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one sphere of influence to the next, and releasing a proton from the water 
molecule to which it arrives and this proton moves on to the next water 
molecule. The idea of mean free path does not fit into such a model too well 
unless it be construed as the length of a “‘canal”’. 

However, the fact that changing the mass of the carrier of positive charge 
did alter the conductivity supports the proton migration theory of conduction 
in hydrophilic films. 
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HALF LIFE, Qs VALUE, AND y-RAY SPECTRUM OF La! 


K. FritzeE, T. J. KENNETT, AND W. V. PRESTWICH? 


ABSTRACT 


The decay of La'* has been investigated. A half life of 14.040.1 min and a 
Qg value of 3.30.1 Mev were found. The y radiation is weak and the spectrum 
is rather complex. 


Information concerning the high-yield fission-product La is quite limited. 
Gest and Edwards (1951) showed by indirect measurement that the lanthanum 
parent of Ce'* has a half life of about 19 min. 


CHEMICAL PROCEDURES 


In order to obtain La" as highly enriched as possible, it was separated from 
fission-product barium rather than directly from irradiated uranium. The 
13-sec Ba! (Nethaway 1959) together with the longer-lived barium isotopes 
was isolated as the chloride from which the lanthanum was separated. Two 
somewhat different methods were used to obtain the lanthanum samples. 


Method 1 

The irradiated uranium solution (1 ml) was added to 2 ml BaCl: solution 
(100 mg Ba) and the BaCl, precipitated immediately by the fast addition of 
15 ml of ice-cold HCl-ether mixture (5+1). The BaCl. was collected on a 
47-mm membrane filter and washed with 5 ml of conc. HCI. (Tracer experi- 
ments with Sr showed that about 15% of it is co-precipitated in the absence 
of holdback carrier and about 10% if 20 mg of strontium is present.) The 
BaCl» was dissolved in 10 ml of hot lanthanum solution (0.5 mg of lantha- 
num) inside the filter funnel. The lanthanum was precipitated immediately 
as the hydroxide and filtered over the same membrane as had been used for 
the BaCls. The La(OH); was ‘‘dry”’ approximately 1 minute after the BaCl. 
precipitation. Therefore, the other lanthanum activities which grow in from 
the 11-min Ba! and the 18-min Ba"! were kept at a minimum. After addition 
of more lanthanum and barium carrier, the La(OH); was reprecipitated and 
dissolved in 10 ml of conc. HNO;. To remove yttrium which had grown in 
from the co-precipitated strontium this solution was extracted three times 
with 10 ml of tributylphosphate. This extraction was found to give a decon- 
tamination factor of about 10‘. The nitric acid phase contained about 50% 
of the purified lanthanum which was eventually converted to Laz(C2O4)3. 
The samples were ready for counting approximately 25 minutes after the 
end of the irradiation. 


Method 2 
In his thesis, Nethaway (1959) showed that the exchange of carrier-free 
radiobarium with pre-precipitated BaSO, is very fast. To reduce strontium 
1Manuscript received January 6, 1961. 
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contamination a similar procedure using BaCl, was developed for the investi- 
gation of La'*, The irradiated uranium solution was added to a suspension of 
BaCl, in 15 ml conc. HCl contained in the funnel of a large membrane filter 
assembly. Immediately following this, a vacuum was applied to the filter 
flask and the “‘dry’’ BaCl: washed with 5 ml of conc. HCI. (Separate experi- 
ments showed that the co-separation of carrier-free strontium is only 0.1% 
in this procedure.) Following this the procedure was identical with that of 
method 1 until after the La(OH); reprecipitation. The La(OH)3 was then 
dissolved in 10 M HCl and passed over a small Dowex 1 column. The lan- 
thanum was finally obtained as the oxalate. Counting could usually be started 
15 minutes after the irradiation. 

Generally samples of 100 or 200 ug of U** in 1 ml of water were irradiated 
for 10 seconds in a flux of 3X10" n/cm? sec using the pneumatic rabbit facility 
of the McMaster Reactor. The BaCl: separation was completed 10 to 12 
seconds after the end of the irradiation. 


HALF LIFE 


The La samples were 6 counted with a 2X2 in. plastic scintillator. The 
output of the detector was fed to a DD2 linear amplifier followed by a single- 
channel pulse-height analyzer. This analyzer was used to select energy regions 
in the 6-ray spectrum. The output was recorded in a 256-channel time analyzer, 
using a channel width of 1 minute for the first 190 minutes. After the data 
were read out, the channel width was changed to 10 minutes in order to obtain 
detailed information on the long-lived background. 

In order to minimize the contribution of other components, the single- 
channel analyzer was used to discriminate against La’? and La"! and also to 
eliminate any contribution from the ingrowing Ce™*. Three separate experi- 
ments employing different energy requirements and chemical procedures were 
conducted to determine the La'*® half life by direct measurement. For runs 
1 and 2 lower levels were set at 1.0 and 1.5 Mev respectively, while for run 3 
a window covering an interval from 2.2 to 3.56 Mev was used. These energy 
requirements, combined with the chemical enrichment, reduced the con- 
tribution of other lanthanum activities to less than 2% of the initial activity. 
The conditions and results for these three measurements are summarized in 
Table I. 

The third half-life run was analyzed using a least-squares program. The 
background associated with this run was determined by observing the long- 


TABLE | 


Summarized results of the La!* half-life measurement 


Initial 





activity Half life 
Run Chemistry Counting (c.p.m.) Analysis (min) 
1 Method 1 B>1.0 Mev 49000 Graphically 14.1 
2 Method 1 B>1.5 Mev 16000* Graphically 14.0 
3 Method 2 2.2>8>3.5 Mev 30000 Least-squares 14.0 
fit 
Average 14.0+0.1 


*Loss of lanthanum during purification. 
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lived tail of the decay curve and analyzing it by means of a least-squares 
program. This tail was resolved into 81-minute, 3.7-hour, and constant com- 
ponents. These components were then subtracted from the original data. 
The contribution of the three background components were 81-minute La™, 
1.3%; 3.7-hour La™!, 0.07%; and room background 0.07% of the initial 
activity. The corrected data were then analyzed for a region covering 10 half 
lives (140 channels). The half life was found to be 14.03 min with a standard 
deviation of 0.03 min. 


COUNTS PER MINUTE 





2 os. 
so 100 200 
TIME (MINUTES) 
Fic. 1. The decay curve of La'® obtained in run 3. In this case only 8 rays whose energy 
fell between 2.2 and 3.5 Mev were recorded. The contribution from other activities amounted 


to 1.5% of the original counting rate. After this small correction was made, a least-squares 
analysis of the data from initial counting until 140 minutes gave a value of 14.03+0.03 min. 


In view of the marked discrepancy between the half life reported here and 
the previously accepted value of 19 min it was decided to perform a similar 
experiment to that of Gest and Edwards. The Ce" was separated from aliquots 
at 6.5, 13.5, 20.5, 27.5, 171.5, and 215.5 minutes after the preparation of a 
stock solution of fission-product lanthanum using the Ce‘+-hexone extraction 
(Glendenin et al. 1955). This method, which is described in more detail by 
Gest and Edwards (1951), gave a half-life value of about 12 min. This value 
is in agreement with the direct measurement. Finally it can be said that the 
15-min rare-earth activity observed by Hahn and Strassmann (1943) un- 
doubtedly must have been La’, 


8B ENERGY AND Q VALUE 
To obtain the 8-ray spectrum of La™ a 2X2 in. plastic scintillator (Nuclear 
Enterprise NE 102) coupled to a 9536B EMI photomultiplier was used. The 
output of this detector was fed to a 256-channel pulse-height analyzer. The 
data stored in the analyzer were transcribed on paper tape in order that com- 
putation for a Fermi analysis of the spectrum might be easily performed. 
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The La" source was prepared using the BaCl, exchange reaction (method 2) 
and 3 mg of lanthanum carrier. The final step in source preparation consisted 
of ashing the oxalate precipitate and mounting the La,Q; on a thin mylar 
film. This source was counted 18 minutes after the irradiation for a 20-minute 
period. 

The spectrometer was calibrated before and after the analysis of the La! 
sample using Sr9°-Y*% (2.26 Mev), Ru!*-Rh! (3.53 Mev), Cl®* (4.9 Mev), 
and the Compton edges of Cs'7 and Zn®. Fermi—Kurie plots were obtained 
and the end points determined. The calibration indicated good linearity of 
the spectrometer. The end point for La was found to correspond to an energy 
of 3.30+0.08 Mev. 

An experiment conducted to determine the average number of y rays asso- 
ciated with the decay of La™* gave a value of less than }. This low y-ray yield 
suggests that-the end point observed in the 8-ray spectrum, 3.3 Mev, probably 
corresponds to the Q value for La’. This value is not inconsistent with the 
decay-energy systematics in this mass region. 


La’ y-RAY SPECTRUM 


The characteristic y radiation from La'* was investigated using scintillation 
spectroscopy. An unshielded 3X3 in. Nal (Tl) detector was employed in 
conjunction with a double delay-line amplifier and a 256-channel analyzer. 
As mentioned above, this particular nuclide was found to be a relatively 
weak y-ray emitter, so that both the contribution from La! and ingrowing 
Ce! was expected to be important. The contribution from La"! is negligible. 
Thus in order to obtain the pure La y-ray spectrum, the spectral contri- 
bution of La'?+Ce'* had to be determined and removed. To reduce the 
effects of 8 radiation a 1 g/cm? polyethylene absorber was placed on the face 
of the crystal. 

A source was prepared using method 2 and the lanthanum oxalate sample 
mounted directly on the polyethylene absorber. The initial counting rate was 
~ 3000 counts/sec and resulted in a dead time of about 4% in the analyzer. 
From previous experience, rates of this order were known to produce no dis- 
tortion in the data. Spectra were obtained beginning at 14 minutes, 162 
minutes, and 18.3 hours after the end of the irradiation, for counting periods 
of 8 minutes, 32 minutes, and 1 hour respectively. The first spectrum was 
composed of La'® (14.0 minutes), La'? (81 minutes), and the Ce'* (33 hours) 
which had partially grown in. The second spectrum, taken after the complete 
decay of La", consisted only of La’? and the Ce! which is fully grown in 
at this time. The third spectrum was solely that of Ce™*®. The spectra were 
read out on tape, allowing reduction of the data using a Bendix G-15 computer. 
From a knowledge of the counting times, and the half lives involved, the 
interfering La'#?+Ce'* components were removed. The contribution of the 
long-lived activities was found to be 20%. The resultant La™* spectrum is 
shown in Fig. 2. 

The energy calibration was obtained both with standard sources and known 
lines in the La'*?+-Ce' spectra as given by Schuman, Turk, and Heath (1958) 
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Fic. 2. The y-ray spectrum of La'’, The spectrum shown has been corrected for the 
contribution from La™? and Ce, The region around 0.300 Mev is uncertain because of the 
strong interference caused by the 0.290-Mev line in Cel. 


and Martin et a/. (1956). Table II lists the observed lines in order of increasing 
energy. These intensities have been corrected for detection and photofraction 
efficiency of the Nal (TI) crystal using the results of Miller, Reynolds, and 
Snow (1957). 

The 0.265- and 0.320-Mev lines shown in Fig. 2 are not given in Table I] 


TABLE II 
y-Ray lines for La! 





Relative Relative 





Energy (Mev) intensity Energy (Mev) intensity 
0.200+0.01 0.08* 1.58+0.04 0.28 
0.440+0.02 0.13 1.70+0.04 0.19 
0.625+0.01 1.0t 1.98+0.03 0.35 
0.800+0.01 0.44 2.22+0.03 0.06t 
0.915+0.02 0.08 2.46+0.03 0.13 
1.07+0.04 0.26 2.56+0.03 0.27 
1.17+0.02 0.57 2.85+0.04 0.15 


*The shape of this line indicates that it is real; however, some contribution from 
backscattering is to be expected. 

tNormalized to unity. 

tThis line may possibly arise from the summing of the 0.625- and the 1.58-Mev 
transitions. 


because their existence is somewhat in doubt. This ambiguity arises from 
interference of the intense 0.290-Mev line in Ce™*. No similar problem arose 
in the subtraction of the La? spectrum. 
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NOTE ADDED IN PROOF 


Recent measurements in this laboratory indicate that the half life of La! 
is 92 min. In view of this discrepancy in the La™ half life, all the data pre- 
sented in this paper have been recalculated. It was found that this change in 
the La! half life did not alter any of the data quoted for La! within the 
assigned errors. 








THE EFFECTIVE DIRECTIONAL SENSITIVITY OF 
COSMIC-RAY NEUTRON MONITORS! 


S. M. LAPoInTE? AND D. C. KosE 


ABSTRACT 


The direction of maximum sensitivity of a neutron monitor is ‘calculated 
numerically for a set of points on the same geomagnetic meridian but extending 
in latitude from the equator to the pole. This leads to two master curves, one 
for the longitude, the other for the latitude of this direction. From these curves 
this direction is obtained in geographic co-ordinates for some 20 cosmic-ray 
stations. The method of calculation is described taking into account atmospheric 
absorption and the energy spectrum of the incident particles. The aperture of 
the sensitive cone, or source width, is also calculated. Finally the accuracy of 
the results is discussed and the application of the concept of effective direction 
is described. 


1, INTRODUCTION 


Neutron monitors, usually consisting of a layer of paraffin surrounding 
neutron counters, then a layer of lead, and then an outer layer of paraffin, 
have been used for several years as detectors of cosmic rays at ground level. 
They are particularly important in studying cosmic-ray intensity variations, 
because atmospheric effects can be eliminated by a simple correction for 
changes in the barometric pressure. 

In problems relating to intensity changes due to solar activity or due to 
the earth’s rotation it is desirable to know the direction of travel of the 
particles being studied at some distance away from the earth, that is, before 
they are deflected by the earth’s magnetic field. It is also important to know 
whether a detector at ground level sensitive locally to particles in all directions, 
and accepting all particles above a minimum energy, can in fact be said to 
be measuring incoming particles from a unique direction in space or at least 
an effective direction that has some meaning in examining intensity variations 
that may be related to directional anisotropies in space. 

The passage of the particles through the atmosphere, interactions with air 
nuclei, and the generation of secondaries which result in a strong collimating 
effect on the particles detected, is reasonably well understood and the behavior 
of charged particles in the earth’s magnetic field has been investigated suffi- 
ciently to give reasonably accurate information on the above problem. In 
this paper we are reporting calculations based on available information which 
will give an effective direction of the primary particles responsible for the 
intensity recorded by a conventional neutron monitor. 

Actually a neutron monitor at sea level responds to primary particles in a 
wide range of primary particle energies, from about 1 Bev upwards. It also 
detects particles impinging on the top of the earth’s atmosphere from a wide 
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range of directions. Primary particles of different energies and different impact 
directions on earth come from different asymptotic directions far away from 
the earth ($10 earth radii) outside the deflecting geomagnetic field. The 
neutron monitor’s sensitivity is different for different primary energies or 
different impact directions on earth. From available information a weight 
can be given to each asymptotic direction, using the known energy spectrum 
of the incident radiation. In this manner an effective direction can be calculated 
for each neutron monitor station on earth. The concept of effective direction 
was introduced by Dorman (1957). It has been applied to the determination 
of onset times of Forbush-type decreases by Fenton, McCracken, Rose, and 
Wilson (1959) and the results of the present calculation have been used by 
Rose and Lapointe (1961) in the analysis of the recovery period after a series 
of Forbush-type decreases. 

In what follows, a method for calculating numerically the effective direction 
of a neutron monitor is described in detail taking various factors into account. 
The calculation was performed for neutron monitor stations at all geomagnetic 
latitudes between the equator and the pole by steps of 5° or 10°. The results 
are presented in the form of a curve. From this master curve the effective 
direction can be obtained easily and rapidly for a neutron monitor in any 
geographic location. 

Finally the accuracy of the results is discussed and the application of the 
concept of effective direction is described. 


2. THE METHOD OF CALCULATION 


The direction of the particle’s trajectory at infinity can be specified by two 
angles: the asymptotic latitude ¢x and the asymptotic longitude yz (Brun- 
berg 1953). The dependence of these angles on energy, impact latitude (geo- 
magnetic), and impact direction, already published, has been used (Jory (1956) 
and Liist (1957) for low energies; Brunberg (1953) and Brunberg and Dattner 
(1953) for high energies). The problem then is to establish an effective weight 
at the detector for each energy and each zenith direction, that is, the zenith 
direction at impact on the atmosphere. 


(a) Zenithal Weights 

The mechanism giving rise to decreasing detector sensitivity with increasing 
zenith angle at the point of impact is atmospheric absorption. It has been 
experimentally established by Simpson, Fonger, and Treiman (1953), and by 
Simpson and Fagot (1953) that the counting rate in a neutron pile monitor 
decreases exponentially with increasing atmospheric depth. The attenuation 
length for the secondary nucleonic cascade in the atmosphere used in the 
present calculation is 150 g/cm*. This leads to a zenithal weight 


volts) 


neglecting the curvature of the earth where h is the thickness of the earth’s 
atmosphere ~1000 g/cm?, L is the attenuation length 150 g/cm?, and Z is 
the zenith angle. 
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The size of the solid angle at a given zenith angle must also be taken into 
account. The method introduced by Brunberg and Dattner (1954) was used 
here. Namely, the sky above the observer is divided first in annular regions 
with boundaries Z = nX8° (n = 0, 1, 2, 3, 4, 5, 6). Then these rings of equal 
widths are cut in four quadrants N, S, E, and W by two appropriate per- 
pendicular planes. In the N quadrant, for instance, there are six regions with 
solid angle proportional to sin Z,, where Z, = (n+3)X8°, (m = 0, 1, 2, 3, 
4, 5). 

The procedure followed to construct a histogram for asymptotic longitude 
is then the following: particles of a single fixed energy are now considered. The 
asymptotic longitudes are taken from the references above for particles of 
that energy having impact directions in the N-S plane and making zenith 
angles of 0°, 8°, 16°, 24°, 32°, 40°, and 48° N. Above 48° the weight is suffi- 
ciently low to neglect the larger zenith angles. These asymptotic longitudes 
are marked on the abscissa of the N histogram. Between the marks for 0° 
and 8° N a rectangular surface is plotted on the histogram with area equal to 


noe Ge) 
sin 4 exp | t ia al 


In other words, asymptotic longitudes for impact directions between zenith 
angle 0° and 8° N are given equal weights totalling to the weight of the region 
centered on 4° N and so on for other zenithal regions. The same thing is done 
for the S, E, and W quadrants. Then all these four histograms are added 
together to give the final histogram for asymptotic longitude for incident 
primary particles of that chosen energy. The same thing is done for asymptotic 
latitude. 


(b) Energy Weights 

The histograms in Section (a) give an answer to the following question: if 
primary cosmic-ray particles of a given energy fall isotropically on a station, 
from what direction do the most efficient particles come? Or, in other words, 
what is the direction of maximum sensitivity of the neutron monitor for 
particles of that energy? 

Now a further question must be answered: what is the direction of maximum 
sensitivity of a given neutron monitor with respect to a mixed beam containing 
particles of many energies? Obviously, many histograms of the type described 
in Section (a) must be obtained for many different energies; each one must 
be multiplied by an appropriate energy weight representing the abundance 
of particles of that energy in the incident beam. Then all the resulting histo- 
grams added together give the final one. 

However, the energy weight is not given simply by the relative abundance 
of the various energies or, in other words, the differential energy spectrum 
of the incident beam. It must be remembered that one is dealing here with a 
detector at sea level and that one is looking for the direction in space which 
gives rise to largest counting rate in such a detector. Therefore, the energy 
weight is simply the primary differential spectrum times the specific yield 
function (Simpson, Fonger, and Treiman 1953). 
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Equivalently one can use the formalism of the coupling constant (Dorman 
1957) and for practical reasons this will be done in what follows. Following 
Dorman (1957) one can write the following equations. The total counting rate 
at geomagnetic latitude \ and atmospheric depth ho is: 


®%m~ 


Ny(ho) = J D(e)m(e, ho)de 
€,min 
where emin is the geomagnetic threshold energy, D(e) is the differential energy 
spectrum of the primary cosmic rays, and m(e, ho) is the multiplicity or specific 
yield function at depth Ao. Considering fixed atmospheric and geomagnetic 
conditions the intensity variation is given by: 


®*a 


6Ny (ho) = J 6D(e)m(e, ho)de 
€,min 
and obviously the relative variation is: 
iNalhe) _ f° BDO meh g, 
Ny (ho) €\min D(e) Ny (ho) 


%m» 6D 
i . Ry Wyle, ho)de. 


The last equation defines the coupling constant W,(e, io). In this formalism, 
the effective angle described above can obviously be defined: 





J or (€) Wy (e, ho) 6D (e) ‘ee 
dr = — QO De) 
%@ : sD(e) 
a W ale, ho) D(e) de 


This represents the effective angle of a given variation in the intensity of the 
primary particles 6D(e)/D(e). The present calculation is carried out for an 
intensity variation having the same spectrum as the ordinary cosmic rays, 
that is, a power law spectrum; this is the case, for instance, of the Forbush 
decrease. Therefore, one can write 6D(e) = CD(e) where C is a constant, and 


J — gale)Wy(e, ho)de 


ie) €,min 
8 eri sere nici netstat 


| . Weleda 


€,min 


In the present paper, ¢(e) is not a single angle but a whole histogram. However, 
this makes no difference to the above expression when it is integrated numeri- 
cally in the manner described below. A point of practical interest is the following: 
an examination of Dorman’s curves for the coupling constant for neutron 
monitors shows that, though it is a function of latitude, its variation with 
energy is the same for different Jatitudes. To a close approximation in this 
application one can write: 
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Wy(e, ho) 
Weo (€, ho) 


where C,(ho) varies only with latitude. Therefore, 


= Cy (ho) 


J . __ 6.(€) Wat (6, he) de 
Dr a €,\min “ 


J Wso° (e, ho)de 
€\min 


The coupling constant for the neutron component at sea level used by Dorman 
was derived from the measurements of Rose and Katzman (1956). To perform 
the above integration the energy axis is divided in intervals of various lengths. 
The interval is chosen of such a length that the asymptotic angle over that 
length does not vary too much. Then the energy in the middle of each interval 
is chosen to work out a zenithal histogram; that histogram is then multiplied 
by the following energy weight: length of the energy interval times the value 
of the coupling constant at the middle energy. This is done for energies 
extending from a practical cut-off, where particles from more than half the 
directions are admitted, up to energies where an additional interval adds a 
negligible weight to those already accumulated. All the resulting weighted 
histograms are added together to give the final one. 

From this final histogram the mean effective angle is calculated, and also 
the root mean square deviation giving a measure of the width of the sensitivity 
aperture of the neutron monitor. 

3. RESULTS 

The above described calculation has been performed for geomagnetic lati- 
tudes: 0°, 10°, 20°, 25°, 30°, 35°, 40°, 45°, 50°, 55°, 60°, 65°, 70°, 75°, 80°, 
90°. For 25°, 35°, and 75°, interpolation of the raw data had to be used. The 
results are given in the form of two curves. Figure 1 presents the longitude 
shift of the source to the east of the observer’s meridian Wy, as a function of 
the geomagnetic latitude of his station. Figure 2 shows the northern latitude 
of the source ¢x. Both asymptotic angles yg and @y are expressed in geo- 
magnetic co-ordinates. To find the asymptotic effective angle at a given cosmic- 
ray station one proceeds as follows: the longitude shift Y_ read off the master 
curve is added to the station’s geomagnetic longitude ¢,. This, together with 
gn read off the second master curve, gives the source’s geomagnetic co-ordi- 
nates. Using the Smithsonian Physical Tables one can convert these to 
geographic co-ordinates. This has been done for some 20 cosmic-ray stations 
and the result is shown in Table I. The present calculation agrees within 
1 or 2 degrees with the four longitude shifts calculated independently by 
Fenton, McCracken, Rose, and Wilson (1959). 

The root mean square deviation o of the final histogram of the process is 
plotted in Fig. 3 for the longitude shift yz, and in Fig. 4 for the latitude ¢y, 
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Fic. 1. East longitude Yr of the effective direction of the monitor vs. station latitude, 
both in geomagnetic co-ordinates. Station longitude taken to be 0°. 
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Fic. 2. Latitude ¢y of the effective direction of the monitor vs. station latitude, both in 
geomagnetic co-ordinates. 
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TABLE I 


The first two columns give the station’s location, the last two columns 
give the direction of maximum sensitivity of a neutron monitor there 


Geographic co-ordinate 
Geographic co-ordinate of direction of maximum 








of station sensitivity 

Station Lat. Long. Lat. Long. 
Sydney —34 151 E 0 230 E 
Berkeley 38 122 W —16 300 E 
Climax 39 106 W —14 310 E 
Mt. Wellington —43 147 E —4 205 E 
Chicago 42 88 W 0 315 E 
Ottawa 45 76 W 3 320 E 
Deep River 46 78 W 6 320 E 
Sulphur Mt. 51 116 W 5 272 E 
College 65 148 W 28 228 E 
Churchill 59 94 W 31 285 E 
Uppsala 60 18 E 27 72E 
Mawson —68 63 E —44 TOE 
Zugspitze 47 l1E 7 85 E 
Munich 48 1ZE 8 82E 
Herstmonceux 51 0 14 60 E 
Leeds 54 2W 20 53 E 
Ellsworth —78 41W —4d4 343 E 
Mt. Washington 44 71W 2 328 E 
Rome 42 IZ E 6 95 E 
Pic du Midi 43 0 5 78 E 
Mt. Norikura 36 138 E -—5 209 E 
Lincoln 41 97 W —10 310 E 
Resolute 75 96 W 64 263 E 





Thule 77 69 W 74-84 


both as a function of station geomagnetic latitude A. This deviation o gives 
an idea both of the uncertainty of the calculated mean effective angles and 
also of the width of the sensitive aperture of the neutron monitor. 


4. ACCURACY 


Unfortunately, it is not possible to estimate accurately the uncertainty in 
these calculations. Therefore, some factors affecting the uncertainty will be 
discussed here in an attempt to give an idea of the magnitude of the possible 
errors. 

(a) The accuracy of the fundamental data used, namely, the asymptotic 
directions of cosmic-ray trajectories as a function of energy and geomagnetic 
latitude, are estimated by Brunberg (1953) not to exceed 2° to 3°. On the other 
hand the agreement between the experimental points of Brunberg and Dattner 
and the calculated ones of Jory and Liist is also within about the same limits 
of a few degrees. Therefore, in the present calculation, all the angles were 
read to the nearest 5°: histograms were, therefore, obtained for 5° intervals. 
This is likely to be as good as the basic data. In averaging over a range of 
angles wide compared to 5°, the inaccuracy on the resulting average angle is 
probably less than the one on individual points. The smoothness (within 1°) 
of the curves in Figs. 1 and 2 seems to corroborate this argument. 
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Fic. 3. Source width in longitude vs. station latitude, both in geomagnetic co-ordinates. 
Fic. 4. Source width in latitude vs. station latitude, both in geomagnetic co-ordinates. 


(6) The shape of the curve at low latitudes in Fig. 1 where there is a double 
reverse in curvature deserves special comment. The variations in Wx, in this 
region amount to between 3° and 7°. All the neighboring points in that region 
have been treated in the same fashion, and therefore, if any, they should all 
bear the same bias and not vary up and down. Perhaps, therefore, this fluctua- 
tion is significant. The only test that can be applied in the present method is 
to break up an energy interval in a critical region (usually low energy where 
the asymptotic angles change rapidly) into two or three narrower intervals 
and with this to calculate the new mean angle in such a finer grid. This was 
tried in the low latitude region and it is felt that a further subdivision would 
not bring about a 1° change. If these variations are not significant the reason 
should probably be sought elsewhere, such as in the accuracy of the basic 
data, cut-off values, the approximation in the use of a single attenuation 
length, or perhaps other factors more critical in this region of higher magnetic 
field. It should be noted that the r.m.s. deviation plotted in Fig. 3 has its 
largest value at low latitudes, therefore, indicating greatest uncertainty in 
that region. This might explain the irregular curve. However, the r.m.s. 
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deviation curve has itself the same irregularity in that region, indicating that 
the width of the sensitive cone fluctuates. This fluctuation in source width 
is quite real in regard to these calculations but its physical significance is 
not clear. 

(c) Another source of error is the attenuation length L used in the exponential 
term of the azimuthal weight. The same value of 150 g/cm? was used for all 
stations. However, it is well known that this length varies with altitude and 
latitude. This constant value can only be considered a rough and practical 
average value. 

5. DISCUSSION 

In conclusion, perhaps the more exact meaning of effective direction can 
be reviewed in the light of its derivation above. The method of calculation 
described here assumes an isotropic beam (at least within a minimum cone 
around the station’s vertical) falling on the top of the atmosphere and having 
the same differential energy spectrum as ordinary cosmic rays, that is, an 
inverse power law. Then effective direction is that direction from which the 
particles are coming which are preponderant in producing the counting rate. 
It is obvious that particle trajectories from other directions could be absent 
without affecting the counting rate very much. In other words, if the beam 
is isotropic, by and large the particles counted are those coming from the 
effective direction. If an anisotropic variation comes from the effective direction 
it will be recorded with maximum sensitivity. But the reverse to this last 
statement is only probably true: that is, a recorded intensity variation only 
probably comes from the effective direction. Given any recorded counting 
rate variation it can be caused by a primary variation taking place in the 
effective direction, or by a larger one not in the effective direction. But if 
the intensity variation recorded is of short duration, for instance a peak above 
the neighboring background in the intensity-time curve, and if it occurs at 
many stations, and above all if all these stations’ peaks and effective directions 
point at the same direction in space for the source, then it is very probably 
meaningful. 
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MODELS OF AURORAL IONIZATION 


PART I. AURORAL IONIZATION MODELS 
AND THEIR RADIO-REFLECTION CHARACTERISTICS! 


D. R. Moorcrort 


ABSTRACT 

Although radio observations of aurora contain information about the nature 
of the reflecting ionization, the use of dissimilar models of auroral ionization 
has led different workers to widely differing conclusions. In. this paper several 
general models of auroral ionization are developed. By considering the ionization 
as an assembly of individual scatterers, it has been possible to include a unified 
treatment of both weak scattering and critical reflection. This treatment should 
provide a basis for resolving some of the difficulties in the interpretation of 
auroral radio observations. In a second paper the available experimental 
evidence is examined in the light of this theoretical treatment. 


INTRODUCTION 

Over the past thirteen years radar has been used extensively to study auroral 
ionization, but the distribution and density of the ionization responsible for 
these radar reflections is still the subject of much speculation. Perhaps the 
most fundamental question to be resolved is that of the density of this ioni- 
zation. In terms of radar reflection this is a question of the relative importance 
of two basically different reflection mechanisms: critical reflection and weak 
scattering. The first of these mechanisms (sometimes known as ionospheric- 
type reflection) was proposed originally by Lovell, Clegg, and Ellyett (1947) 
to explain echoes observed by them on 46 Mc/s. More recently, other workers 
(Forsyth 1953; Bullough and Kaiser 1954; Kaiser 1956; Forsyth and Vogan 
1957; Collins and Forsyth 1959) have suggested that critical reflection is 
important in explaining auroral echoes, at least on lower frequencies. Some 
of the workers (Herlofson 1947; Aspinall and Hawkins 1950; Booker, Gartlein, 
and Nichols 1955; Booker 1956; Presnell et al. 1959; Flood 1960) believe 
that all echoes are due to irregularities of small strength (i.e. electron density) 
in the normal ionosphere, the mechanism of reflection being partial reflection 
from these irregularities. Also in doubt are the size of the irregularities and 
the degree of elongation along the magnetic field. Estimates of the radius 
of the irregularities range from a fraction of a meter (Booker 1956) to 100 
meters (Kaiser 1956); estimates of their elongation range from an axial ratio 
of 4:1 (Forsyth and Vogan 1957) to 40:1 (Booker 1956). Part of the reason 
for this controversy is the lack of a sufficiently general model with which the 
experimental results may be compared. The most developed theory of auroral 
reflections (Booker 1956) considers only weak scattering. While many authors 
have discussed the general characteristics of critical reflection, there is no 
detailed theory for the reflection of radio waves from critically reflecting 
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irregularities, nor of the transition from critical reflection to weak scattering 
that will take place with increasing radio frequency. A satisfactory theory 
should provide a thorough and unambiguous basis from which to examine 
the existing experimental data, and should suggest still other types of experi- 
ments needed to prove the existence of particular models of auroral ionization. 

The purpose of this paper is to develop a general theory which encompasses 
both critical reflection and weak scattering, and includes several different 
distributions in the size of the scatterers. The results of the theoretical develop- 
ment are examined in Part II (Moorcroft 1961) in the light of the available 
observational evidence. 

In order to include several different scattering mechanisms, a rather versatile 
model is required. One quite useful method has been used by Booker (1956) 
in the development of his model of auroral ionization. He described the dis- 
tribution of the ionization in terms of an autocorrelation function associated 
with variations in the dielectric constant (or electron density), and derived 
the scattered field in terms of the parameters of this autocorrelation function. 
This approach excludes consideration of critical reflections, and is therefore 
not satisfactory for a general treatment of reflection processes. In this paper, 
a family of models, consisting of assemblies of individual, specified scatterers 
is considered. The scattered field from the assembly of scatterers is obtained in 
terms of the scattered field from each of the individual scatterers. This permits 
the inclusion of critical reflections and allows the model to be more directly 
related to the physical situations which are believed to exist at auroral levels. 
During the development of the models a considerable number of approxi- 
mations are introduced. These approximations can be justified because they 
permit theoretical consideration of otherwise intractable problems, and the 
results based on the approximations are still of greater precision than the 
observational material with which they are to be compared. 

All the proposed models are assumed to consist of assemblies of scatterers 
having an electron density distribution 


(1) 
wherea < ¢, &, n, and ¢ are spatial co-ordinates in meters, and N is the electron 
number density in m~*, These may be called prolate-spheroidal, Gaussian 
distributions; for convenience they will be called ‘‘gaussoids” in this paper. 
This distribution can arise by diffusion from a point concentration of ioniza- 
tion in the presence of a magnetic field; this process may occur at auroral 
heights in the ionosphere. 


- —[(E249)2 /a2-+(22 se2 
N = Me ((€2+)? /a2+ (2 /c?)] 


’ 


THE SCATTERED FIELD FROM A GAUSSOID 
The generalized behavior of a gaussoid in the presence of an electromagnetic 
wave is very difficult to determine. Several limiting situations can be con- 
sidered, for each of which a different, special technique is available. This 
approach is analogous to the very successful treatments of radio-wave reflection 
from meteoric ionization (Eshleman 1952; Kaiser and Closs 1952). The 
limiting situations, amenable to theoretical treatment are weak scattering, 

critical reflection (ray theory), and small, strong scatterers. 
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Weak Scattering 

If the electron density is sufficiently small or the frequency of the electro- 
magnetic wave sufficiently high, the wave is essentially unattenuated by the 
ionization, and the Born approximation is valid. This implies that the electrons 
are far enough apart to be unaffected significantly by reradiated fields. Only 
the effect of the field of the incident wave needs to be considered. The scattered 
field from such an assembly of electrons is given by the sum of the fields 
reradiated by each electron separately, taking into account their relative 
phases. When these approximations are valid the reflection is said to be due 
to weak scattering and the resulting treatment is called scatter theory. If 
N(é, n, £) is the electron density in m~*, (/,, m1, 1,) the direction cosines 
specifying the direction of the incident wave normal, (/2, m2, m2) the direction 
cosines specifying the direction to the field point, k=2z/\ the wave number 
in m=, E, the incident field, 7) the distance to the field point from the origin 
of the co-ordinate system (é, 7, ¢), assumed to be within the scattering volume, 
a the angle between E, and fo, and r.=e?/4am the classical electron radius, 
the scattered field E, at the field point is 


3 resinak; x 86 ade E tee se 
(2) E, = a ettto Sf xe ” rye Btn 12) &+(m1i—m2) 9+(n1 ™ld¢dnde. 
( 


It is assumed that the scattering volume has dimensions small compared to 
ro. Substituting (1) in (2) gives for the scattered field from a gaussoid 


+ ~ 3/2 
(3) E. = re Sin OE ym CN get? / Pla tt 12) 2402 (mma) 24+ 2 (m—na)?) 


ro 
A scattering coefficient a, defined as the power scattered per unit solid angle, 
per unit incident power density, is given by 
(4) o = r'(E;/E)}). 
The total number of electrons in a gaussoid is given by 
(5) q = Na’cr®” = Nahr”. 
It follows that the scattering coefficient can be expressed by means of (3) as 


y— (a?k? /2)[(li— 12) 2+ (m1 —m2)*+h? (ni—n2)?) , 
’ 


, 198 
(6) o¢=rgq sin ae 
where h=c/a is the axial ratio of the distribution. For the case of back- 
scattering 
lo 2 2 —2a°k?(cos*¢+h? sin?) 

(7) o=rqe a ti alia! 
where ¢ is the angle that the z-axis of the distribution makes with the incident 
wave front. This is the case of particular interest for observations of radar 


reflections from aurora. 


Critical Reflections (Ray Theory) 

If the electron density No increases, or the frequency of the radio waves 
decreases, until the index of refraction of the medium becomes negative, the 
incident wave no longer passes through the scatterer. In an ionized medium 
the index of refraction is related to the electron density by the expression 
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(8) uw = 1—(4ar,/k’)N. 
The index of refraction becomes zero when 
(9) k” = 4ar.No. 


If effects due to collisions are ignored, it can be shown that at normal incidence 
effectively all of the energy in the wave is reflected at the surface defined by 
(9). This is called critical reflection, and the surface will be referred to as the 
surface of critical reflection (N, is called the critical density). If this critically 
reflecting surface has dimensions large compared with the wavelength, the 
propagation of a wave can be described with reasonable accuracy by geo- 
metrical optics, or ray theory. 

For simplicity only the case of backscatter is considered. As a first approxi- 
mation the scatterer can be replaced by a metallic surface coincident with 
the critically reflecting surface: the rays are then all straight lines specularly 
reflected from this surface. The scattering coefficient from a smooth, curved 
surface of this sort is obtained using similar considerations to those of Aspinall 
and Hawkins (1950); 


ar sae Rero ) 
(10) ae A( Rite R2+1ro/’ 


where 7% is the distance from the point of reflection to the field point (receiver), 
and R, and R: are the principal radii of curvature of the surface at the point 
of reflection. Using the results of differential geometry, it can be shown that 
the scattering coefficient of a prolate spheroid which is small compared to 
ro is 





si raceiainsos 
(11) re 4(cos’ @+h° sin’ )*’ 
where p is the semiminor axis of the spheroid, h is the axial ratio, and ¢ is 
the angle that the axis of the spheroid makes with the incident wave front. 
For a gaussoid, p=a[In(No/N,)|* and the scattering coefficient becomes 
a’h” In(No/Ne) 








(12) es 4(cos’ +h’ sin’ @)* ° 

Actually the rays are not straight lines as assumed in this derivation, but 
are curved due to the effect of the electrons outside the critically reflecting 
surface. The computation of the actual ray paths in a gaussoid is quite difficult. 
The effect of refraction can be approximated by using the results of a similar 
problem relating to a cylindrical, Gaussian distribution treated by Manning 
(1953) in connection with overdense meteor trails. He showed that the electrons 
outside the critical surface made the distribution behave as a specularly 
reflecting surface having an effective radius, pe, that can be expressed in terms 
of r,, the radius of the true critically reflecting surface: 


(13) Pe = ref(Ne No). 
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The function f(N./No) is shown in Fig. 1 (after Manning). The scattering 
coefficient for a gaussoid can be approximately corrected for this effect by 






° o-2 0-4 0-6 0-8 0 


Fic. 1. Ratio of the effective radius to the critical radius for a cylindrical, Gaussian dis- 
tribution (after Manning 1953). 


multiplying each of the principal radii of curvature by f(N,/No). Thus 


(14) a‘h'f'(Ne /No)ln(No/Ne) 
<™ 4(cos +h’ sin’ ¢)” * 


The electrons outside the critically reflecting surface also have an effect 
on the angular dependence of the scattered field. Equation (14) is expressed 
in terms of the radii of curvature of the critically reflecting surface at the 
point where the normal to the surface makes an angle ¢ with the plane normal 
to the axis of the gaussoid. The ray which is normal to the critically reflecting 
surface at this point is, in general, curved, and the angle, ¢’, that it makes 
with the plane normal to the axis of the distribution at a large distance from 
the gaussoid is less than ¢. The scattering coefficient can be expressed in 
terms of the angle ¢’ (which will provide the required correction to the angular 
dependence) by replacing ¢ in (14) by ¢’+6(¢’, h, N./No), where 6=o—¢’. 
The definition of ¢’ now corresponds with the original definition of ¢, so the 
prime can be dropped to give 


a*h’f?(Ne/No)ln(Ne/No) _ 
4[cos’(¢+8)-+h’ sin?($+6)]° ” 


Without detailed ray tracing, all that can be said about 6 is that it is zero 
when ¢ is 0° or 90°. 





(15) c= 


Small, Strong Scatterers 

If the gaussoids are small compared with the incident wavelengths, and 
the electron density is large, interactions between the electrons can modify 
the scattered field, giving rise, on occasion, to strong electrostatic fields and 
to resonances for certain polarizations. Even though the peak electron density, 
No, is greater than the critical density, the wave may still penetrate the 
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scatterer if it is too thin to attenuate significantly the evanescent wave which 
exists within the surface of critical reflection. 

If the wave penetrates the scatterer, it acts on the electrons with a force 
which is nearly equal across the scatterer. The resulting electrostatic fields 
can increase the amplitude of the reflected wave to a value greater than that 
obtained from scatter theory. For a Gaussian, cylindrical distribution (Kaiser 
and Closs 1952; Herlofson 1951) the scattered field may be as much as twice 
the value predicted by scatter theory for polarization transverse to the cylinder 
axis, but no increase occurs for parallel polarization. A similar study for a 
spherical distribution (Kaiser and Closs 1952) gives essentially the same value 
of the field as obtained from scatter theory. Although no calculations are 
available in the literature for a gaussoid, it is likely that the amplitude of the 
scattered field will be less than twice the value given by scatter theory and 
probably close to it. Since this study is concerned only with the general char- 
acteristics of the amplitude of the scattered field and not with polarization 
effects, it appears sufficient to use scatter theory to describe reflection from 
small, strong scatterers which are penetrated by the wave. 

If the electron density is sufficiently high to prevent the wave from pene- 
trating the scatterer, the electrostatic forces produced by the motion of the 
electrons are now large enough to keep the amplitude of this motion very 
small; in other words, the scatterer becomes highly conducting. An approxi- 
mate indication of the behavior of a gaussoid in this situation can be obtained 
by considering the behavior of a conducting spheroid. The particular size 
of spheroid that should best represent the gaussoid is not evident. A spheroid 
having the dimensions of the surface k = —1 will be used. This is similar to, 
although more simple than, the criterion chosen by Kaiser and Closs (1952), 
in their discussion of meteor reflections. Assuming that this surface has 
dimensions that are small compared to the wavelength, and that the trans- 
mitted and received polarizations are identical, the following expression is 
obtained for the scattering coefficient: 


(16) o = k’a‘h'In®(No/2N,)[g, sin’ 6+(g, sin’o+g, cos’ ¢)cos’ 6]* 

Here ¢ is the angle that the axis of the distribution makes with the plane of 
the wave front, and @ is the angle that the incident field vector makes with 
the plane which contains the normal to the wave front and the axis of the 
distribution. The functions g, and g, are given approximately in Fig. 2. 
These are the principal dipole moments of spheroids (parallel and transverse, 
respectively, to the axis of the spheroid) of different axial ratios, normalized 
to the dipole moment of a conducting sphere. 


The Regions of Validity of the Different Mechanisms 

The reflection of radio waves from gaussoids for various limiting conditions 
has been discussed. It is now necessary to examine the ranges of size and 
frequency for which these different reflection mechanisms apply. There are 
three conditions separating the mechanisms. The first is the one separating 
critical from non-critical reflection, 


(17) No=N.. 
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Fic. 2. Approximate, normalized, principal dipole moments of a spheroid as a function 
of its axial ratio, h. 


The second distinguishes between large and small scatterers with greater 
than critical electron densities. This can be written 


(18) kr,=1, 


where r, is the radius of the critically reflecting surface. The third condition 
distinguishes between small, strong scatterers which are and are not fully 
penetrated by the wave. This is given approximately by equating the depth 
of penetration of the evanescent wave to the size of the scatterer. Following 
Kaiser and Closs (1952) the depth of penetration is assumed to be equal to 
the skin depth for a plane dielectric. This is given by 


(19) d=1/kV/—x. 


The skin depth is the depth at which the wave amplitude falls to 1/e of its 
original value. Assuming that the gaussoid behaves as a uniform slab of 
thickness ‘a’, and electron density No, 


2 





and the condition becomes 


(20) ka=1l, 
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where k,=(4mr,N,)? is the critical wave number for a scatterer. By defining 
normalized wave numbers and sizes we can represent all the cases on a single 
diagram: 


(21) k'=k/ke, a’=k,a. 
Using these new variables, (17), (18), and (20) become 
k’=1, 
(22) (k’a’)'In(1/k”) =1, 
a’=1 


The resulting diagram (Fig. 3) shows the regions of validity for all gaussoids 
having a given peak electron density (which fixes k,, and thus the scale of 
the diagram). This diagram indicates three different types of frequency 
dependence, depending on the size of the scatterer relative to the wavelength 
of the critical frequency for the scatterer. When ak,<1, scatter theory is 
valid throughout the entire frequency range (keeping in mind the possibility 
of a slight resonance effect at the frequencies below the critical frequency); 
when 1 < a’ < 1.65, the small conductor theory becomes important at the 
lower frequencies; when a’ > 1.65, there is a transition from scatter theory 


Ray Theory 
(kr >!) 


\ Transition 
Small ‘ Conductor Theory : Region 


ta (Nee) mk =-| 


Scatter Theory 
(a < skin depth) 





0 0-2 0-4 0-6 0-8 1-0 2 
k' Critical 
Frequency 


Fic. 3. Regions of validity of the different reflection mechanisms (a’ =ak, is the normalized 
scatterer “‘size’’", and k’=k/k, is the normalized wave number of the radio wave). 
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to ray theory near the critical frequency, and again from ray theory to small 
conductor theory at a lower frequency (just how low depends upon the magni- 
tude of a’). 

Not all of these dividing lines can be given equal importance. The k’=1 
line almost certainly separates scatter theory from ray optics, particularly for 
large values of a’. The boundary a’=1, corresponding to the condition that 
a be equal to the skin depth, is very approximate and indicates only that a 
transition of this sort takes place in that region of the diagram. Similarly, 
the line separating ray theory from small conductor theory could probably 
be represented just as well by the curve kr,=}, which, as can be seen in the 
diagram, is widely separated from kr,= 1, indicating a broad transition between 
these two limiting conditions. 


The Frequency Dependence of a Gaussoid 
An examination of the expressions (7), (15), and (16) for the scattering 
coefficients for the different mechanisms reveals that they can be expressed 
in terms of the normalized wave numbers and sizes as given by (21). With 
the introduction of a normalized scattering coefficient, o’ =k?0/h’, (7), (15), 
and (16) become, respectively, 
ra’ —2k’"a’” (cos’-+h2sin2¢) 


2 —— 
(23) os 16 \ ’ 





5 a’’f’(k’”)In(1/k’”) 
(24) oh = a Si ae Cee Dy 
4[cos’ (¢+6)+h° sin“(¢+6)] 

(25) of = ka"In*(1/2k’”)[g, sin’ 0+(g, sin? +g) cos’ ¢)cos’ 6)’. 

These expressions were used to plot the curves in Fig. 4, which show the fre- 
quency dependence of the normalized scattering coefficient for ¢=0, and for 
an axial ratio of 5. The curves would be the same for all but the small conductor 
theory for any other value of h. For the small conductor theory the average 
of the scattering coefficients for the two orthogonal polarizations, @=0° and 
90°, has been plotted. From Fig. 3, one would expect that below a’=1, the 
scatter theory would be valid throughout the entire frequency range. Examples 
of this type of frequency dependence are shown for a’=0.5 and a’=1. When 
1 <a’ < 1.65, there should be, on the basis of Fig. 3, a transition from scatter 
theory to small conductor theory. The curve for a’ =1.5 illustrates this type 
of behavior. The curves for a’=2, 3, show al] three mechanisms becoming 
important in different frequency ranges. The scattering coefficient for the 
scatter theory decreases very rapidly for values of a’ larger than 2. 

Also of interest is the angular dependence of the scattering coefficient. 
As a measure of the angular dependence or aspect sensitivity (as this charac- 
teristic of radio-auroral reflections has come to be called) a normalized scatter- 
ing coefficient, r=0(¢)/a(0), can be defined. The different mechanisms yield 
the following expressions for this aspect sensitivity factor: 


oR —2a’ 2k’ 2(h2—1) sin? 
(26) Tr =e : —. (scatter theory), 


(cos’[¢+6(¢, h, k’)| +h’ sin*[¢+5(¢, h, k’?)])~?, (ray theory), 


ll 


(27) tr 
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Normalized Scattering Coefficient 


Fic. 4. Normalized scattering coefficient, o’ =k2c/h?, as a function of frequency, for a 


gaussoid (h=5, ¢=0). : F - 
Fic. 5. Aspect sensitivity factor as a function of frequency, for a gaussoid (h=5, ¢6=10°). 
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r,=l1, 
: (small conductor theory). 
28) y= [1-1 g/g sin’, 


All of the variation of rg with frequency is due to the effect of 6(¢, h, k’’), 
the behavior of which is not completely known. It is likely that 7, should 
have some similar modification, but just how this is to be done is not clear. 
Figure 5 gives an indication of the probable type of frequency dependence 
of r for h=5 and ¢=10°. 


ASSEMBLIES OF GAUSSOIDS 


The next step in developing a model of auroral ionization is to consider 
the nature of the field scattered from an assembly of gaussoids. If, as seems 
likely, the scatterers are distributed randomly in space, the power scattered 
from each individual scatterer will add. This is equivalent to adding the 
scattering coefficients. The use of such a procedure assumes that the scatterers 
are all at distances from the field point that are large compared to the dimen- 
sions of the region enclosing them, that they are sufficiently numerous, and 
that the scatterers are neither sufficiently strong nor sufficiently closely packed 
so that multiple scattering becomes important. It will be assumed in this 
section that such a situation prevails. 

The simplest assembly of gaussoids is obviously one in which all of the 
scatterers are identical. For such an assembly a scattering coefficient per 
unit volume, o,, can be defined, which is related to that for the individual 
scatterers by u, the number of scatterers per unit volume: 


(29) o,=Ne. 


Evidently the characteristics of the assembly are essentially the same as 
those of the individual scatterers. It is undoubtedly naive to use such an 
assembly as a model of auroral ionization, yet it is interesting to note that it 
includes as special cases most of the previously suggested models. 

It is of particular interest to compare such an assembly of very weak 
scatterers with the model developed by Booker (1956). From an autocorrelation 
function for the electron density of the form 


—4h[(¢2 )2 fa2+ (£2 /e2)] 
(30) AS, rr 


Booker derived a scattering coefficient per unit volume which can be written 


(31 ) o, = r'(AN)-a‘h (2 Sei Sealine da tama ; 


where (AN)? is the mean square deviation of the electron density. Equations 
(7) and (29) give for an assembly of weak gaussoids, 


2 2 —2%2q2 26+h? sin? 
(32) =e eee. 


The mean square deviation for the assembly can be easily determined if it is 
assumed that the scatterers are sufficiently far apart so that in the volume 
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1/n around one scatterer the electrons due to the other scatterers are negligible. 
Then for the distribution of (1), the mean square deviation is 


(33) (AN)? = nq’ /a*h(24)*”. 


Substitution of (33) in (32) gives (31), which shows that the two approaches 
are essentially similar. The need to require that the scatterers be widely 
separated is clear, since it is only then that the autocorrelation of the medium 
will be that given in (30). In spite of the fact that the two approaches give 
the same expression for the scattering coefficient, there are interesting and 
useful differences in the two results. Since the Booker expression is derived 
on the basis of an autocorrelation function, it describes the scattered field 
from any distribution of ionization having an autocorrelation function of the 
form given by (30), of which media there are an infinite number. However, 
the problem of specifying these media is difficult. On the other hand, the 
expression derived from an assembly of similar gaussoids, (32), is known to 
describe the field scattered from a certain family of scattering media; namely, 
assemblies of similar, weak gaussoids. From (33), any assembly of similar, 
weak gaussoids which leaves the product nq? constant has the same scattering 
coefficient per unit volume. 

A model such as that proposed by Kaiser (1956), which consists of an 
assembly of large, similar, critically reflecting, spheroidal scatterers, is 
adequately described by an assembly of similar, critically reflecting gaussoids: 
in addition the gaussoidal assembly includes the effect of diffuse boundaries, 
and is therefore perhaps more reasonable from the physical point of view. 

This type of assembly also provides a quantitative basis for considering 
a model such as that suggested by Forsyth and Vogan (1957), which involves 
scatterers critically reflecting at one frequency, and weakly scattering at a 
higher frequency. 

Although an assembly of similar scatterers is useful, it is likely that other 
classes of scattering assemblies will resemble the physical situation more 
closely. It therefore seems profitable to consider the physical processes that 
are likely to control the distribution of ionization in the aurora and the sort 
of assemblies that they suggest. 

The irregularities in auroral ionization may be formed in several different 
ways, but it seems certain that the sharpness of the boundaries of such irregu- 
larities will be limited by diffusion. The earth’s magnetic field will cause 
electrons moving across the field lines to gyrate around them; the diffusion 
will therefore be slower across the field lines than along them, although this 
anisotropy will disappear at heights where the collision frequency becomes 
comparable with the gyro-frequency. A rigorous examination of this problem 
requires consideration of the ambipolar effects in the diffusion. Apparently 
such an examination has not been carried out, but Kaiser (1953) has given 
a simplified treatment which suggests that the magnetic field probably intro- 
duces anisotropy into the diffusion at heights near 100 km and above, the 
anisotropy increasing with height. The distribution of ionization arising due 
to diffusion from a point in the presence of a magnetic field is a gaussoid. 
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Although many pictures of auroral ionization would directly imply diffusion 
from a point, and hence gaussoidal scatterers, it seems likely that even irregu- 
larities produced by other processes will, under the action of diffusion, acquire 
characteristics quite accurately represented by gaussoids. 

It is possible to think of several basically different physical processes con- 
trolling the formation of auroral ionization, each of which leads to a different 
type of scattering assembly. The ionization might be produced in small 
clumps or columns by high-velocity, incoming particles continually entering 
the atmosphere. The resulting scatterers might well al] have the same total 
number of electrons, and their subsequent diffusion would lead to a regime 
of scatterer sizes which was unchanging with time. On the other hand, it 
might be more appropriate to think of these scatterers as maintained against 
debilitation by recombination and diffusion by a continual regeneration of 
ionization, so that the resulting scatterers would have about the same peak 
electron density, and have a range of sizes introduced by the varying equili- 
brium conditions between the production and decay of ionization. Again, the 
auroral ionization might be produced by the incoming particles as a large, 
rather uniform blanket, the irregularities in the ionization being subsequently 
introduced by the action of turbulence or wave motion of the atmosphere. 
This would likely produce scatterers with about the same peak electron 
density, and with a distribution of scatterer sizes which is characteristic of 
turbulence. Almost certainly the actual situation is a combination of these 
various production mechanisms, the resulting assembly having a distribution 
of both size and peak electron density of the scatterers. 

Consider first an assembly of scatterers produced by the continuous for- 
mation of very small clumps of ionization which subsequently diffuse, the 
number of electrons in each clump being the same, and their occurrence being 
random both in space and time. It is assumed that the process has been going 
at a constant rate for a time sufficiently long to make the distribution of 
scatterers with size constant with time. Then the number of scatterers per 
unit volume per unit range of scatterer sizes is 


(34) n(a) = Ma/2D, 


where M is the number of scatterers formed per m* per second, and D is the 
diffusion coefficient in m?/sec. If a scatterer of size a has scattering coefficient 
a(a), the scattering coefficient per unit volume for the assembly is 


(35) o,= J m(aro(arda = “ J aotaaa, 


where dp is a limit put on the initial size of the scatterers by the physical 
processes involved. No irregularity in the auroral ionization can be smaller 
than the electron mean free path, or the electron gyro-radius (which ever is 
the smaller) for directions transverse to the field lines. It appears likely that 
this will limit a) to values greater than about 0.3 meter in directions trans- 
verse to the magnetic field, and to values greater than about 0.3 to 1.0 meter 
for directions parallel to the magnetic field at ionospheric heights of about 


100 km. 
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To evaluate (35) consider the case where ¢=0, and ignore the small con- 
ductor theory, since it is valid only for very small, strong scatterers, which 
are not important physically. If ao is sufficiently large, the scatterer will never 
be thicker than the skin depth at any frequency, and o(a) will be given entirely 
by (7). From (5) and (20), this will occur if 


: 4req _ 
(36) ayo > ay fe 7. 





For smaller values of a, o(@) will still be given by (7), if the frequency is 
sufficiently high so that the peak electron density is less than the critical 
density. From (5), (9), and (17), this occurs if 


1/2 


(37) ao<vy, k> (y/ai) 
At lower frequencies both mechanisms are valid at different sizes, the size 
at which the transition occurs being given by 
(38) a = (y/k*)'”. 
Using these results, and numerically integrating the ray theory expressions 


it is possible to evaluate (35) for different values of ao as a function of fre- 
quency. Figure 6 shows the resulting normalized scattering coefficient per 


Normalized Scattering Coefficient, = 
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Wave number, k 


Fic. 6. Normalized scattering coefficient, 2=Do,/h?M, as a function of frequency, for 
an assembly of diffusing scatterers (6 =0, y =0.64a,°). 
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unit volume for various values of @ and y (proportional to the total number 
of electrons in each scatterer). The values used were chosen so that for each 
curve the ray theory made a contribution only for wave numbers less than 
0.8 m—. This leads to a transition at this wave number from a region where 
ray theory is important to one where scatter theory is important. The curve 
for d)=1 m illustrates a case where the scatter theory is valid throughout 
the frequency range. As the minimum size increases the sharpness of the 
transition also increases, but a comparison with Fig. 4 indicates that the 
transition is much less abrupt than for an assembly consisting of similar 
scatterers, all of the minimum size. 

The other class of assemblies considered in this paper is one consisting of 
scatterers of different sizes, all having the same peak electron density. This 
is equivalent to integrating (with a weighting factor) along a vertical line on 
Fig. 3. If k” < 1, there is a transition from scatter theory to ray theory at a 
value of a’ of approximately 2. If k’ > 1, scatter theory is valid throughout. 
The scattering coefficient can be written 


(39) oy J n(a)o,(a)da, =i, 


ac 


a? 


| eae’ J n(a)on(a)da, k<, 


Consider a size distribution of the form 


(40) n(a) = nae? (Ain, 

This distribution has a maximum value at dy, falling to zero both at a=0 and 
a=, Several examples of this family of distributions are shown in Fig. 7. 
As p approaches infinity the distribution becomes a very narrow peak about 
@=Q,, and thus approaches the case of an assembly of similar scatterers. For 
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comparison, it is of interest to rewrite these size distributions so that they all 

have the same total number of electrons per unit volume, Q: from (5) and (40), 
(p+4) /2 a 

(41) n(a) = write _pitcas:,,gugeni : 


rigornay( P44) Nohany* 
2 





where I(x) is the gamma function. Here the minimum size, a, has been 
assumed to be zero. The scattering coefficients from these assemblies can be 
obtained by substituting (23), (24), and (41) in (39). In general the cases for 
k’ < 1 can only be solved by numerical integration. However, in practice, 
if m is sufficiently large, and a,k, is greater than 2, the scatter theory term 
can be neglected. Under these circumstances, little error will be introduced 
by assuming the minimum size to be zero. The resulting expressions can be 
integrated to yield the following scattering coefficients: 


(42) ex —Oihseds —-, k >, 


4a’2h” , ‘ (p+7) /2 
we 14 S88 ( -os* ¢-+h* sin” 6) | 





Qhk Byf”(k’*)In(1/k" 





(48) ~ Nodalcos’(¢+6) +h sin aE] re 
where 
(e) 
(44) hig ae oc eeaett, 
air girprap(2+4) 
prar(e+3) 
(45) ge eee 
ergy (044) 


and where a’=k,a and k’=k/k,, as before. As p approaches infinity the 
assembly described by (41) approaches an assembly of similar scatterers, all 
of size d,. Using the expansion of the gamma function for integral values of 
p, Wallis’ product for 7, and the limit expression for e”, it can be shown that 
A, approaches 1/167? and B, approaches 1/42*” and the scattering coefficients 
given by (42) and (43) become the expressions for an assembly of similar 
scatterers of size d,, having a total of Q electrons per unit volume, just as one 
would expect. Figures 8 and 9 show the normalized scattering coefficient per 
unit volume as a function of wave number for different values of p, for two 
values of d,, and for ¢=0. The distribution of sizes has little effect on the 
nature of the critically reflected coefficients, but changes markedly the fre- 
quency dependence of the weakly scattered waves. The lower values of p 
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Normalized Scattering Coefficient, =’ 





o-2 0-6 1-0 14 


Normalized Wave number, k 


Normalized Scattering Coefficient, =’ 





O-2 0-6 1-0 1-4 
Normalized Wave number, k 


Fic. 8. Normalized scattering coefficient, =’=Noo»/Qhke, as a function of frequency, for 
an assembly of gaussoids with a distribution of sizes given by (41) (am’=4, ¢=0). 

Fic. 9. Normalized scattering coefficient, =’ =Noo,/Qhke, as a function of frequency, for 
an assembly of gaussoids with a distribution of sizes given by (41) (am’=3, ¢=0). 
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decrease the frequency dependence over that for the assembly of similar 
scatterers. For k’=1 this effect is reversed if aj, is less than 1.3. Decreasing 
the size of » then makes the scattering coefficient fall more rapidly with fre- 
quency. Such a low value of aj, does not allow the complete integrals from 
which (43) was derived to be used, and therefore cases such as this are not 
shown in these diagrams. 

CON CLUSIONS 

Scattering coefficients have been determined for several families of models 
of auroral ionization. These models appear to encompass many of the possible 
types of auroral models. The assembly of similar gaussoids is basic to any 
consideration of auroral ionization, covering the whole range of sizes, shapes, 
and frequencies. The effects of distributing the sizes and strengths of these 
scatterers are illustrated by the other assemblies considered. It should be 
noted that none of these assemblies involve a distribution of sizes with maxima 
at more than one size. (A simple example of a distribution with maxima at 
two sizes would be one consisting of -scatterers of two different sizes.) The 
characteristics of such distributions can be quite different from those con- 
sidered in this paper. 

In this paper equations and curves have been presented, representing the 
properties of the scattering coefficients of these different auroral models. An 
experimental investigation does not usually provide a direct measurement 
of the scattering coefficient of auroral ionization, so the results of this paper 
must be interpreted in terms of the conditions of each particular experiment. 
In Part II (Moorcroft 1961) the available experimental data are examined, 
particularly with regard to their suitability for interpretation in terms of models 
of auroral ionization, and some of the results are used to draw conclusions 
concerning the nature of auroral ionization. 
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MODELS OF AURORAL IONIZATION 
PART II. APPLICATIONS TO RADIO OBSERVATIONS OF AURORA! 


D. R. Moorcrort 


ABSTRACT 


In this paper the available experimental evidence concerning radio reflections 
from aurora is examined in relation to the reflection characteristics of the models 
of auroral ionization discussed in Part I. 

The existence of critical reflection from auroral ionization at frequencies 
between 30 and 50 Mc/s appears to be established. This implies electron densities 
as great as 3X10" electrons/m*. It is shown that the observations are consistent 
with a model consisting of irregularities of ionization elongated parallel to the 
earth’s magnetic field in a ratio of between 5 and 10 times, and having sizes 
transverse to the field lines of the order of a few meters. Some of the observations 
require the irregularities to be distributed in size. It is clear that there is a need 
for more relevant measurements on radio reflections from aurora to specify the 
characteristics of the ionization more precisely. 


INTRODUCTION 


In Part I (Moorcroft 1961) several general models of auroral ionization 
were presented, and their reflecting characteristics for radio waves were dis- 
cussed. These models all include consideration of both critical reflection and 
weak scattering of radio waves, and, taken together, provide a general treat- 
ment, covering a range of differing pictures of auroral ionization as seen by 
means of radar techniques. A re-examination of the relevant experimental 
evidence in the light of this general treatment should provide information 
regarding the nature of the auroral ionization, and should indicate other 
types of experiments that would distinguish the various reflection mechanisms. 
It will be shown that such an examination settles the question of the existence 
of critical reflections in auroral reflections. 

A description of auroral ionization can be given conveniently in terms of 
the size and distribution of sizes of the individual irregularities, their shape, 
their strength (i.e. peak electron density), the sharpness of their boundaries, 
and the number of scatterers per unit volume. In the literature there is con- 
siderable disagreement upon many of these aspects of the model. The estimates 
of size range from a fraction of a meter (Booker 1956) to hundreds of meters 
(Kaiser 1956). Some types of auroral ionization have been suggested to consist 
of isotropic scatterers (Collins and Forsyth 1959), while estimates of the 
elongation of anisotropic scatterers range from an axial ratio of 4:1 (Forsyth 
and Vogan 1957) to about 40:1 (Booker 1956). The most fundamental and 
important difference relates to the strengths of the scatterers. From the 
beginning of these studies, when Herlofson (1947) suggested that it would be 
physically unreasonable to interpret the auroral reflections obtained by Lovell, 
Clegg, and Ellyett (1947) as due to critical reflections, with the higher electron 
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densities that that interpretation implied, the role of critical reflection in 
producing auroral echoes has been widely discussed. More recently, Booker 
(1956) has suggested that weak scattering is the reflection mechanism at 
frequencies as low as 30 Mc/s, while Kaiser (1956) has visualzed critical 
reflection being important up to frequencies as high as 100 Mc/s. Forsyth 
(1953), Forsyth and Vogan (1957), and Collins and Forsyth (1959), on the 
basis of multiple-frequency measurements, have concluded that for frequencies 
between 30 and 100 Mc/s, there are times when the frequency dependence is 
such that a transition from strong to weak scattering must take place within 
this frequency range. They also suggest that at other times the whole fre- 
quency range may be reflected by either weak scattering or critical reflections. 


EXPERIMENTAL EVIDENCE 

Not all of the experimental evidence on radio-auroral reflections is directly 
relevant to a discussion of the properties of auroral ionization. Many studies 
have been undertaken to reveal the diurnal, annual, and solar variations of 
radio-aurora, and its correlation with various other geophysical phenomena: 
these are not concerned with the process of reflection. Other observations 
which do relate to the reflection mechanism (and thus to the nature of the 
ionization responsible for the reflection) are complicated by the superposition 
of effects due to the properties of the radio system being used, and other 
effects resulting from the spatial distribution of auroral ionization. 

In particular, there are several limitations imposed on the interpretation 
obtained from a radar system which, while perhaps obvious, are nevertheless 
worthy of mention. Since the pulse length of a pulsed radar is finite, all echo 
information is integrated over distances equal to the spatial length of the 
pulse. Similarly, the finite repetition rate means that changes in the scattering 
region which take place in times shorter than the repetition period cannot be 
detected. The antenna pattern of the radar is affected to some extent by the 
terrain over which it is directed. For this reason it is sometimes difficult to 
specify theoretically the precise antenna pattern. Even experimental determi- 
nations of the antenna pattern are not usually sufficiently accurate to allow 
accurate information to be obtained about the height of the reflecting region 
(this is particularly true at low frequencies). 

The relevant experimental evidence available is now examined in the light 
of the considerations mentioned above, under the following five headings: 
(i) classification of auroral echoes; (ii) height of reflecting regions; (iii) echo 
amplitude; (iv) aspect sensitivity; (v) frequency dependence. 

Almost all workers, in discussing the nature of echoes received, have classi- 
fied the reflections into two general categories, usually termed ‘‘diffuse’’ and 
“discrete”. The diffuse echoes are long lasting and widely spread in range 
(several hundred kilometers). Discrete echoes exhibit structure which is 
characteristically tens of kilometers in size, sometimes limited by the finite 
pulse length, which is typically in the range 12 km (for a pulse length of 
40 usec) to 600 km (for 2 msec). Some workers (Booker, Gartlein, and Nichols 
1955) do not find any sharp demarcation between short and long spreads in 
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range. Aspinall and Hawkins (1950), Hellgren and Meos (1952), Currie, 
Forsyth, and Vawter (1953), and Meos and Olving (1958) suggest two quite 
distinct types of echoes. Bullough and Kaiser (1954, 1955) obtained both 
types of echoes but suggest (Bullough and Kaiser 1954) that at least on one 
occasion the diffuse echo was the coalescence of a large number of discrete 
echoes. Presnell et al. (1959), operating a narrow beam antenna on 218, 398, 
and 780 Mc/s, find two quite different types of auroral echoes, and associate 
them with quite different distributions in space. The ‘‘diffuse’’ scattering 
centers seem to be distributed in an almost horizontal layer of wide extent, 
while the “‘discrete’’ scattering centers seem to lie along lines of the magnetic 
field. They found that the diffuse type of echoes occurred during the daytime. 
There may be no relation between the “‘diffuse’’ and ‘‘discrete’’ classifications 
of Presnell and those of the earlier workers. It is tempting to associate the 
diffuse radio-aurora of Presnell et al. with the type A; radio-aurora observed 
by Collins and Forsyth (1959) on their bistatic system. Of the three types of 
auroral echoes observed they found only A; occurring in the daytime. How- 
ever, it appeared to be located in a layer rather low in the ionosphere (that is, 
probably below 100 km), and was apparently not aspect sensitive, which 
contrasts with the observations of Presnell e¢ al., who found the diffuse radio- 
aurora to be quite aspect sensitive, and located at about 100-km height. It is 
not clear how all these different observations are related, but there seems to 
be evidence of at least two distinct types of auroral ionization, which may 
have quite different characteristics. 

The ionospheric parameters used in conjunction with a model of auroral 
ionization depend very markedly upon the height at which the reflecting 
ionization is located. These heights are very difficult to determine accurately 
with a v.h.f. radar alone, since usually the width of the antenna pattern and 
the uncertainties in it are so large as to make height estimates ambiguous. 
Studies by various workers suggest that most radio-auroral reflections come 
from ionization located between 100 and 120 km. The most accurate determi- 
nation of height available is that of Unwin and Gadsden (1957), who used an 
antenna with a large number of well-defined lobes (produced by interference 
between the direct signal and the signal reflected from the sea, over which 
the antenna looked from a high cliff). They obtained very consistent results, 
indicating reflection from a layer many hundreds of kilometers in extent and 
often varying in height by only 4 or 5 km along its extent. They concluded 
that the reflecting layer was less than 2.5 km in thickness, and had a height 
of about 110 km, with individual heights ranging from 100 to 120 km. How- 
ever, recent observations of the Stanford group (Schlobohm ef al. 1959) at 
106 Mc/s suggest that at low latitudes radio-auroral echoes can be obtained 
from much greater heights, perhaps as high as 300 km. 

Although, in principle, amplitude measurements of radio-auroral echoes 
contain much useful information about the mechanism of reflection, the 
information is very difficult to separate from the many other factors that 
contribute to the echo amplitude. Finally, one would like to deduce the scatter- 
ing coefficient per unit volume, whereas a radar gives the total scattering 
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coefficient. To obtain the latter from the former one needs to know the location 
and extent of the scattering volume, and its depth. Because of these com- 
plexities, it is only possible to indicate by amplitude measurements that radio- 
aurora is quite transparent to v.h.f. radio waves. If the total scattering 
coefficients observed by different workers are expressed in square kilometers, 
the extent to which radio-aurora is transparent will become evident. The 
total scattering coefficient is easily expressed in terms of the radar parameters 


(Bowen 1954): 


_ (an) P 
( 1 ) = Gy P. ’ 


where G is the antenna gain (with respect to an isotropic radiator), P; is the 
transmitted power, P, is the received power, r is the range, and A is the free 
space wavelength. (Note that ¢ is defined by equation (4) of Part I. The radar 
cross section often used in the literature is just 4m times this value.) Using 
(1) it can be seen that Aspinall and Hawkins (1950) detected a total scattering 
coefficient of 0.36 km? at 72 Mc/s. Forsyth (1953) found values of ¢ as large 
as 4.5 km? on both 56 and 106 Mc/s. Little, Rayton, and Roof (1956) quote 
Dyce as saying that at 50 Mc/s the equivalent target area is about 30 km? for 
strong aurora, which corresponds to ¢=2.4 km’. It is obvious that the value 
obtained depends to a considerable extent on both the width of the beam, and 
the integration introduced by the pulse length. Another way of illustrating 
the over-all transparency of radio-aurora is to define an effective reflection 
coefficient, defined as the ratio of the echo strength received to that which 
would be received from a specularly reflecting metal sheet located at the 
radio-auroral range. Little, Rayton, and Roof (1956) give a value of 10-4 for 
the effective reflection coefficient of strong aurora at 100 Mc/s. These measure- 
ments indicate that, on the average, radio-aurora is transparent; they do not 
imply necessarily that small regions of it can not be highly reflecting. 

One of the more controversial aspects of the radio-auroral literature is the 
interpretation of the observed angular dependence of auroral echoes. Part 
of the difficulty is that several different factors contribute to the observed 
distribution and their relative importance is not easily determined. The term 
aspect sensitivity is generally used to imply an angular dependence of echoes 
produced as a result of the nature of the reflection process itself. 

Before the observed angular dependence can be interpreted in terms of 
aspect sensitivity, the other factors contributing to the angular dependence 
must be removed. It is reasonable to expect that radio-aurora is distributed 
geographically much like optical aurora. If this is the case, then simply from 
the geographic situation one would expect that for low latitude stations the 
echoes would come from the north at low elevation angles. Bullough and 
Kaiser (1954) and Kaiser (1957) have shown that their echoes lie distributed 
along parallels of geomagnetic latitude, suggesting that the dominant effect 
in determining the echo location is the geographic distribution of the ionization 
(presumably closely associated with the optical aurora), and not an aspect 
sensitive reflection mechanism. Echoes may also be restricted to the northern 
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azimuths as a result of the actual volume of ionization presented to the antenna 
beam as it looks in different directions (Forsyth 1960). Suppose that the 
ionization is distributed along lines of geomagnetic latitude. Then Fig. 1 
shows how the difference in volume scattering operates to restrict echoes to 


North 





¥ 
pulse width 


Fic. 1. The effect of the volume distribution of auroral ionization on the angular dependence 
of radio reflections (after Forsyth 1960). 


the north, even for stations which are quite close to the auroral zone. Forsyth 
has shown that this can account for all of the azimuthal aspect sensitivity 
observed. Although it is apparent that the ionization is not likely to be dis- 
tributed strictly in arcs along the geomagnetic latitudes, this effect is never- 
theless likely to be of some importance in determining the final angular 
dependence of the auroral echoes. This same effect (to a reduced extent) is 
likely important in explaining the variation with elevation angle. 

As an alternative to aspect sensitivity for the explanation of the limitation 
of echoes to low elevation angles, it has been suggested (Currie, Forsyth, and 
Vawter 1953) that the ionization responsible for radio reflections is underlain 
with absorbing ionization. It is sufficiently patchy so that from a distance the 
reflecting ionization is visible between the patches of the absorbing layer. At 
greater elevations the reflecting ionization is obscured by absorption. There 
is some question as to how effective this would be. It would require a very 
particular distribution of ionization below the radio-aurora to extinguish 
echoes at elevations as low as 15 degrees, and still leave them visible at 12 
degrees. 

In the light of the above discussion, most of the low latitude observations, 
which localize radio-auroral occurrence to the north and to low elevation 
angles, can be explained without any recourse to aspect sensitivity. One excep- 
tion to this statement is the low latitude, high elevation observations of 
Schlobohm et al. (1959), who detected echoes on 106 Mc/s at heights of 200 
to 300 km from a latitude of 43°. They suggest that the only reason these 
echoes could be seen was that the observations were made normal to the lines 
of the magnetic field, along which they assume the ionization to lie. Dyce 
(1955a) has operated an auroral radar on 51.7 Mc/s at Point Barrow, which is 
about 200 km north of the visual auroral zone maximum. He found that more 
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than 90% of the echoes were obtained from directions to the north of east 
and west. He presents this as strong evidence for the importance of aspect 
sensitivity in determining the azimuthal distribution. Recently McNamara 
(1959) has studied the azimuthal distribution of auroral echoes on 48 Mc/s at 
Resolute Bay, Baker Lake, Saskatoon, and Ottawa, using four 90° corner 
reflectors at each station. The results he obtained are shown in Fig. 2. The 
relative total occurrence at the different stations were in the ratios of 1:8:24:8 
for Resolute Bay, Baker Lake, Saskatoon, and Ottawa, respectively. For 


2 


Resolute Bay Baker Lake- 


Saskatoon Ottawa 


Fic. 2, Azimuthal dependence of auroral echoes (after McNamara 1959). 


Resolute Bay and Baker Lake, both of which are far inside the auroral zone 
the fewest echoes were found to come from the north, although the total 
number of echoes received was considerably less than at Saskatoon. This 
suggests that although aspect sensitivity plays some part in determining the 
azimuthal distribution of echoes, it is not enough to counterbalance the effect 
of the geographic distribution of ionization at high latitudes. More significant 
estimates of aspect sensitivity can be made using the variation of echo occur- 
rence with elevation angle, since the effects of geography are not as important 
as for the azimuthal distribution. Apparently echoes have not been seen at 
elevation angles above 30°, regardless of the location of the radar. The effect 
of the volume distribution of scatterers is not likely to be strong enough to 
cut off echoes at such low elevations, although the strength of the effect will 
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depend very much on the thickness of the scattering region. Although absorb- 
ing regions may be capable of explaining observations of this type at fre- 
quencies below, say 100 Mc/s, they will have little effect on high frequencies. 
It appears likely that the most important factor in limiting auroral echoes to 
low elevation angles, particularly at high latitudes, is aspect sensitivity. 
However, if aspect sensitivity is responsible for the elevation angle cutoff, 
then this maximum angle should be a sharp function of the azimuth, an 
effect which does not appear to have been observed. Presnell et al. (1959) 
have reported observing echoes at College, Alaska, on 218, 398, and 780 Mc/s 
at angles up to 12, 10, and 6 degrees, respectively, off perpendicularity with 
the magnetic field lines. These measurements were not simultaneous, were 
not made with antennas of similar beam widths, nor were the different fre- 
quencies all operated for the same length of time, so it is difficult to draw any 
reliable conclusions about the frequency dependence of aspect sensitivity 
from these results. Recently Greenhow, Neufeld, and Watkins (1960) have 
reported measurements of the angular dependence of auroral echoes made 
with the Jodrell Bank radio telescope on 36 Mc/s. In analyzing the results 
they made a correction for the probable geographical distribution of the 
ionization by using a distribution for optical aurora given by Vestine. They 
deduced the angular dependence of the scattering from the number of echoes 
received at different aspect angles, correcting for the relative times spent by 
the beam in different parts of the sky. The observations and corrected results 
which they obtained are shown in Fig. 3. 


Number of echoes 





g in degrees 


Fic. 3. Aspect sensitivity of auroral echoes observed at Jodrell Bank on 36 Mc/s (after 
Greenhow et al, 1960). The histogram indicates the experimental results; the dots are corrected 
values. 
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Further qualitative, if not quantitative evidence for the existence of aspect 
sensitivity is provided by measurements using bistatic radar systems. Because 
of the different geometry of the bistatic system it is possible to reduce the 
effect of the geographical distributions of ionization, and the effect of the 
volume distribution of the scatterers. That is, a bistatic system can ‘‘see”’ 
scattering centers over a sufficiently wide area that the conditions related to 
the location of the ionization (near to the auroral zone, in long arcs aligned 
along parallels of geomagnetic latitude, perhaps) can be satisfied and still 
leave a wide range of different conditions of aspect sensitivity to be fulfilled. 
In the absence of aspect sensitivity, the most favorable condition for forward- 
scattering is straight along the path, since isotropic scatterers will scatter the 
greatest amount of energy in the forward direction. If the scatterers are aspect 
sensitive, echoes could be detected from directions off the path of the system 
if the aspect conditions were more favorable there. Collins and Forsyth (1959) 
have reported a detailed investigation of radio-aurora using several bistatic 
systems. The three types of radio-aurora they distinguished were found to 
have quite different characteristics, the type A; (mentioned before in the 
discussion of the classifications of auroral echoes) being apparently isotropic, 
while types A2 and A, showed evidence of progressively more aspect sensitivity. 

It is generally recognized that the only satisfactory way to study the fre- 
quency dependence of radio-auroral echoes is to use scaled systems. That is, 
to use antenna systems on the different frequencies which have as nearly as 
possible identical beams, to use the same pulse lengths and repetition rates, 
and to have the same sensitivity on all frequencies. Further, it is necessary for 
the observations to be carried out simultaneously if the information is to be 
useful, since observations at different times may be of quite different strengths 
of radio-aurora, and involve different mechanisms of reflection. 

Simultaneously, multiple-frequency observations which have been made to 
study the frequency dependence of echo strengths have been made by Forsyth 
(1953, see also Currie, Forsyth, and Vawter 1953), Harang and Landmark 
(1954), Dyce (1955)), Forsyth and Vogan (1957), and Collins and Forsyth 
(1959), and Flood (1960). The experiments of Harang and Landmark, and 
Dyce were both on two frequencies, with quite different antenna patterns for 
the two frequencies. The work of Forsyth (1953) was on 56 and 106 Mc/s. 
The pulse widths were different, and the 56 Mc/s beam was twice as wide as 
that on the higher frequency. In analyzing the data Forsyth attempted to 
correct for these differences, but this involves certain assumptions, since the 
correction depends upon the distribution of the ionization and the mechanisms 
of reflection. The effect of the different pulse widths on the two frequencies 
was corrected for by defining a total scattering coefficient per kilometer range 
in terms of amplitudes normalized to have a value of unity at the mean noise 
amplitude. Out of a total of 1204 periods of observation during which echoes 
occurred, 922 corresponded to echoes on 56 Mc/s only, 23 to echoes on 106 
Mc/s only, the remaining 259 being simultaneous occurrences on the two 
frequencies. For the simultaneous observations the amplitude was most 
frequently greatest for 56 Mc/s, although the ratio of amplitude took on all 
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measurable values. It appears likely that those observations which indicated 
a smaller signal on the lower frequency were affected by the presence of 
selective absorption. Both Forsyth and Vogan, and Collins and Forsyth 
operated exactly scaled systems on 32.22, 39.22, and 48.82 Mc/s, using forward- 
scattering. Forsyth and Vogan illustrated their three-frequency observations 
by means of what they called an XY diagram. If the amplitudes and wave- 
lengths of 32.22, 39.22, and 48.82 Mc/s are designated by Aj, Ao, A3, and 
Ai, Az, Az, the two ratios X and Y are given by 


— (A2/d2) , — (Ai/ri) 

(A3/Xs3)’ (Ao/X2) © 
The division of each amplitude by the wavelength gets rid of the wavelength 
dependence in the radar equation (1), so that (2) becomes 


rh rely 
(3) X = oe Y= oak? 


where oj, a2, and o3 are the scattering coefficients for 32.22, 39.22, and 48.82 
Mc/s, respectively. Forsyth and Vogan plotted X vs. Y (XY diagram) for 
observations obtained during two nights of radio-auroral activity, and these 
results are shown in Fig. 4. Small values of both X and Y indicate a gradual 
frequency dependence, whereas large values of both X and Y indicate that the 








(2) X 





Fic. 4. XY diagram showing the results of Forsyth and Vogan (1957). 
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scattering coefficient is getting rapidly smaller with frequency. A large value 
of X together with a small value of Y indicates that the difference between 
Az and A; is much greater than that between A; and Az», while a large value 
of Y with a small value of X indicates the reverse. The observations of Collins 
and Forsyth (1959) are essentially the same as those of Forsyth and Vogan, 
except that they find that type A, echoes seem to give a large value of Y and 
small X, while type A; has both X and Y small. 

Flood (1960) has reported results obtained from a simultaneous, three- 
frequency study of radio-aurora on 49.17, 143.5, and 226 Mc/s. The two lower 
frequencies were operated on essentially scaled systems, while the highest 
frequency used an antenna with a considerably narrower beam, but which was, 
as nearly as possible, arranged to illuminate the same portion of the ionosphere 
at the 100-km level. The ratio of total scattering coefficients (as defined by 
(1)) on 49.7 and 143.5 Mc/s varied from values less than 1 toas large as 1700 (or 
more), with the over-all average ratio being 41. Echoes on 226 Mc/s occurred 
only with strong echoes on the other frequencies. On these occasions, the ratio 
of the scattering coefficients on 143.5 Mc/s to that on 226 Mc/s varied from 
about 2 to 230, with the average ratio being about 19. It is interesting to note 
(Flood 1959) that the strongest echoes on 49.7 Mc/s were not accompanied 
by echoes on 226 Mc/s. 


RADIO-AURORAL REFLECTION MECHANISMS 


These experimental results can now be examined in relation to the models 
of auroral ionization developed in Part I. Let us first see how well the results 
can be explained by a model of auroral ionization consisting of an assembly 
of similar gaussoids. Of particular interest is the comparison of the frequency 
dependence of the waves scattered from such a model with that found experi- 
mentally, since it will likely provide evidence relating directly to the relative 
importance of weak scattering (Booker’s theory) and critical reflection. 


Frequency Dependence 

In view of the variety of theoretical curves which can be constructed, it is 
necessary to have simultaneous measurements with scaled systems on as 
many frequencies as possible. The most satisfactory results for this purpose 
are those of Forsyth and Vogan (1957), and Collins and Forsyth (1959). It 
should be borne in mind, however, that these experimental results were 
obtained using forwardscatter systems. Whether or not the following analysis 
applies in detail to forwardscattering, it nevertheless indicates a way of inter- 
preting any future three-frequency backscatter experiments. The formulae 
developed in Part I can be used with (3) to give expressions for X and Y. 
These will be different, according as the critical wave number, k,, for the 
scatterer is less than k;, between k; and k», between ky and k;, or greater than 


kz. From Part I if ke < ky, 


(4) Xi ca e-364a? 


(5) Y= 0-220? 
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If ky < ky <& ko, 





(6) Xp= coset, 
\ 
Ado) a 
(7) oa ar if J 0.6760? 


(If Ye < Yi, use Y;. This is necessary, because at frequencies just below the 
critical frequency, the critically reflecting surface is too small to prevent the 
wave from penetrating, and the scatter theory then becomes valid. It corre- 
sponds to the transition region indicated in Fig. 3 of Part I. The other 
provisional cases which are given below are also a result of this effect.) 
Hikes << Re < hy, 


RB? 2\ 13 
af) [in e 
(8) X3 = Be oe 
(#) (#) 
s( 2 72 
(9) es (i) n(f) . 
(#2)\in(#3) 
(#3) ™ (#) 
(If X; < X, use X1; if Ys > Yo use Yo; if Y2 < Vi, use Yi.) 
If ks < ke, 
2 2 } 
yet 
(10) So wid | _ Ve! | 
CoG) 
ti) (#)\ 
@ In\ 32) 
(11) (8) in(#) * 
“\es NV Ni 


(lf X, > Xz3use X3; if X; < X, use X;; if Ye > Yo use Y2: if Y; < VY; use Y;.) 
The resultant values of X and Y (which are functions of a, the scatterer 
size) are plotted in Fig. 5 for four different strengths of scatterers, illustrating 
each of the four cases given above. The curve for k,=0.818 m~ (critical 
frequency 39.0 Mc/s) represents the upper limit for values of Y associated 
with an assembly of similar gaussoids. Comparison of these curves with the 
experimental results obtained by Forsyth and Vogan (1957) as shown in Fig. 
4 reveals many experimental points which are not encompassed by these 
theoretical curves. No points which lie above the broken line shown in Fig. 4 
can be explained by an assembly of similar scatterers. In particular, values of 
Y greater than X still require explanation. 
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Fic. 5. XY diagram for the frequencies used by Forsyth and Vogan (1957), showing 
curves for various assemblies of similar gaussoids. Values of a range from 1.0m at X=1.4 
to 2.5m for X =10, 


There are several ways in which this discrepancy might be explained: 
(i) the results are for forwardscattering, while the theory is for backscattering; 
(ii) the scatterers are not accurately represented by gaussoids; (iii) some other 
effect such as multiple scattering or absorption plays a role in the scattering 
process; (iv) the scatterers are not similar, but their distribution rather 
involves various sizes and strengths. Each of these possibilities will be con- 
sidered in turn. 

Although it is difficult to develop any accurate theory for forwardscattering, 
it is not difficult to make a rough estimate of the resulting modifications. To 
a first approximation, the scattering coefficient for a critically reflecting 
scatterer reflecting at an angle of 90°, say (a likely angle for the experiments 
of Forsyth and Vogan), will not be substantially different from that for back- 
scattering. There may be a slight increase in the value of o due to the greater 
size of the critically reflecting surface at oblique incidence. This also results 
in the critical reflection being the mode of reflection at higher frequencies 
than for backscatter. For weak scattering, the formulae of Part I include 
the case of forwardscattering. For scattering in the plane normal to the 
scatterer axis, (6) of Part I becomes 
(12) aie men er. 
where @ is the angle between the incident and received directions. If @= 180°, 
we have backscattering, and (12) becomes 

c= ia<n'. 


as before. If 6=90°, then 


(13) 
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Thus, approximately, all the expressions (4) to (11) can be corrected for 
forwardscattering by dividing the exponent by 2. The resulting expressions 
all have exactly the same value, for a given a, as the backscatter expressions 
with a value of a smaller by a factor 2—”?, with the exception of (7) and (8), 
both of which have values which are one-half of the corresponding values 
for 2-”2a. This has the effect of sliding the rising curves further to the right 
of the XY diagram. It is thus clear that this correction provides no basis for 
an explanation of the experimental results above the Y=X line. 

It is possible that in using the gaussoid physically unreasonable characteris- 
tics have been introduced into the model, and for this reason the model is 
not applicable. If this is the case, the difficulty could involve only the diffuse- 
ness of the boundary of the scatterer, since the shape will have little effect on 
the frequency dependence of the scatterer, and the size and strength can 
already be freely varied. It appears likely that for scatterers, which are only 
a few mean free paths in size, the nature of the electron density distribution 
will almost certainly be determined by diffusive process, leading to a gaussoidal 
distribution. Larger irregularities may be formed by other processes, and 
result in boundaries which are less diffuse, relative to the size of the irregu- 
larity, than a gaussoid. However, the sizes shown in Fig. 5 indicate that only 
the small size scatterers are of importance for values of X and Y less than 5, 
and hence it is unlikely that there is any basic error in assuming the scatterers 
to be gaussoids. 

In the treatment of scattering given in Part I no consideration was given 
to the possible influence of multiple scattering. Recent calculations (Moorcroft 
1960) suggest that this may be an important effect, particularly in reflection 
from weak scatterers. The effect of this appears to be to decrease the frequency 
dependence of weak scattering, and this might justify an upward shift of some 
of the curves into the region above Y=X on the XY diagram. 

Another effect which has not been taken into account is absorption of 
the radio waves. There is considerable evidence for this in the observations 
of frequency dependence (Forsyth and Vogan 1957; Flood 1960). The effect 
of this on a point on the XY diagram is to move it downwards to the left 
along a line that makes an angle of 45° or more with the X-axis. Such an 
effect cannot move the curves of Fig. 5 above the Y=X line. 

It seems likely that the correct way out of this difficulty lies in choosing an 
assembly of scatterers having a distribution of sizes and strengths. There are 
several different types of distribution which we could consider, involving 
strength alone, both size and strength together (as in the assembly of diffusing 
scatterers, which was discussed in Part I), or sizes alone (as in the various 
other assemblies considered at the end of Part I). To explain the observations, 
a satisfactory distribution must give values of Y greater than X, for small 
values of both X and Y. This implies that the lowest frequency is reflected 
much more strongly than the other two frequencies. Distributions such as 
those considered in Part I have this type of frequency dependence only when 
radio waves of the lowest frequency are critically reflected. However, none 
of the assemblies considered in Part I had size distributions with maxima at 
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more than one size. There are assemblies of this type which have the observed 
frequency dependence without requiring critical reflection on any of the 
frequencies. Suitable distributions of this kind appear unlikely to occur in the 
ionosphere, and will not be considered in this paper, but they cannot be 
excluded as possibilities by the observations considered here. 

Consider an assembly with scatterers of the same size but having different 
peak electron densities. With a little thought it can be seen that this will 
smooth out the transition between the two mechanisms, since now the number 
of scatterers critically reflecting at any given frequency will be a decreasing 
function of frequency, instead of dropping abruptly to zero at one particular 
frequency. Such a distribution will be no improvement over the assembly of 
similar scatterers. 

The assembly of diffusing scatterers is a physically reasonable example of 
an assembly involving a distribution both of size and strength. Values of X 
and Y for the three frequencies used by Forsyth and Vogan for assemblies 
of this type are given by curve (6) shown in Fig. 6. The parameter which 





Fic. 6. XY diagram for the frequencies used by Forsyth and Vogan (1957), showing curves 
for various assemblies of gaussoids. (a) Scatter theory for an assembly of similar gaussoids. 
(6) Assembly of diffusing scatterers. (c) Assembly of gaussoids having a size distribution given 
by (41) of Part 1; an=4.3 m, ke=0.7 m=. (d) Assembly of gaussoids having a size distribution 
given by (41) of Part I; an=5.8 m, kp =0.7 m7. 
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varies along the curve is do, the minimum scatterer size. For comparison the 
scatter theory curve is also shown (curve (a)). This distribution is of little 
help in extending the theoretically predicted regions of the XY diagram. 

The distributions considered at the end of Part I involve a distribution of 
sizes only, the scatterers all having the same peak electron density. This 
promises to be useful, since, roughly speaking, the weak scattering will tend 
to be characteristic of the small scatterers, while the nature of the critical 
reflection will be largely determined by the presence of larger scatterers, result- 
ing in a sharper transition at the critical frequency. The two curves marked 
(c) and (d) in Fig. 6 have been calculated for groups of models of this type. 
As anticipated, this type of assembly leads to a considerably sharper transition 
than the assembly of similar scatterers and thence to the larger values of Y, 
just as required. It is not surprising to see that the curve representing assem- 
blies with a greater proportion of larger scatterers leads to the larger value 
of Y. The different points along the curve correspond to different values of p 
(see Fig. 7 of Part I for the properties of the size distribution). These curves 
are both drawn for assemblies with the same value for the critical frequency 
of the scatterers, 33.2 Mc/s. Increasing the critical frequency (up to the limit 
of 39.22 Mc/s) has the effect of raising these curves with very little shift in 
the X-direction. It appears likely that smaller values of aj, will give lower, 
flattened curves, covering the major part of the XY diagram. As explained 
in Part I, these curves have not been computed because they involve more 
effort, in the form of numerical integration, than the problem warrants. Thus 
any point in the region to the right of X =2 can be filled by any of an infinite 
number of assemblies having a certain combination of p, k,, and af. 

There are still no theoretical curves lying to the left of X =2. Since a value 
of Y greater than X is required, the reflection on the two higher frequencies 
must be due to weak scattering. Using this fact, a theoretical lower limit can 
be placed on X. An assembly of similar scatterers of size a=1 m gives a value 
for X of 1.44. Any smaller values of X could only be obtained with smaller 
scatterer sizes. If some of the scatterers in this assembly are replaced by 
larger and thus more strongly frequency-dependent scatterers, the frequency 
dependence of the assembly as a whole will inevitably increase, and with it 
the value of X. Thus, for physically reasonable sizes of scatterers, the value of 
X will be larger than 1.5 and usually larger than 2.0. Figure 3 shows that many 
echoes have values of X less than 2. These echoes can be readily accounted 
for by moving the curves in Fig. 6 downward to the left. As mentioned earlier, 
this is precisely the effect to be expected in the presence of absorption. This 
appears to be the correct explanation of the frequency dependence of such 
echoes. 

Although the results reported by Flood (1960) are not in a suitable form to 
examine them in the detail with which we have examined those of Forsyth and 
Vogan, it is useful to represent them on an XY diagram (Fig. 7). The marked 
dot indicates the over-all average values of X and Y obtained; the broken 
line indicates the values of X and Y inside of which fell all of the observations 
(excluding those which did not include echoes on the higher frequencies). 
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Fic. 7. XY diagram for the frequencies used by Flood (1960), illustrating his results. 
The curve gives values of X and Y for weak scattering from an assembly of similar gaussoids. 


The curve shown in Fig. 7 is for weak scattering, and represents theoretical 
values of X and Y from assemblies of similar scatterers having the sizes 
indicated. If the reflections are to be explained on the basis of this curve, the 
scatterers must have values of a as low as .35 or .40 m. Such sizes are equal 
to or less than the values of the mean free path and gyro-radius to be expected 
at heights of 100 or 110 km. Irregularities of such a small size are not likely 
to persist longer than the collision time, which is of the order of 1 microsecond. 
It therefore seems physically improbable that the reflection could be caused by 
such small irregularities. In addition, the position of the average values of 
X and Y lies significantly above the curve. There is no obvious way in which 
this observation can be explained in terms of weak scattering alone. Intro- 
ducing absorption will again shift the curve for weak scattering downwards to 
the left, making the explanation of the observations even more difficult. It 
appears likely that these results can be explained by the introduction of 
assemblies similar to those used to explain the results of Forsyth and Vogan, 
involving a distribution of sizes of scatterers and critical reflection at 49.7 
Mc/s. This is not to say that the observations of Flood contain no instances 
of weak scattering, but that the general nature of the frequency dependence 
suggests that critical reflection plays an important part in the explanation of 
the observations. The curves for such assemblies would be located at very 
large values of X, well off the region shown in the diagram. The effect of very 
strong absorption would be to bring the points down into the region of the 
XY diagram containing the experimental points. This interpretation is sup- 
ported by the observation by Flood (1959) that the strongest echoes on 49.7 
Mc/s are not accompanied by echoes on 226 Mc/s. Such an observation could 
be explained if 226 Mc/s echoes occurred only when the level of ionization 
was so intense that the 49.7 Mc/s echoes began to suffer considerable absorp- 
tion. The form of the results does not allow anything more specific to be said. 
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From time to time experimental observations which appear to be in agree- 
ment with a theory of weak scattering have been taken to imply that critical 
reflections are excluded from auroral reflections, and that the maximum 
electron density in auroral ionization is thus less than the critical density for 
the propagated waves: it is usually considered to depart only slightly from 
the ionospheric background. It is now clear that such observations do not 
exclude the presence of auroral irregularities having electron densities con- 
siderably in excess of the critical density for radio waves in the 30 to 50 Mc/s 
range. For example, the theoretical XY curves for an assembly of diffusing 
scatterers and for an assembly of similar, weakly scattering scatterers (Booker 
theory) are nearly identical for moderate values of X and Y (see Fig. 6). 
The first of these models includes strong critically reflecting scatterers, while 
the second includes none. 

Aspect Sensitivity 

Until now this discussion has not included any mention of the shape of 
auroral scatterers. In this regard, there is a great deal of evidence to support 
the existence of anisotropy of auroral.scatterers, but the actual extent of this 
anisotropy is by no means clearly established. As was pointed out in the 
discussion of the experimental evidence, both the azimuthal and elevation 
distributions of echo occurrence, particularly at high latitudes, and in spite 
of the complications already mentioned, definitely indicate the presence of 
aspect sensitivity. Also mentioned as evidence of aspect sensitivity, and hence 
of anisotropy of the scatterers, were the forwardscatter results of Collins and 
Forsyth (1959). Nevertheless a precise estimate of the anisotropy of auroral 
ionization is difficult to make. It is perhaps possible to put a rough upper 
limit on the degree of anisotropy to be expected by assuming that all the 
angular dependence of auroral echoes is due to the reflection mechanism. 
Presnell et a/. (1959) made this assumption to deduce (on the basis of weak 
scattering theory) from the maximum observed deviation angles from per- 
pendicularity to the magnetic field lines which they observed on 218, 398, and 
780 Mc/s, a “length” of 3.2 meters, corresponding in the notation of this 
paper to a value for ‘“‘ah’’. This implies, from the physical limitations on a, 
an axial ratio of 4 or 5. This result could be considerably in error, since the 
data used by Presnell et al. were not taken simultaneously on the different 
frequencies, nor with scaled systems. 

More recently, Greenhow, Neufeld, and Watkins (1960), using results 
obtained on 36 Mc/s with the Jodrell Bank radiotelescope, and using a theory 
of weak scattering, deduced a value for ah of about 26 meters. Their measure- 
ments are shown in Fig. 3. This gives a maximum value for h of about 30, 
considerably different from the value suggested by the results of Presnell et al. 
Greenhow et al. suggested that this discrepancy resulted from the lack of small 
aspect angles in Alaska, which would not give a proper base against which to 
judge the results. It is interesting to compare the results obtained by Greenhow 
et al. with theoretical curves based on the models of Part I. This has been done 
in Figs. 8a, 8b, and 8c, where the corrected results of Greenhow et al. are 
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compared with theoretical curves for weak scattering, ray theory, and weakly 
scattering assemblies with distributed sizes, respectively. It should be observed 
that the curves for the ray theory are uncorrected for the effect of the function 
5(¢, h, k’”) which was discussed in Part I. Preliminary studies of the nature 
of the function 6(¢, h, k’”) indicate that its inclusion would markedly change 
the shape of the curves in Fig. 7c, probably enough to move the A=5 curve 
down to or past the indicated position of the h=10 curve. It seems that all 
three theories are about equally capable of explaining the results of Greenhow 
et al. These experimental observations do not decrease as rapidly with in- 
creasing angle as any of the theoretical curves. The curves for the ray theory 
are the most similar in shape. The scattering medium may consist of a mixture 
of critically reflecting and weak scattering irregularities (quite possible, since 
these results are a statistical distribution over many different events), and the 
experimental points may not be satisfactorily corrected for such effects as the 
geographical distribution of the ionization and the volume of ionization 
intercepted by the beam as it looks in different directions. Both these effects 
may contribute to the observed difference between the experimental points 
and the theoretical curves. However ‘the results are interpreted, they imply 
a value for h between 5 and 10. 

On the basis of the discussion of the effect of the magnetic field on diffusion 
given in Part I, there is no reason to expect that the degree of anisotropy 
should be constant over the scattering volume (or, in view of the variation of 
the collision frequency with height, that the length of irregularities should be 
constant throughout the scattering volume). This will depend, of course, on 
the range of heights over which the echoes are obtained, but there appears 
to be no reason why the scatterers at, say, 100 km, could not be isotropic, 
while those at 110 or 115 km might be very much elongated along the field 
lines. The bistatic observations of Collins and Forsyth (1959) suggest two or 
three different types of aspect sensitivity. 


CONCLUSIONS 


In spite of the relatively few experimental observations which are relevant 
to the determination of the characteristics of auroral ionization, it has been 
possible to arrive at several important conclusions concerning this ionization. 
The existence of critical reflection at frequencies between 30 and 50 Mc/s 
appears to be established. This implies electron densities which at times reach 
values as large as 3X10" electrons/m*. This value does not seem to be in- 
compatible with densities deduced on the basis of the intensity of auroral light. 
There is unquestionable evidence that the scatterers composing the auroral 
ionization are elongated along the lines of the magnetic field, but the extent 
of this elongation is by no means settled. Although there is no reason to expect 
that the degree of anisotropy is constant throughout the auroral ionization, 
the characteristic axial ratio of auroral scatterers appears to be between 5 and 
10. The scatterers which are important for the reflection of radio waves have 
sizes of the order of a few meters (by size is meant that distance measured 
from the axis of the scatterer, in a plane normal to it, at which the electron 
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density is 1/e of its maximum value), the experimental results being consistent 
with a distribution of sizes as given by (41) of Part I, having the values of 
ad range from 2 to 5 m or more, and those of p range from 1 to 10. 

It is clear from the nature of the experimental observations used to obtain 
them, that the characteristics given above for auroral ionization can not be 
considered to be final. The purpose of this paper has been simply to indicate 
a broader point of view in considering the possibilities for an auroral model, 
and in so doing to indicate the lack of adequate observational evidence on 
which to base a model of radio-aurora. It is evident, for example, that none 
of the observations made to date are really suitable for studying the aspect 
sensitivity of the reflection mechanism. What is needed is an experimental 
arrangement which will measure the strength of the signal reflected from the 
same volume at different angles. Similarly, the frequency dependence of auroral 
reflections will not be free of ambiguities until scaled systems are operating 
on three or more frequencies simultaneously. It appears that it would be 
particularly valuable if these frequencies were chosen within about 10 mega- 
cycles of each other in the frequency range of 50 to 100 Mc/s. Because of the 
nature of the differences between the predicted frequency dependence of the 
various auroral models, as well as a superimposed effect due to absorption, 
distinctions are better made with closely spaced frequencies than widely 
separated ones. 

It is hoped that this paper has also illustrated the usefulness of the particular 


approach to models of ionization used in Part I. This approach allows inclusion 
in one treatment of both critical reflection and weak scattering and enables 
one to discuss the experimental results in terms of what seem to the author to 
be physically significant parameters. 
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SCATTERING OF RADIO WAVES BY AN IONIZED GAS IN 
THERMAL EQUILIBRIUM IN THE PRESENCE OF 
A UNIFORM MAGNETIC FIELD! 







J. A. FEJER 







ABSTRACT 


In an earlier paper (Fejer 1960) a theory was developed for the scattering of 
radio waves by the electron density fluctuations that exist in an ionized gas 
in thermal equilibrium. The theory treated only the extreme cases where the 
“‘characteristic scale’ of the scattering irregularities is either very much larger 
or very much smaller than the Debye length. The presence of only one type of 
singly charged ion was considered and the ion and electron temperatures were 
assumed equal. The effects of an external magnetic field were not taken into 
account. 

These earlier limitations are removed in the present paper and the effects of 
an external magnetic field are taken into account. 

It is shown that the total power is independent of the magnetic field and an 
expression for the frequency spectrum of scattered power in the presence of a 
uniform magnetic field is obtained. Useful approximations to this expression 
are derived for various limiting cases of interest. 

It is concluded that the magnetic field need not be taken into account in the 
interpretation of past observations by Bowles (1958, 1959) and by Pineo, Kraft, 
and Briscoe (1960). In future experiments, however, particularly at great heights, 
the effect of the magnetic field could be considerable. 



























1. INTRODUCTION 


Gordon (1958) suggested that radio waves incident on the ionosphere at 
frequencies well above thecritical frequency are scattered by irregularities in 
electron density which are present even under conditions of thermal equili- 
brium. His simple theory predicted that the scattered power would be pro- 
portional to the electron number density, and that the width of the frequency 
spectrum would be determined by a Doppler spread caused by the thermal 
velocities of the electrons. 

Observations of this phenomenon by Bowles (1958, 1959) and by Pineo, 
Kraft, and Briscoe (1960) have approximately confirmed Gordon's predictions 
about the scattered power but showed the frequency spectrum to be much 
narrower than anticipated. 

Bowles (1958) has explained this discrepancy by showing that, under the 
conditions of the experiments so far performed, the Doppler spread is deter- 
mined by the thermal velocities of the ions although the electrons are respon- 
sible for the scattering. Subsequently several authors (Fejer 1960; Dougherty 
and Farley 1960a; Renau 1960; and Salpeter 1960a, 6) have developed the 
theory of this phenomenon, and using a variety of methods have arrived at 
identical conclusions. The aim of the present paper is first to remove certain 
limitations from one of these developments (Fejer 1960), and then to extend 
it to include the effects of the magnetic field. 

One of the most serious limitations of the earlier treatment (Fejer 1960) 
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was the assumption that the Debye length was either very much greater or 
very much smaller than the ‘‘characteristic scale’ of the scattering irregu- 
larities. In removing this limitation and generalizing the treatment to several 
types of ions and to unequal values of the ion and electron temperature no 
fundamentally new results are obtained which were not already present, 
implicitly or explicitly, in two of the earlier treatments (Dougherty and Farley 
1960a; Salpeter 19606). There are, however, two reasons for devoting Sections 
2-5 to the generalization of Fejer’s earlier work. The first reason is that the 
generalized theory is needed for the treatment of magnetic effects in later 
sections. The second no less important reason is that the new treatment makes 
fewer assumptions and is therefore simpler and of a more general applicability. 

The basic assumptions are that in the absence of Coulomb forces, the 
particles are randomly distributed, and that their velocity distribution is 
Maxwellian. (In principle a different velocity distribution could be assumed 
and would then lead to different results.) For purposes of the present paper 
the particles are said to be randomly distributed if their probability density 
in configuration space is independent of the co-ordinates of all the particles. 

It is further assumed that the electric field acting on a particle is, for the 
present purposes, well approximated by the Coulomb field of the other par- 
ticles. It is difficult to justify this assumption fully without including the 
effect of induced electric fields in the theory first and then showing that their 
effect on the result is small. It is easily shown that, if the random thermal 
motion were unperturbed by forces acting on the electrons and ions then the 
electric field caused by the thermal motion would be very much smaller than 
the Coulomb field for thermal velocities much smaller than the velocity of 
light. It is further possible to show that for the parameters of interest the 
electric fields, induced by the magnetic field of currents due to the unper- 
turbed random thermal motion, are insufficient to perturb that motion appre- 
ciably. Such considerations make it seem likely that the perturbation of the 
random thermal motion and the resulting electron densities are mainly deter- 
mined by the Coulomb forces. 

In Sections 6-8, the effects of a uniform magnetic field are taken into 
account in the calculation of the spectral power density. The total scattered 
power could be calculated from the results of Section 7 by integrating the 
spectral power density over the whole frequency spectrum. An alternative 
method of calculation for the total scattered power was given in the Appendix 
of the earlier treatment (Fejer 1960). That treatment is not altered by the 
presence of a magnetic field if the electric field energy is taken to be only 
that contained in the Coulomb field. Therefore the results derived there for 
the total scattered power are not affected by the magnetic field within the 
limits of the present approximation. The author would, however, like to 
take this opportunity to point out that the result derived in that paper for 
the total scattered power was in error. The nature and cause of this error 
was explained in an Erratum in the September issue of Canadian Journal of 
Physics. The correct form of Fejer’s (1960) equation (8) (also derived by 
Renau (1960) and by Salpeter (19600)), is 
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[4alp sin (8/2)]’+2? 
[47lp sin (6/2) ]?+2n? 


as already given in the above-mentioned Erratum. 

It is of interest to mention that Dougherty and Farley (19606) have also 
extended their earlier theory (Dougherty and Farley 1960a) to take the effect 
of the magnetic field into account. They first derive an exact expression for 
the spectral power density without approximating the forces acting on the 
electrons by the Coulomb forces. They then show by contour integration of 
this expression that the total scattered power is not affected by the magnetic 
field and is given by an expression identical with the above corrected form 
of Fejer’s equation (8). In subsequent parts of their paper, Dougherty and 
Farley (19605) make the same approximations as the present paper in order 
to simplify their rather cumbersome exact expressions for the spectral power 
density. They derive a sufficient condition for the validity of this approxi- 
mation but do not show that the condition is also a necessary one. Their 
condition is satisfied in most ionospheric applications and thus their treatment 
provides further justification for the approximations made here. Results 
similar to those of the present paper were also obtained separately by Renau, 
Salpeter, and Hagfors (private communication). 

The result for the spectral power density is given in two forms in Section 7, 
as an integral and as an expansion into tabulated functions. Either of these 
forms could be used for a computation of the frequency spectrum. 

A better insight into the effects of the magnetic field is gained by considering 
certain extreme cases where the results can be expressed more simply. Two 
such extreme cases are considered in Section 8, where it is assumed that the 
gyro-radius is either very much greater or very much smaller than the ‘‘charac- 
teristic scale’’ of the scattering irregularities. The approximations obtained 
for these two extreme cases are probably sufficiently accurate for most iono- 
spheric applications and in addition provide a good insight into the nature 
of the spectrum even in intermediate cases. Computations would be needed 
for more accurate results. 


o =o0,V 


2. THE TOTAL SCATTERED POWER FOR RANDOMLY 
DISTRIBUTED ELECTRONS 
It may be shown (for example, from equation (4) of Fejer 1960) that for an 
ionized gas, for sufficiently high frequencies, for weak scattering by irregu- 
larities in electron density, and for variations of electron density in time 
which have dominating periodic components much longer than the period 
of the incident radio wave 





(1) o=0,V 


J AN,(r,t) exp[i(k,—ko).r]d*r 





where o, the scattering cross section of the gas, is defined as the power per 
unit incident power density, unit solid angle, and unit scattering volume, and 
where V is the scattering volume, A.V, is the deviation of the electron number 
density from its average value V,, where the suffix 1 denotes electrons, Ko is 
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the propagation vector of the incident wave, k, is a vector of the same mag- 
nitude as kp which points in the direction of scattering, r is the radius vector 
to a point in the scattering volume, d*r is the volume element in the scattering 
volume, and o, = (ue? sin y/4rm,)? is the Thompson scattering coefficient 
of an electron. In the latter expression, uo is the permeability of free space, 
—e is the electronic charge, y is the angle between the electric field vector 
of the incident wave and the direction of scattering, and my, is the electronic 
mass. The bar in equation (1) indicates either averaging in time or an ensemble 
average. It is convenient to introduce the vector k = k,—Kp. 

It is relatively simple to show (Kahn 1959) that for randomly distributed 
stationary particles 


ee 
iS AN, exp(ik.r)d’*r| = NiV 
V 


where the bar now indicates an ensemble average. Substitution of (2) into 
(1) yields 


(3) ¢= Nice. 


Equation (8) expresses the result originally stated and proved by Rayleigh 
(1871) that the power scattered by m random scatterers in any direction 
other than the direction of the incident wave is on the average u times the 
power scattered by a single scatterer. Rayleigh’s proof is based on the fact 
that the fields due to randomly distributed scatterers have random phases. 

The above step leading from equation (2) to equation (3) could be reversed 
in order to derive the mean square value of the spatial Fourier component 
of electron density from a knowledge of the scattering cross section. Such a 
reversed procedure is useful for some purposes and is followed in Sections 3 
and 4. In Section 3 the frequency spectrum of scattered power for randomly 
distributed electrons is derived and exemplified by application to a Max- 
wellian gas. In Section 4 the space-time Fourier component of the number 
density for random scatterers is then derived from the frequency spectrum 
of scattering. 


3. THE FREQUENCY SPECTRUM OF SCATTERED POWER FOR RANDOMLY 
DISTRIBUTED ELECTRONS IN THE ABSENCE OF COLLISIONS 

The frequency spectrum of the scattered power may be calculated by 
dividing the electrons into groups which have almost the same velocity com- 
ponent in the k direction and therefore produce almost the same Doppler 
shift. The scattered power in each group is given by equation (3), with NV, 
replaced by the mean number density of electrons within the group. If f(z) 
is the normalized velocity-component distribution in the direction of the k 
vector, then this mean number density is N,f(v)dv for the group that has 
velocity components between v and v+dv. The Doppler shift for this group 
of scattering electrons is 
(4) Ws—wo = kv 


where wo is the angular frequency of the incident wave, and w, the angular 
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frequency of the scattered wave. Abbreviating w = w,—wo, one obtains for the 
scattering cross section o,dw of this group with the aid of equation (3) 


(5) owiw = N,f(v)dve, 

where o, satisfies 

(6) o = foudw. 

Equations (4) and (5) yield the expression 

(7) oe = Nif(w/k)oe/k 

for the spectrum function oj For a Maxwellian distribution 
(8) f(v) = (m,/2nKT))' exp(—my?/2KT;) 


where 7) is the temperature of the electrons and K is Boltzmann’s constant, 
and then 


(9) a» = (Ni0./k) (mi/2rKT;)* exp(—myw?/2K Tyk?). 


4. THE SPACE-TIME FOURIER COMPONENT FOR RANDOM SCATTERERS 


The number density deviation AN from the mean may always be expressed 
as a Fourier integral in four dimensions 


(10) AN (r,t) = fs(w,k) exp[i(wt—k.r) ]dwd*k. 


A simple and general relationship between the space-time Fourier com- 
ponent, defined by (10), and the frequency spectrum o, will be derived in 
this section. 

For the mathematical convenience of keeping s(w,k) finite the deviations 
AN are assumed to vanish outside a finite time interval to which can be made 
arbitrarily large compared with the period of interest. Similarly the deviations 
AN are assumed to vanish outside a finite volume V. 

In equation (10) the deviations AN are decomposed into travelling density 
waves with different frequencies w and wave numbers k. Density waves of 
angular frequency w (and the correct wave number) are then responsible for 
scattered radio waves of angular frequency wo+w. 

The inverse form of equation (10) with respect to k is 


(11) J? sok) exp(twt)dw = (any f AN exp(ik.r)d’r. 
aoe Vv 


Equations (10) and (11) are valid for any type of particle. If the particles 
are electrons and if only those density waves, whose angular frequency lies in 
the small but finite interval between w and w+Aw are considered, then the 
density deviation A.V; corresponding to this group of density waves satisfies 


w+ Aw 
(12) J $1(w,k) exp(iwt)dw = eny*f AN, exp(ik.r)d’r . 
V 


@ 


This group of density waves is then responsible for scattered radio waves 
with angular frequencies between wo+w and wo+w+Aw. The power scattered 
by these waves is then from equation (1) 
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w+ Aw So) ok eee d aie 
(13) f oddw = o,V* f AN, exp(ik.r)d’r 
o Vv 


or by combination of (12) and (13) 


2 
= oAw 








2 


w+ Aw 
J $1(w,k) exp (iwt)dw 








(14) o,Aw = o,V‘(27)° 


where the bar indicates averaging in time over the time interval between 
t = 0 and ¢ = to which is approximately the time interval over which the 


function u(t) defined by 
w+ Aw 
(15) u(t) = ar) f $1(w,k) exp (twt)dw 


has non-vanishing values, if to >> (Aw)~!. If equation (15) is satisfied then by 
a well-known theorem on Fourier transforms 


oo a“ w+ Aw es 5s 
J |u|“dt = J \s1(w,k) |"dw & tolu|* 
and a combination of this equation with (14) leads to 
(16) = (29)'oeV~ to '|s1(w,k) |? 
where 


ne Al w+ Aw 
|s1(w,k) |* = (a) f |s1(w,k) | "dw. 


The bar in equation (16) must therefore be interpreted as an average over 
a small but finite angular frequency range Aw in the vicinity of w. The quantity 
|s1(w,k)|? varies rapidly and irregularly with frequency in the interval between 
w and w+Aw but its average value is well defined. 

The mean square value |s,;|? of the space-time Fourier component of ran- 
domly distributed electrons may now be obtained from a combination of 


equations (9) and (16) 
(17) |se1(w,k)|? = (24) ~"VtoN R*(m1/22KT})! exp(—mywo?/2k°KT}) 





where the mean value of |s,1(w,k)|? is formed over a small frequency range 
which, however, is large compared to f;'. Expression (17) is by no means 
restricted to electrons and is equally valid for the Fourier component of the 
number density for any type of particle as long as the velocity distribution 
is Maxwellian and the positions are random. An alternative direct derivation 
of an equation similar to (17), which does not use the concepts of scattering 
and Doppler shifts, is given by Dougherty and Farley (1960a) in their 
Appendix. 

The subscript 7 in s,; has been used to distinguish this space-time Fourier 
component, which would exist if the electrons were randomly distributed in 
space, from the actually existing space-time Fourier components 5s}. 
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5. FREQUENCY SPECTRUM OF SCATTERED POWER FOR A 
COLLISION-FREE PLASMA 

In the derivation of equation (17) it was assumed that the particles move 
in straight lines at constant velocity. This can be true only if there are no 
forces acting on the particles. The particles of a plasma are charged, however, 
and there are electrostatic forces acting on them. No external magnetic field 
is assumed in this section and the magnetic field of the particles themselves 
is neglected. The electrostatic forces modify the motion of the particles and 
therefore also modify the spectrum of scattered power. Fejer (1960) has 
shown how this modification to the frequency spectrum of scattering can be 
introduced by a correction factor. His treatment, which was restricted to 
very small values of the Debye length, is generalized here. 

Let the mean number density of electrons in the plasma be ,, and let 
the charge of the jth type ion be Z,e, and its mean number density 1, so 
that 


(18) 2 Z,;N;=0 
j=l 


where Z; = —1 (for electrons). 

If, in the absence of electrostatic forces, the particles were randomly dis- 
tributed with a Maxwellian velocity distribution, the mean square space-time 
Fourier component of the electron number density would be given by (17). In 
general, the mean square space-time Fourier component of the number density 
for the jth type ion would be 


(19) |5,)(@,K) > = (20)~’VtoN km ;/24KT;)* exp(—mywo"/2kKT)). 


The real problem, however, is the determination of the mean square space- 
time Fourier component |s;(w,k)|? for electrons in the presence of the electro- 
static forces, having obtained a knowledge of |s,;(w,k)|? for the jth type of 
charged particle in the absence of electrostatic forces. If this more realistic 
quantity |s,(#,k)|? were obtained, the spectrum function of scattering 4 
would be given by equation (16). 

The meaning of the quantities s is made clear by equation (10); they are 
associated with elementary number density deviations of the form s(w,k) 
exp[i(wt —k.r)]dwd*k. In the following calculation, the quantities s are taken 
to be these elementary number density deviations, although the factors 
exp[i(wt—k.r)]dwd*k are omitted. The quantities s,1, sp2,..., Sr are given 
by equation (19), and represent number density deviations in the absence of 
electrostatic forces. A new set of quantities S41, 512,..., St are introduced 
and are taken to represent the elementary number density deviations which 
would arise if the electrostatic forces acting on the particles were taken to 
act on a uniform distribution of particles. The elementary number density 
deviations 51, 52,...,5, actually present may then be equated, within the 
accuracy of a linear approximation, to the sum of the deviations s,; which 
would have occurred in the absence of the fluctuating electric field, and the 
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deviations s,;; which would have been caused by the fluctuating electric field 
if it had acted on a uniform distribution of particles. Thus: 


(20) Sy = SpgHStj. 
The fluctuating electric field has the charge density 
e b Z 583 
j=l 


as its source, and so is proportional to this quantity. Similarly the quantities 
Stj must be proportional to the electric field and therefore may be expressed 


(21) Sty = §; x Z Sq 
ee 


where the é are factors of proportionality which remain to be calculated. 
Multiplication of equations (20) by Z; and summation over j results after 
substitution of s;; from (21) in 


(22) L158; = 2X, Zisrt( 3 ze) Ze.) 


or 


(23) dX Z 483 = (= Zs) J (- X Z4t)). 


Therefore, from (20), (21), and (23) 


1-2 2; nh Of 
(24) iy a teh Bila 
j= 
ee oe 2 
j= = 


Since the s,; in equation (24) are statistically independent and of random 
phase 

| n 2 

es 2. eae E Oe 

(25) |ss[° = |—__=— sal + |———— a Z3\seil's 


7% Zs;| }1- D a 


This is the required general equation for |s;(w,k)|?. It is equally valid in the 
presence or in the absence of the magnetic field. 

For the present section, then, there remains only the calculation of the 
factors of proportionality £; in the absence of an external magnetic field. In 
this calculation one uses Liouville’s equation 


(26) dW’/dt+v.grad,W’+F,.grad,W’ = 0 
7-8 


| 
| 


where W’ is the distribution function of the particle type considered, and 
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F,, is the force per unit mass. The direction of the x-axis will be taken parallel 
to k for convenience and a force per unit mass proportional to exp[i(wt—k.r)] 
will be assumed. 

The electric field E responsible for the force per unit mass F; = Z,eE/m, 
has its source in the charge density 


n 
‘7. aay 
j=1 


and therefore 
(27) eo divE = —ikeE =e >> Z3s; 
j=1 


where €o is the dielectric constant of free space. 
If the velocity distribution of the unperturbed gas is assumed Maxwellian, 
one may write 
(28) W’ = (1+W) (m,/2rKT;)3” exp(—myv?/2KT,) 
where W is a perturbation which is a function of velocity, position, and time. 
If (28) is substituted into (26), one obtains after linearization, 
(29) OW /dt+v,0W/dx— (m;/KT;)vzF; = (). 


Substituting the force per unit mass from (27) into (29) gives 


(30) iwW —ikv,W+(my,Z je"/KT jikeom;) >, Z;s; = 0. 
j=1 
Solving (30) for W 


n 
Ze Z 245; Ae 
j=1 a 
RKTj¢9 w—kv," 
The perturbation (m,/2rKT,)*/?W exp(—myv?/2KT;,) in the distribution 
function W’ must be integrated over velocity space to obtain the fractional 
perturbation s;;/N, in the number density. Using (31) this leads to 


( m, Se x a ( my ‘ 
9 seit id cinch ne gigeas cea re ; 
(63) Sta ™ \OEKT, RKT 40 a ON oT 

The integrand in (32) has a pole at v, = w/k. The path of integration near 
the pole can be decided by letting w have a very small negative imaginary 
part so that the waves exp[i(wt—k.r)] represent waves growing in time 
which at ¢ = — © had negligible amplitude. (The more rigorous method of 
using Laplace transforms has been described by Landau (1946).) 

The integral in (82) may then be replaced by the principal value minus 
mi times the residue of the pole. This leads to 


(33) sty = z( . Z4s4) We" /eaKT)| -1+2(0/2,) exp(—w"/2,") 


(31) W= 


w/Q; 
J exp y'dy+in'(w/2;) exp(—o'/a | 
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where 
(34) Q, = k(2KT,/m,)}. 








Combination of (33) with (21) gives, 





w/Q; 
(35) & = Z| -1+2(0/2,) exp(—w"/9)) J exp y'dytin' (w/2,) 





exp(—w"/ a | 





where 
(36) lp; = (K T €0/ N ye?)* 


is the Debye length for the jth particle. 
Combination of (16), (17), (25), and (35) yields for the scattering cross 
section 

















n 2 
, Z 33 £ 2 
«= |—= 1 40;Nyo, exp(—w*/0}) + | —# — 
~ X Z 5; Soe X Z sk; 
j= j= 





Zz nr *03'N jo5 exp(—w’/9}). 


j=2 






In the absence of the magnetic field, equation (37) is the formal solution to 
the problem. The quantity o, represents the power scattered per unit incident 
power density, unit solid angle, unit scattering volume, and unit width of 
angular frequency. 

Examination of equation (37) shows that the first term on the right repre- 
sents a relatively wide spectral distribution of scattered power; the scattered 
power contained in the succeeding terms is spread over a much narrower 
range of frequencies. It is convenient to refer to the first term in equation (37) 
as the electronic term and to the succeeding terms as the ionic terms. This 
terminology will also be used for the terms of equation (25) and will be more 
fully justified in Section 8. 

In the special case nm = 2, Z; = —1, Z. = 1, T; = T2 equation (37) agrees 
with the corresponding equations of Dougherty and Farley (1960). Similarly 
Salpeter’s (1960) and Fejer’s (1960) results can also be derived from equations 
(35) and (37). 

Since certain special cases of this equation have already been discussed in 
some detail in the above papers the discussion here will be restricted to some 
additional computed results shown by Figs. 1-3. 

Figure 1 refers to the case of m = 2, Z; = —1, Z2 = 1, 7; = T2, m, K mz. 
The curves on the right are identical with the curves of Salpeter’s (19606) 
Fig. 2. They represent o,92;/N as a function of w/Q, for different values of the 
parameter a = (kp)! for large values of w. The curves on the left represent 
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71G. 1. The frequency spectrum for different values of the parameter a@ = f-/=1 T 
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positive ions are assumed to be atomic oxygen. The ion and electron temperatures are assumed 
equal. 
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Fic. 2. The frequency spectrum for very large values of the parameter a = k~/p! and for 
different ratios T./T; of the electron temperature to the ion temperature. 
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Fic. 3. The frequency spectrum for a plasma containing both atomic oxygen and hydrogen 
ions for different values of the ratio of their number densities No and Ng, for large values of 
the parameter a. 


o422/N as a function of w/Q: and thus represent the same curves drawn 
on a different scale for a different frequency range in order to show the 
behavior for small values of w. 

These two sets of curves of Fig. 1, show that for small values of the para- 
meter a the spectrum function is very broad (the electronic part of the spectrum 
predominates) and of Gaussian shape; for large values of a the spectrum 
becomes narrow (the ionic part of the spectrum predominates) with the 
exception of a small amount of power scattered at frequencies differing from 
the frequency of the incident wave by the plasma frequency. 

Figure 2 shows the spectrum functions o,Q2/.V of a two-component plasma 
for a> 1 for different ratios 7,;/T2 = T./T; of the electron and ion tem- 
peratures as a function of w/Q2. It is seen that for low values of this ratio a 
Gaussian but narrow spectrum is obtained as predicted by Bowles (1959). 
The curve for 7, = 72 has been discussed previously. For large values of 
T,/T2 the total scattered power is greatly reduced and the spectrum has 
two sharp peaks which correspond to scattering by acoustic waves in the 
plasma involving the movement of both ions and electrons (Greene 1959). 
The Landau damping of these acoustic waves is low when the electron tem- 
perature is much greater than the ion temperature. 

The curves of Fig. 2 show that given a knowledge of the ionic mass, the 
shape of the spectrum function not only gives an indication of the electron 
temperature, but also of its ratio to the temperature of the positive ions. For 
example, observations by Pineo ef al. (1960), interpreted with the aid of the 
present theory, tend to show that the temperature of electrons is approximately 
equal to the temperature of positive ions in the F region and immediately 
above it. 

Figure 3 is drawn for a three-component plasma consisting of electrons 
and positive atomic oxygen and hydrogen ions for a> 1, 7; = T; = T;3 for 
different values of the number density ratio V3/N2 = Ny/No of hydrogen 
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to oxygen ions. Figure 3 thus shows the transition from a narrow spectrum 
corresponding to oxygen ions to a somewhat broader spectrum corresponding 
to hydrogen ions. 


6. FREQUENCY SPECTRUM OF SCATTERED POWER FOR RANDOMLY 
DISTRIBUTED ELECTRONS IN A UNIFORM MAGNETIC 
FIELD IN THE ABSENCE OF COLLISIONS 

The calculation carried out here uses a technique similar to that employed 
by Fejer (1960) in the calculation of the frequency spectrum for randomly 
distributed electrons with a high electron-neutral collision frequency. 

The distribution of displacements /, after time ¢ in a direction parallel to 
the wave number vector k = k,—Ko must be determined first. A new Car- 
tesian co-ordinate system is introduced for this purpose with the z axis in the 
direction of the magnetic field and the k vector in the x,z plane. Consider 
the motion of an electron which, starting from the origin, reaches the point 
X14,¥ 12, after time ¢. The distance R, = (x?2+-y2)}, calculated from the equations 
of motion is then 


(38) R, = (w/Qz)[2(1—cos Qzt)]} 


where Qy is the electron gyro-frequency (the suffix 1 for electrons is dropped 
in this section) and w is the magnitude of the initial velocity component 
perpendicular to the magnetic field. The distance z, travelled in the z direction 
in time ¢ is given by 


(39) 3, = ut 


where u is the velocity component in the z direction. 
The Maxwellian velocity distribution, written in terms of the velocity com- 
ponents wu and w, is given by 


(40) df = (m/2xKT)*” exp[—(m/2KT) (u?+w?) ]2rwdwdu 


or rewritten with the aid of equations (38) and (39) in terms of x,, 4, 2: 


m \*? Qn 
(41) df= (55) 2(1—cos Qut) 


it Saito) |t 
exp stn | H+ tee Out) AF edede. 


Let the unit vector in the direction of k have the Cartesian components 
sin 6, 0, cos 6. Then the distance /, in the direction k is given by 


(42) l, = x,sin 0+2, cos 6. 


The distribution may now be expressed in terms of x, y,, and J;: 


Mr 





(43) df ( oF seiin Mca — 
Y  \OaeKT 2(1—cos ryt) cos O° PC oni t cos 6 


Gilet V dea 
31 —cos Qué) dx, dy, dl,. 
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The distribution in terms of /, alone is obtained by carrying out the inte- 
grations in x, and y,. The result, without the ome factor, is 


(44) f« exo] — 5 QKT f° cos’ 6+203%"(1—cos Qyt) sin“d dl. 


It is assumed with Fejer (1960) that the particles moved with hypothetical 
velocities v such that vf = /, and it is further assumed that Doppler shifts 
@s—wo = w = kv result from these velocities so that the frequency spectrum 
of the scattered waves is proportional to 





(45) ex = (ws—wo) *1 */k? } 
“Pl “OKT Fcos “+204 (1—cos Qyt) sin’6 
The Fourier transform of this hypothetical spectrum gives a hypothetical 


autocorrelation function 


(46) pi(r) = exp) —iaor—1" At [t°cos"0+20%°(1—cos Qyt) cin’e)} ‘ 
The value of p;(7) is applicable to the actual problem only for 7 = ¢ 


(47) p(t) = exp{ —twot— 30" [t’ cos’ 0+20%"(1—cos Qyt) sin*6]} 
where Q is defined by equation (34). Then in terms of the Wiener—Kintchine 
theorem (Wiener 1930) the relative power spectrum Q(w,) is given by 


Olen) = (2n)* fp) expling) 
where 
J QGosrdeng = 1 


The quantity «, is then obtained in terms of equation (3) from oj = Q(ws) Noe 
where JN is the electron number density. Thus 


(48) o= (No,/2n) f exp {iwt — 32° [t* cos” 0+20%"(1—cos Qyt) sin? 6]} dt 
Comparison of this expression for the power spectrum o, with equation (16) 


vields for the mean square space-time Fourier component |s,|* of randomly 
distributed electrons with a Maxwellian velocity distribution 





(49) [s-(w,k) |? = (2e)*VuN f exp {iwt — 30°[t’cos’ 0 


+207" (1—cos Qyt) sin” 6]} dt. 


This is one of the quantities which have to be determined before equations 
(16) and (25) are applied. Equation (49) may be rewritten to give |s,,]? for 
any particle 7 by attaching the suffix 7 to V,, and Qg on the right. 

It may be convenient sometimes to expand the factor exp(}2?Q;? sin? 6 
cos Qgt) in the integrand of (49) with the aid of the identity (Watson 1948) 
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(50) exp(x cosa) = Io(x)+2J;(x) cosa+2I2(x) cos 2a+.... 


where J,(x) is the modified Bessel function of order p. The integration indi- 
cated in (49) can then be carried out, leading to 


4 — var = 2 ee 
. 7 = (2 i a - ——F ————— 
(51) Ise(wk)|" = (24) TaaQ al exol 202 ut 202 


p=—c 


exp] — et 609)* | 


Q*cos” 6 


The spectrum o, or the corresponding quantity |s,(w,k)|? may then be 
regarded as the sum of a number of Gaussian functions of w, each centered 
on one of the harmonics (positive or negative) of the gyro-frequency of the 
particle. These Gaussian functions degenerate into 6 functions when the wave 
number vector is perpendicular to the magnetic field. 

An equivalent result has also been obtained by Laaspere (1960), who used 
a different method of calculation. 

It follows from equation (49) that for Q cos 0/Qq > 1 the integrand vanishes 
for large values of Qgt. For small values of Qgt one may replace 1—cos Qgt 
by 302t? and therefore the effect of the magnetic field can be neglected if 
2 cos 0/Qy > 1. 


7. CALCULATION OF THE FACTORS ¢; IN THE PRESENCE 
OF THE MAGNETIC FIELD 

In the previous section the quantities |s,,;|? have been determined in the 
presence of the magnetic field. They represent the mean square space-time 
Fourier components which would exist in the absence of Coulomb inter- 
actions. Before equations (16) and (25) can be applied to the calculation of 
the spectrum, the factors £; defined by equation (21) have to be determined. 
These factors have already been calculated in Section 5 for the field-free case 
from Liouville’s equation. The present extension of the calculation closely 
follows Bernstein’s (1958) analysis. In the presence of the magnetic field of 
induction B an additional force term (Ze/m)(vXB) has to be included in 
Liouville’s equation which takes the form 


(52) OW’/dt+v.grad,W’+ (Ze/m)(E+vXB).grad,W’ = 0 


where the subscripts j are again omitted. Substitution of expression (28) for 
W’ vields for the perturbation W after linearization 


(53) iwW —i(k.v) W+ (Ze/m) (vX Bo). grad,W— (Ze/KT)(E.v) = 0 


where By is the induction of the uniform dominating magnetic field. Since 
the electric field E is taken to be the Coulomb field, it must be parallel to 
the vector k, which points in the direction of the wave normal of the appro- 
priate electron density waves. The electric field vector may therefore be 
expressed as E = Ek/k where k/k has the components sin 6, 0, cos 6. After 
introduction of the cylindrical velocity co-ordinates w, ¢, u equation (53) may 
be written in the form 
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(54) dW/d¢+0%'[—iw+ik(w sin 0 cos ¢+u cos 6)]W 
+(E/Bo)(m/KT)(w sin 6 cos ¢+u cos 6) = 0 
where Qy = ZeBy/m. If w is again assumed to have a small negative imaginary 


part, the required solution of this first-order linear differential equation, which 
must be periodic in @ with a period of 27, is 


(55) W = —(me/KTBo) exp{Qy'[(iw—iku cos 0)¢—ikw sin @ sin ¢]} 


> 
(w sin 6 cos ¢’+u cos 6) exp{Qy'[(—iw+iku cos 0) ¢’+ikw sin @ sin ¢’]}d¢’. 
steo 


The sign of the lower limit of the integral is taken to be the same as the 
sign of the charge. After introduction of the new variable 8 = (¢’—@)/2 in 
the place of ¢’ and of the notation 8-+@¢ = a equation (55) may be rewritten 
in the form 

too 
(56) W = (2mE/KTB,) J [w sin 6(cos a cos B—sin a sin B)+u cos 6] 

p= 


exp{Qz'[(—iw+iku cos 0)28+2ikw sin @ cos a sin B]} dB 


where the sign of the upper limit is taken to be the same as the sign of the 


charge. 
In this form the periodicity of the solution (55) for W is clearly seen. In 
view of the discussion following equation (31) 


(57) sp = N(m/2eKT)*” J wiw f du { “W expl—m(u?+w*)/2K TMs 
0 —c 0 


where W is given by equation (56). The expression on the right of equation 
(57) is then an integral in the four variables w, u, ¢, 8. It is convenient to 
start with the integration over ¢ and to change from ¢ to the new variable 
a = ¢+6 already introduced in equation (56), with a corresponding change 
in the limits of integration. Further in terms of equations (21) and (29) one 
can write 


(58) E = tes;/RkeoE 


and combination of this equation with equation (57) results in 


toa oo oo 
(59) & = (2QieN/keoBo)(m/KT)(m/22KT) wf dB f wdw | du 
0 0 —ow 


B+ 


exp[—m (u*+w’) /2KT+280q'(—iw+iku cos 0)] (w sin 8 cos B cosa 
8 


+u cos @) exp (2i2% kw sin 6 sin 6 cos a) da. 


The integrations over u are elementary. The integrations over a@ are readily 
carried out, using (Watson 1948) 


Qn 


: exp(ix cos ¢) cos(p¢) d¢. 


(60) Iyle) = 2a? 
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The result is 


sto © 
(61) & = (ieN/keoBo)(m/KT)? 7 dp J w exp(—mw’"/2KT —2iBw/Qy 





—2KTB°R’ cos’ 6/mQz) [2iw sin 6 cos BJ;(2kw sin 6 sin B/Qyz) 






+4i0%'Bk cos’0(KT/m) Jo(2kw sin 6 sin B/Q,)] dw. 


The integration over w may be carried out, using (Watson 1948) 







(62) Joes (aw) exp(—w’p’)dw = a’(2p’)~”* exp(—a’/4p’). 





The result is, using the notation of equations (34)—(36) without the sub- 
script j 






too 
(63) &= — Zi5%#-*(0"/08) J (sin’@ sin 28+28 cos’ 6) exp[—2i(w/Qy)B 





— (2°/2%) (B" cos” 6+sin’B sin’ 6)] dB. 


In equation (63) the factor before the exponential differs only by the term 
2iw/Qy from the derivative of the exponent with respect to 8. This makes a 
further simplification of the integral possible. The resulting equation 







soo 
(64) 1+k pt/Z = (2ia/a) J exp[—2i(w/Qz) 8 — (2° /M%) (B" cos” 6 





+sin’B sin® @)] dB 


and equation (49) or (51), rewritten with the subscript j, yield all the quanti- 
ties required for the calculation of the spectrum function o, from equations 
(16) and (25). 

In the special case when k is parallel to the magnetic field, equation (64) 
for £ may be shown to be equivalent to equation (35) for &; in the absence of 
the magnetic field, by execution of the indicated integrations. 

Equation (64) remains nearly equivalent to equation (35) as long as 
Q cos 6/2 >> 1, because then sin 8 can be replaced by 8 in the integrand 
for small values of 8. For large values of 8 the integrand is vanishingly small. 

A sometimes more convenient form of equation (64) is obtained by the 
substitution 28 = y and subsequent expansion of the factor exp(#2?Q,* 
sin? 8 cos y) with the aid of equation (50). This results in 















ctoo 
(65) 1+kpEt/Z = (iw/Qz) exp(—0"/208) f exp(— wy y — 4y72°Qq’ cos” 8) 
0 





: I,(2? sin” 6/22%) cos(ny)dy. 


n=—co 


The integrals occurring in equation (65) may be evaluated, using the identity 







co a/b 
(66) f exp(—iax—1b°x’) dx = b* exp(—a"/8*)( wt—2i f exp y'dy) ‘ 
0 0 





Ci 
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The resulting expansion is given by 


2,2 Je 2” sin’ 6 
(67) I+kint/Z = Q cos 0° exp( 2) ee 4 (ear? ) 


(o= = ila) alg)" } | ‘eae ; | 
exp| - 7 cos? @ IL*” +2 ‘. exp y dy |. 


If this expansion converges rapidly, it may be more useful for numerical 
calculations than equation (64) because tables for the functions f ; exp y*dy 
and J,,(x) exist. 


8. FREQUENCY SPECTRUM OF SCATTERED POWER FOR A 
COLLISION-FREE PLASMA IN A UNIFORM MAGNETIC FIELD 

The frequency spectrum is formally given by equations (16), (25), (49) 
or (51), and (64) or (67) in the presence of a uniform magnetic field. It may 
be necessary in some cases to use these equations for the evaluation of the 
spectrum but more often certain simple approximations are sufficient. 

Before a description of any particular approximation is presented, a justi- 
fication of the terminology introduced in Section 5 for the terms of equation 

25) is given. It is shown that the spectrum function o, can often be separated 

into two terms, the electronic and ionic term, since the mass of an ion is very 
much greater than the mass of an electron. For simplicity consider an ionized 
gas with only a single type of singly charged positive ion. In that case £ 
may be sometimes replaced by its limiting value for high frequencies in the 
first term of equation (25) and similarly £, may be replaced by its limiting 
value for low frequencies in the second term. The justification for these 
approximations is the fact that |s,,|? usually represents a very much wider 
spectrum than |s,2|? and similarly the characteristic variation of £ occurs 
over a much wider range of frequencies than that of £5. 

It is clear from equation (67) that for finite values of cos @ the right-hand 
side vanishes for w — 0 and thus lim, ,o 1 = R~*/5%. It is similarly clear from 
equation (63) that lim,,..f2 = 0. The simplified form of equation (25) be- 
comes 


|seal? |seo[? 
(68) Is]? = (1+é,|? IT TTR Ip —Rlpta|” 


where the first term depends only on the electronic parameters and the second 
term only on the ionic parameters. Equation (68) shows therefore that within 
the limits of its validity the effect of the magnetic field on the electronic and 
ionic spectra may be considered separately. 

The approximations of equation (68) break down when Qg; > 1, Qy2 K Qe 
and at the same time cos @ K Q2/Q,; because then |s,:|? represents a narrower 
spectrum (width of order Q, cos @) than |s;2|? (width of order Q2). Details of 
this case will be discussed later. 

Since both &, and £: are proportional to k~*/;? the second term vanishes 
for klp > 1 and the first term vanishes for klp K1 in equation (68). With 
these results in mind three particular approximations valid in different limiting 
cases will now be considered. 
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(a) Strong Magnetic Field (Qy > Q) 

If the magnetic field is sufficiently strong, so that Qg >> for both ions 
and electrons, then equation (64) becomes equivalent to equation (35) on 
replacement of Q,; in the latter equation by Qcos@ since the factor 
exp (— 2°07? sin? 6 sin? @) can then be replaced by unity. Similarly equation 
(49) becomes equivalent to equation (19) on replacement of 2, = k(2KT/m)! 
by Q cos @ in the latter. 

Thus, if Qy > (i.e., if the radius of gyration is much smaller than k-!) 
then the effect of the magnetic field is a reduction in the width of the spectrum 
by a factor cos 6; the effect in this case only depends on the direction of the 
field. In the extreme case of cos @ = 0 (i.e., when k is perpendicular to the 
magnetic field) the spectrum consists of only a single line at the frequency 
of the incident wave. The total scattered power is not changed by the magnetic 
field. 

In ionospheric applications the condition Qy > Q is sometimes satisfied by 
electrons but not by ions. In such a case the above approximation applies 
only to the electronic spectrum. If at the same time &/p < 1, then the power 
in the electronic spectrum is negligibly small and the approximation is of no 
interest. 


(b) Weak Magnetic Field (1.e., 9g K 2) and k Sufficiently Far from Perpendicu- 
larity to the Magnetic Field (i.e., cos 0 > Qy/Q) 

It has already been shown in Sections 6 and 7 that, for cos @ > Qg/Q2 equa- 
tions (49) and (64) for |s,|? and é in the presence of the field are equivalent 
to equation (19) and (35). The magnetic field thus has no effect on the 
spectrum under these conditions. 

If the condition cos @ > Q,/Q is fulfilled only for ions, then the present 
approximation applies only to the ionic part of the spectrum. If at the same 
time k/p > 1, then the power in the ionic spectrum is vanishingly small and 
the approximation is of no interest. 


(c) Weak Magnetic Field (i.e., M94 KQ) and k Nearly Perpendicular to the 
Magnetic Field (i.e., cos 0K Qy/Q) 

Laaspere (1960) has already discussed the case for k/lp > 1 and has shown 
that a line spectrum is obtained if @ = 2/2. The envelope of the line spectrum 
has the same shape as the spectrum in the absence of the magnetic field. The 
frequencies of the spectral lines differ by exact multiples of the electron 
gyro-frequency from the frequency of the incident wave. The spectral lines 
spread out if the value of cos @ is increased from zero, and for values of @ not 
very far from 7/2 the lines are smeared out altogether so that approximation 
(6) applies. 

The purpose of the present subsection is to show that, if equation (68) is 
valid, a similar approximation exists for other values of k/p. The spectral 
lines do not differ from the frequency of the incident wave by exact multiples 
of the electronic or ionic gyro-frequency, however, and there are two sets of 
lines instead of one. The separation between the lines is close to the electronic 
gyro-frequency in one of these sets, and to the ionic gyro-frequency in the 
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other set. The envelopes of the two sets of lines have the respective shapes of 
the electronic and ionic spectrum that are found in the absence of a magnetic 
field. 

The approximation is derived here only for kip K 1, which is the solution 
of greatest interest in ionospheric applications. A more general approximation 
may be derived by the same techniques for arbitrary values of klp, for both 
the electronic and the ionic parts of the spectrum, but the additional algebra 
is of no particular interest. 

For klp K 1 equation (68) may be — by 


(69) 


For sufficiently small values of cos @ equation (67) is well approximated by 


22 


(70) 1+k pf = exp(—23/ 2072) . Fal QF 208){ ‘sige Mle 


2 292, 
n=—co o—n Que 


+in* (w/M cos 6) exp| - (o- =nQy2) |} 


03 cos’ 6 


If @ = x/2, then the imaginary part of the right-hand side of equation (70) 
vanishes for w # nQy2. The expression on the right of (69) therefore has 
infinities for all those values w,, of the frequency w for which M = Re(1—k? 
i?&2) = 0, where from equation (70) 


9 


(71) M = 2—exp(—93/20}32) D> I, (03/2052) ———-- 
n=—co wo—n ee. 

If 0 is very nearly but not quite equal to 7/2, then in the vicinity of one 
of the roots w, of MZ = 0 one may write approximately 
(72) 1—k' bjt: = Pm(w—wm) +iQm 
where Pm = (0M/0) um and Qn = Im(1—Rl2£2)s-0,- A comparison of 
(51) and (70) shows that 
(73) (|S: 2|”)onwm = —2- © VtoNw 'Qm 
and therefore 
(74) OG» = —t!Now'Qn/|Pm(w— wm) +iQm|? 
from equations (16) and (69). 

If o, is integrated over the frequencies in the neighborhood of the spectral 
line at w = wm, using equation (74), the result is 


wm + Aw 
(75) J ow = Now '|Pnl. 


m— Aw 


The roots w, of the equation M = 0, with M given by (71), have to be 
determined before the application of equation (75), which gives the scattered 
power contained in the spectral lines. 
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The summation indicated in equation (71) may be replaced by integration 
for those terms for which 


(76) |nQy2— | > Oye 


since Qy2 K Qe. In such an integration it is not necessary to exclude those 
terms for which condition (76) is not fulfilled, since they in any case do not 
contribute to the integral. 

The result M’ of the integration does not account for all the terms on the 
right of (71). The summation of terms in the vicinity of w may, however, be 
carried out separately and yields M”, so that M = M’+M”. It will be seen 
after carrying out the integration and summation that the result M’ of the 
integration is a relatively slowly varying function of w, while the result M” 
of the summation is a far more rapidly varying function if Qy2 «K Q:. For 
this reason the full expression M = M’+M” will have to be used only in 
the determination of the roots w,, of M = 0, while in the calculation of the 
derivative Pm = (0M/0w)o-v0, it is sufficient to set M = M”. 

The integration indicated above is carried out first. From equation (71) 


co 2 
(77) M' = 2—exp(—03/20%:) f I,(23/2Q}72) — 


on Oye 


For large values of the argument one may approximate the modified Bessel 
function (Watson 1948) by 


(78) I,(z) & (242) exp(z—n?/2z2). 
Substitution of (78) into (77) and integration leads to 


wo / Qe 
(79) M’ = 2—2(w/Q2) exp(—w'/at) f exp y'dy 


which is the required contribution from the integration procedure. It is con- 
venient to introduce the notation wg = w— (m+4)Qg2 for the summation over 
terms for which nQg2 ~ w where w is taken to lie between mQg2 and (m+1)Qyo. 
Then » may be replaced by m in all the factors under the summation sign in 
equation (71) with the exception of the factor (ow—nQg2)-! = [wa+(m+n—3) 
Qy2|~'. Equation (71) then gives, using mQy2 ~ w, 


it mm 2 2 2 2 = ee eee 
(80)  M” = —exp(—92/203s) In(93/20%2) >- Tita,’ 


n=—oo Wa 


Substitution of J, from (78) and combination of the terms 2 = m+p—1 
and = m—p> leads to 


(81) M” = —27*(Qq2/22) exp(—w?/22) wag z [wi— (+4) Qe}. 


The summation indicated in (81) can be carried out by contour integration 
which yields the identity (Titchmarsh 1939) 


tanz = 22 > [(n+3)?r’—2°J"'. 
n=0 
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Equation (81) then takes the form 


(82) M” = 23(w/Q2) exp(—w"/O3) tan(mwy/Qy2). 


The root wan = &m— (m+4)Qxg2 of the equation M’+M” = 0 is then obtained 
from equations (82) and (79) as 


w/Qe 
(83) tan (mwgm/Qy2) = 2x] — (Q2/w) exp(w’/23) +f : exp yay | 


where the argument of the tangent is taken to lie between —2/2 and 2/2. 
This equation gives the position of the spectral lines for @ = 2/2. (Similar 
equations for the position of both electronic and ionic spectral lines are easily 


derived for general values of lp.) 
The quantity Pn = (0M”/dw¢)eg-wem May now be evaluated from equa- 
tions (82) and (83) with the result 


(84) Pm = 1'Qq}(w/22) exp(—0"/08){ «*+4) — (02/e) exp (w’ /23) 


wo/Qe 2 
+f exp yay | \ 


Substitution of (84) into (75) results in 


@m+4eo 
(85) J ow = x *Noae(Qu2/ 22) 


Wm—Ao 


exp (—w’?/23) 





which proves the original contention that the envelope of the spectral lines 
has the same shape as the spectrum given by equation (51) of Fejer (1960), 
which is found in the absence of the magnetic field. Compariso.: with equation 
(37) of the present paper shows in addition that the total scattered power is 
unchanged, as expected. 

The amount of spread in the spectral lines, as cos @ is increased from zero, 
can also be determined. Equation (74) shows that at half power the spread 


Anw is given by 
(86) Aw = 2\0,./Pal. 


For sufficiently small values of cos @ equations (70) and (78) give 





- Q 2 22 
(87) |Qm| = eo exp[— (w"/23) be (Wam+ 4Qq2)"/Q3c0s 6] 


and P,» is given by (84). 

Although equations (84), (86), (87) give the spread in the spectral lines 
quite accurately for sufficiently small values of cos 6, they are insufficiently 
accurate in the region cos @ ~ Qy2/Q, where the transition from a spiky to 


a smeared-out spectrum occurs. 
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The displacement of the spectral lines, for kip K1, from the positions 
which they would occupy for k/p > 1, is illustrated by Fig. 4 which shows 
(Wam +Qy2/2)/Qu2 as a function of w/Q». 


1.2 16 20° -@A° 2A 
w/Qe0 


Fic. 4. The ratio of the displacement of spectral lines to the ion gyro-frequency as a 
function of w/Q: for klip K 1. 


Three conditions were given for the validity of the approximation for the 
ionic spectrum discussed in the present subsection. These three conditions 
Rly K1, Qys K Qe, and cos 6 K Qy2/Q. are sufficient if at the same time 
Qy1 K Q;. If, however, this latter condition is not satisfied then the additional 
condition cos 6 >> 22/2, = (m,/m:)! also has to be fulfilled because, as has 
been mentioned previously, equation (68) breaks down when Q yg; > Q,, 
Qy2 K Qe, and cos @ K Q2/Q,;. There is usually an intermediate range of angles 
which satisfy the inequalities Qg2/Q:, > cos @ > 22/Q,, and for which the 
approximation of the present subsection still applies even if Qg: > Q,. 

The case Qy; > 21, cos 8 K Q2/Q; is not very important in practice because 
it occurs only over a very narrow range of angles. For atomic oxygen ions 6 
would have to differ by less than one quarter of a degree from 90°. By con- 
siderations similar to those of the present subsection it may be shown that 
in this case practically no power is contained in those spikes for which 
w >> Q, cos @ and thus the spectrum, while it still consists of spikes separated 
approximately by the ion gyro-frequency, is narrower than the ionic spectrum 
in the absence of the magnetic field. For 6 = 7/2 the spectrum consists of a 
single line at the frequency of the incident wave. 


CONCLUSIONS 


It has been shown, taking the magnetic field into account, that the spectrum 
function o can usually be divided into two parts, one of which is determined 
by the electrons, the other by the ions. The former dominates for k/p > 1, 
the latter for klp <1. This result makes a separate discussion of the elec- 
tronic and ionic components of the spectrum possible. It has also been shown 
that the total power scattered is not affected by the magnetic field. 
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If Qy > Q, for either electrons or ions then the magnetic field causes a 
contraction of the corresponding spectrum by a factor cos 6, where @ is the 
angle between k and the magnetic field. 

If Qy KQ for either of the particles, and at the same time cos @ > Q,/Q 
then the magnetic field has practically no effect on the corresponding spectrum. 

If Q9g K 2 for both ions and electrons but cos @ K Qy/Q, then each spectrum 
has sharp maxima at frequencies separated approximately by the gyro- 
frequency of the corresponding particle. For 6 = 2/2 these maxima become 
discreet spectral lines within the approximations of the present theory. The 
frequencies of these spectral lines differ by exact whole multiples of the 
gyro-frequency of electrons from the frequency of the incident wave, for 
klp > 1. For kip K1 the spectral lines are displaced by amounts smaller 
than one-fifth of the ion gyro-frequency from positions which they would 
occupy if they differed by exact whole multiples of the ion gyro-frequency 
from the frequency of the incident wave. If one disregards the fine structure 
of the spectrum as manifested in these closely spaced maxima or discreet 
lines, then it is shown that the distribution of power in the spectrum, for 
Qy K Q, is unchanged by the presence-of the magnetic field. 

If klp K1, Qy2 K Q2 but Qy: > Q, then there is an intermediate range 
of angles, given by the condition (Qgi/Q1) (Q2/21) > cos @ > 22/2), for which 
a spiky spectrum is again obtained in which, apart from the spiky structure, 
the distribution of power is still the same as in the absence of the magnetic 
field. For angles 6 still nearer to r/2, however, the distribution of power in 
the spikes becomes progressively narrower. For 6 = 2/2 all the scattered 
power is contained in a single line at the frequency of the incident wave. 

In ionospheric experiments, at all frequencies likely to be used, the con- 
dition Qg K Q is satisfied for oxygen ions but not necessarily for hydrogen 
ions or electrons. 

In experiments of Bowles (1958, 1959) and Pineo, Kraft, and Briscoe (1960), 
atomic oxygen ions were dominant and most of the energy was contained in 
the ionic part of the spectrum since Alp <1. Further, the inequality 
Qy2/Q2 cos 0K 1 was certainly satisfied by oxygen ions, and therefore the 
effect of the magnetic field on the spectrum was small. 

In other experiments, however, the magnetic field could have an important 
effect on the spectrum as the following examples show. 

For hydrogen ions, Qg2 and Q2 may be of the same order of magnitude. If 
at the same time &/p < 1 and k is nearly but not quite perpendicular to the 
magnetic field, then a large part of the scattered energy is contained in the 
two slightly broadened spectral lines whose frequency differs from the frequency 
of the incident wave approximately by the ion gyro-frequency. The existence 
of these spectral lines has been predicted by Bowles (1959). 

Similarly, if the electron gyro-frequency Qg: is of the same order of mag- 
nitude as Q;, and at the same time &lp > 1 and k is perpendicular to the 
magnetic field, then a large part of the scattered power is contained in those 
spectral lines whose frequency differs from the frequency of the incident wave 
by the electron gyro-frequency. 
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In both the above examples, the amount of scattered power contained in 
the spectral lines could be computed from the results of the present paper. 
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EFFECT OF TWO-DIMENSIONAL MECHANICAL STRESS ON 
THE DIELECTRIC PROPERTIES OF POLED 
CERAMIC BARIUM TITANATE AND LEAD ZIRCONATE TITANATE! 


R. F. BRown 


ABSTRACT 


Studies have been made of the dielectric behavior of several ferroelectric 
ceramic materials when a two-dimensional stress was applied normal to the 
axis of polarization. It has been shown that the dielectric constant decreases 
and the dielectric loss increases with increasing stress, part of the change being 
irreversible and part reversible. Upon application of stress to a sample, the 
dielectric constant did not change instantaneously but appeared to decrease 
linearly with the logarithm of time. Mechanisms are suggested for the observed 
effects. 


INTRODUCTION 


It is well known that the dielectric constant of barium titanate is dependent 
on the direction of measurement relative to the crystal axis. The dielectric 
constant of a pure barium titanate single crystal is about 300 in the direction 
of the spontaneous polarization (the c-axis) but about 4500 normal to the 
direction of polarization (the a-axis). Hence, in multidomain crystals, the 
dielectric constant is dependent upon the domain configuration and the 
orientation of the domains relative to the direction of measurement. If the 
direction of measurement is kept constant and the domain configuration is 
changed, the dielectric constant may be modified. 

One method of altering the domain orientation is by the application of 
mechanical stress. Upon compression, the domains with the c-axis normal to 
the direction of stress have the lowest energy. Thus, the application of com- 
pressional stress to a multidomain crystal of barium titanate favors domains 
with a component of the c-axis normal to the applied stress. The dielectric 
constant normal to the direction of stress would be decreased while that in 
the direction of the stress would be increased. 

Polycrystalline ferroelectric ceramics have a more complex structure than 
single crystals and the c-axes of the individual crystallites are randomly 
oriented, resulting in zero net polarization. However, a strong d-c. electric 
field tends to align the axes of spontaneous polarization of the domains in 
the direction of the external field, and the alignment remaining after removal 
of the field results in a remanent polarization in the ceramic. The application 
of stress to poled barium titanate ceramic would be expected to affect the 
dielectric constant, the magnitude and sign of the change being dependent 
upon the stresses and directions of application relative to the direction of 
polarization. 

The effect of unidirectional stress on the dielectric properties of ceramic 

1Manuscript received November 11, 1960. 
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barium titanate has been reported by a number of investigators in recent years 
(Takagi et al. 1948; Shirane and Sato 1951; Sinyakov and Izhak 1955; and 
Izhak 1956, 1957). In addition to unidirectional stress, the effects of two- 
dimensional stress have been studied by Roi (1955) and Ksendzov and Roten- 
berg (1959). The results presented in these papers show disagreement in the 
observed direction of change of the dielectric constant with stress, but this 
may be due in part to differences in the conditions of measurement such as the 
state of polarization of the samples, their history, and the time of measure- 
ment after the application of stress. 

It has been observed by several investigators that the change in the dielectric 
constant of ceramic barium titanate which occurs after a change in stress is 
time dependent. Izhak (1957) showed that the application of a unidirectional 
compressional stress produced a rapid increase in dielectric constant followed 
by a gradual decrease to a steady value after several minutes. 

This paper describes the effects of two-dimensional compressional stress 
on the low-field dielectric constant and dielectric loss of poled samples of two 
modified barium titanate compositions and two modified lead zirconate 
titanate compositions. Lead zirconate titanate has a crystal symmetry and 
domain structure similar to that of barium titanate and it would be expected 
that the effect of stress on the dielectric properties would be similar. It is 
found for both types of ceramic that, except for small values of stress, the 
dielectric constant in the direction normal to the stress decreases as the 
compression is increased. Most of this change is reversible, but there is usually 
a permanent drop in dielectric constant after the first pressure cycle on a 
fresh sample. Subsequent pressure cycles produce nearly reversible changes. 
The dielectric loss is found to increase as the stress is increased and falls to a 
value slightly higher than the initial value when the stress is released. 

Quantitative measurements have been made of the aging process which is 
initiated by the stress, and it is shown that except for short times of less than 
about one minute, the dielectric constant usually decays substantially linearly 
with the logarithm of time, similar to the aging which follows poling of a 
ferroelectric ceramic. The aging rate is found to depend on the ceramic com- 
position and the history of the sample. 








































STRESS SYSTEMS 


A thin-wall sphere subjected to external hydrostatic pressure offers a 
simple means of obtaining high two-dimensional compressional stress in the 
wall for moderate external pressures. If a sphere of outside radius ‘‘b’’ and 
inside radius ‘‘a’’ is subjected to uniform external pressure Po, the stress at 
any point at radius ‘‘r’’ in the wall of the sphere is fully described by three 
orthogonal stress components which are defined as the principal stresses. 
These are given by Love (1927): 

(1) the tangential stresses 












09 = 0% = 
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b\*| r?—a’® 
a, = Po , Pa 2 


If the wall is thin relative to the radius, the radial component of stress can 
be neglected compared to the tangential components, which become 


(2) the radial stress 


9 = 04 = 8 ahs ‘J 
2 b°—a 
and the stress is essentially two-dimensional. 

A thin-wall circular cylinder with ends closed by rigid caps and subjected 
to external hydrostatic pressure is also essentially a two-dimensional stress 
system. The two principal stress components are given approximately by Love 
(1927): 

(1) the tangential stress 


orn ole], 


(2) the axial stress 


b° 
a, = Po a | = 9/2. 


Both the spherical and cylindrical samples of ceramic used in these experi- 
ments had silver electrodes on inner and outer surfaces and were poled radially. 
The capacity measured between electrodes was thus in a direction normal to 
the two principal stress components, and along the induced ferroelectric axis. 

Although the stress tensors are different for spheres and end-capped cylinders 
exposed to external hydrostatic pressure, in both cases the important principal 
stresses are normal to the axis of polarization of the ceramic. If a; and a» 
represent the principal stresses in both two-dimensional systems, it has been 
found experimentally that the variation of dielectric constant as a function 
of the sum of the principal stresses, o:+¢2, is approximately the same for 
spheres and cylinders of the same material. Most of the measurements reported 
here were made on cylindrical samples since these were more readily available 
than spherical samples. 


DESCRIPTION OF CERAMIC SAMPLES 


The hollow spheres were assembled by cementing together matched hemi- 
spherical shells of barium titanate or lead zirconate titanate ceramic with a 
nominal outside diameter of 3 in. and wall thickness of } in. The shield of a 
coaxial cable was connected to the outer electrode of the sphere and the core 
was connected to the inner electrode via a small high pressure glass-to-metal 
seal cemented in the wall of the sphere. A thin waterproof coating was applied 
to the assembled sphere by dipping it in liquid neoprene. 

The cylindrical test specimens (Fig. 1) consisted of ceramic cylinders of 
1} in. nominal outside diameter, 1} in. length, and 3 in. wall thickness, closed 
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Fic. 1. Cross section of cylindrical ceramic specimen. 


on the ends with rigid plastic caps. A length of high pressure steel tubing was 
moulded into one end cap to allow mounting of the specimen inside the high 
pressure tank with the electrical lead to the inner electrode running through the 
tube and that to the outer electrode connected to the tubing. The cylindrical 
samples were also given a thin waterproof covering of neoprene. 

Four different ceramic compositions were investigated. The first was a 
common commercial material consisting of barium titanate with 4% by weight 
of lead titanate and 2% calcium titanate, which will be referred to as (Ba, 
Ca,Pb)TiO;. The second was (Ba,Ca,Co)TiO;, a composition described by 
Schofield and Brown (1957) consisting of barium titanate with additions of 
5% by weight calcium titanate and 2% cobalt carbonate. The lead zirconate 
titanate compositions were commercial PZT-4* and PZT-5*. All the samples 
were poled several months prior to these measurements, so that they were 
in a quasi-stable state and the changes due to stress application could be 
considered isolated from the changes due to normal aging of the ceramic. 


EXPERIMENTAL METHOD 


The pressure system consisted of a water-filled tank capable of 1000 kg/cm? 
hydrostatic pressure, monitored by a Bourdon gauge which had been calibrated 
against a dead-weight tester. Measurements of the capacity and dielectric loss 
of the samples were made with a General Radio Type 716C capacity bridge 
using a signal of about 1 volt across the samples at | kc/s. From the measured 
capacity and dimensions of the samples, the permittivity or dielectric constant 
was calculated. The stress in the wall of each sample was calculated from the 
measured dimensions of the sample and the known external pressure. Although 
no attempt was made at temperature control of the pressure vessel, the 
temperature remained in the region 21° to 23° C during the measurements, a 
range in which the dielectric constant of all four compositions is essentially 


temperature independent. 
*PZT-4 and PZT-5 are registered trademarks of the Clevite Corporation. 
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Measurements of capacity and dielectric loss were made as a function of 
stress as the external pressure was increased in steps of 70 kg/cm? to a maxi- 
mum of 490 kg/cm? and then reduced to zero in similar steps. The maximum 
pressure was limited to keep the stresses well within the estimated com- 
pressional strength of the materials. Because of the time dependence of 
permittivity after a change in stress, a standard time reference was established 
by taking the measurements 10 minutes after each change of pressure. The 
pressure was cycled at least four times in order to determine the permanent 
and reversible components of change in dielectric constant and the influence 
of each cycle on the following one. At the maximum pressure, the variation of 
capacity with time was recorded for a period of at least 2 hours on the first 
three pressure cycles, and for 2 or 3 days on the fourth cycle. 


EXPERIMENTAL RESULTS 
(a) Effects of Stress 

The dielectric constant of cylindrical samples of each of the four ceramic 
materials is plotted against the sum of the principal stresses, o7 = og+c, in 
Figs. 2 to 5 for the first two pressure cycles. With one exception, the dielectric 
constant on the third and fourth pressure cycles was only very slightly less 
than that on the second cycle, and is not shown in the figures. However, the 
fourth cycle is shown for PZT-4 because there was an appreciable change in 
dielectric constant. The drop in dielectric constant at the maximum pressure, 
which can be seen in the figures, occurred during the aging investigations which 
will be shown in detail later. After the release of pressure at the end of each 
cycle, there was an 18-hour interval before the start of the next pressure cycle 
during which time the dielectric constant of PZT-5 increased slightly while 
that of the other three compositions decreased. This change amounted to 
about 6% for PZT-4 but was less than 1% for the other materials. 

On the fourth cycle, the maximum pressure was held for 4 days while 
measurements of dielectric constant and loss were made on all samples at 
convenient time intervals. It will be seen from Fig. 4 that a relatively large 
change occurred in the dielectric constant of PZT-4 which was not recovered 
after the release of the pressure. The changes which occurred in the other 
materials during this time were much smaller, and during the decrease of 
pressure to zero the dielectric constant followed closely the curves of the 
first two cycles. 

In order to compare the different compositions graphically, the dielectric 
constant measured on the increasing leg of the second pressure cycle has been 
normalized and plotted in Fig. 6 against the sum of the principal stresses in 
the samples. 

The dielectric loss, which is also shown in Figs. 2 to 5, increased with 
increasing stress and returned to a value slightly higher than the initial value 
after the release of stress. Under constant stress, the dielectric loss decreased 
slightly with time. The dielectric loss on successive pressure cycles is not shown 
since there was little change from the first cycle. 

Spherical samples were available in (Ba,Ca,Pb)TiO; and PZT-4 ceramic. 
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Fics. 2-5. Effect of two-dimensional stress on the dielectric constant and dielectric loss 
of ceramic BaTiO; and PZT. 
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Fic. 6. Comparison of the change of dielectric constant with stress for different materials. 


For these, the normalized dielectric constant measured on the second in- 
creasing pressure cycle is plotted in Figs. 7 and 8 against the sum of the 
principal stresses, or = og+o4. On the same graphs, the dielectric constant 
for cylindrical samples in the same cyclic condition is plotted against the 
sum of the principal stresses, 77 = o9+¢,. Even though the ratios of the 
principal stresses are o9/o, = 1 for the sphere and o¢/o, = 2 for the cylinder, 
the curves coincide within 2%. 


(b) Aging Phenomena 

In these experiments it was found that the capacity did not reach a constant 
value within a few minutes after a change in stress, but continued to decrease 
over the period of measurement, which sometimes extended to several days. 
The dielectric constant of the four ceramic materials is plotted against the 
logarithm of time in Figs. 9 to 12. The data for these curves were obtained 
after the external pressure had been raised from 420 to 490 kg/cm? on each of 
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Fics. 7-8. Variation of the dielectric constant with stress for different ratios of the principal 
stresses. 


the four successive pressure cycles which have been described. It should be 
noted that the method of measurement did not allow observation of the 
dielectric behavior immediately after a change in stress, and consequently the 
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Fics. 9-10. 


initial increase in dielectric constant noted by Izhak (1957) could not be 
observed. 

For comparison of the four materials, the dielectric constant measured on 
the fourth time run has been normalized, taking the 1-minute value as K/Ko = 
1.00, and the curves plotted together in Fig. 13. It is evident that the dielectric 
constant for three of the materials continues to fit a linear relationship with 
log ¢ quite closely, but the slope of the curve for PZT-4 increases with time. 
The rates of decrease, expressed as the percentage change in the decade from 
10 minutes to 100 minutes after application of pressure, are tabulated in 
Table I for the various cycles for each composition. It should be noted that 
the values listed for PZT-4 are average values for this decade, since the 
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_Fics, 9-12. Aging of the dielectric constant after four successive applications of two- 
dimensional stress. 


TABLE I 
Rate of stress-induced aging 








Rate of change of dielectric constant (per cent per 
decade) for a change in external pressure from 420 
to 490 kg/cm? 


Ceramic material Ist cycle 2ndcycle 3rdcycle 4th cycle 


(Ba,Ca,Pb)TiO; 2.5 ; 2.0 
(Ba,Ca,Co)TiO; 1.9 4 ; 
PZT-4 4.4 ; 2 
PZT-5 1.9 
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Fic. 13. Comparison of stress-induced aging for different materials. 


decrease is not linear with log ¢. The fact that the rate of aging is higher on 
the first time run than on successive runs indicates a certain amount of stabili- 
zation in the ceramic after the first application of high stress. 


DISCUSSION 
(a) Effect of Stress 


There are two components contributing to the observed dielectric constant 
or permittivity in a ferroelectric ceramic. One is the intrinsic permittivity of 
domains, and the other is a contribution to permittivity due to the oscillation 
of 90° domain walls as described, for example, by Plessner (1956). 

The effect of stress on the intrinsic permittivity of a single domain may be 
predicted from the thermodynamic theory of ferroelectricity. This shows that 
the changes in permittivity which occur when stress is applied normal to the 
c-axis or normal to the a-axis are opposite in sign. In a ceramic, only a portion 
of domains are aligned with c-axes along the poling direction, so that domains 
may be compressed in either the a or ¢ directions by the external pressure. 
The resultant changes in dielectric constant would have opposite signs and 
tend to cancel, so that only a small contribution to the total dielectric constant 
would be expected from changes in the intrinsic permittivity of individual 
domains. 

Probably a more important effect is the fact that compressional stress tends 
to align domains so that their polar axes (c-axes) are normal to the direction 
of applied stress. Since the intrinsic permittivity in the ¢ direction is much 
less than that in the a direction, this domain alignment should result in a 
reduction of the contribution to total dielectric constant from the intrinsic 
permittivity of domains. 
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However, it is also necessary to take account of the contribution from 
domain wall oscillation. In order to explain the experimentally observed 
increase in total dielectric constant which occurs for small compressional 
stresses in some of the ceramic compositions, it is postulated that the appli- 
cation of this stress facilitates the oscillation of domain walls in response to 
a small alternating voltage, that is, that it increases domain wall mobility. 
This would result in an increase in the contribution to permittivity due to 
domain wall oscillation. 

The observed behavior can be explained by assuming that, as the stress is 
increased, greater wall mobility causes the dielectric constant to increase 
initially, but that domain alignment soon becomes the controlling factor 
resulting in a decrease in the dielectric constant. The relative magnitude of 
the two effects depends on the obstacles to domain wall motion and domain 
alignment and can therefore be expected to be composition dependent. The 
process of domain alignment would be expected to continue with increasing 
stress until all domains were aligned with c-axes normal to the stress. This 
would result in the dielectric constant decreasing to a limiting value, and there 
is some experimental evidence in Figs. 2 to 5 that the curves do tend to level 
out for sufficiently high stress. 

The failure of the dielectric constant to return to the initial value after a 
pressure cycle indicates some permanent reorientation of domains. However, 
the greater part of the change is reversible, showing that domain reorientation 
is largely a reversible process. This is in agreement with strain measurements 
made by Subbarao, McQuarrie, and Buessem (1957), which showed little 
permanent domain reorientation after the release of a compressive stress. 

The increase in dielectric loss with stress is explained by the assumption of 
higher domain wall mobility, since domain wall oscillation is inherently an 
energy-dissipating process. The fact that this loss gradually decreases with 
time when the stress is held constant indicates that the wall mobility becomes 
less with time as the walls assume positions of greater stability. 


(b) Aging Phenomena 

The aging of permittivity which occurs after a change in stress is similar to 
that initiated when ceramic samples are freshly poled, and thus seems likely 
to be due to a similar process. The aging which occurs after poling has been 
ascribed to the motion of domain boundaries to more stable positions. By 
analogy with ferromagnetism, McQuarrie and Buessem (1955) and Plessner 
(1956) have shown that, except for very short and very long times after poling, 
the aging process should follow a log ¢ relationship if it is assumed that the 
obstacles to domain wall motion have a random character, resulting in a wide 
distribution of activation energies. However, the aging of some ceramic samples 
has been found to deviate from a linear relation with log t, sometimes increasing 
and sometimes decreasing in slope with time. Therefore it is not surprising 
that the stress-induced aging in PZT-4 is found to deviate from the simple 
log ¢ relation, although the mechanisms are not well understood. A more 
detailed study of the pressure dependence of the aging rate is continuing and 
will be reported separately. 
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INVESTIGATIONS ON THE DIFFUSION OF MINORITY 
CARRIERS FROM A POINT ON SILICON} 


C. H. CHAMPNEss? 


ABSTRACT 


The field-free diffusion of minority carriers injected at a point in silicon has 
been studied and the time from injection to the maximum of the collector signal 
due to the arriving carriers measured for various emitter-collector distances. 
Firstly, it was found that for 2-ohm-cm, 6-microsecond, n-type material the 
time to maximum was proportional approximately to the emitter-collector 
spacing raised to the power 1.2 for flat surfaces and to the power 1.6 for 3- and 
6-degree wedge samples. The power 1.6 was also found for a 17-degree wedge 
of 120-ohm-cm, 100-microsecond, p-type silicon. Secondly, decay times following 
the maxima appeared to increase with emitter-collector spacing. No adequate 
two-surface theory is available but comparison with a one-surface theory revealed 
serious disagreements. It would appear that the results may be explained 
qualitatively by assuming that the effective lifetime is dependent on the excess 
carrier density and decreases as the emitter point is approached. 


INTRODUCTION 


Much of the early work on transistor action and on drift mobility studies 
was done using point contacts on the surface of a semiconductor. There appears, 
however, to have been no special experimental study of the concentration of 
carriers with distance and time as they diffuse out from a point contact with 
no applied electric field. No doubt this is because it was thought that the pro- 


cess of field-free diffusion is a straightforward one which could be accurately 
described by a radial solution of the diffusion equation with a simple recombina- 
tion term. Such a solution assumes only small concentrations. However, 
actually the problem involves a high concentration of extra carriers near the 
point of origin. This could affect both the lifetime and the diffusivity. The 
present experiments were therefore carried out to test the validity of the 
conventional diffusion theory under the circumstance of very heavy localized 
injection at a point contact. 

If minority carriers are injected at a point they will diffuse outwards radially 
and give rise to a collector current through a reversed biased collector point. 
This collector current is approximately proportional to the concentration of 
excess minority carriers near the point (Shockley, Pearson, and Haynes 1949). 
The collector current will increase with time as more carriers reach the point. 
However, for pulsed injection, as was used in the present experiments, it will 
reach a maximum and begin to decrease because of the decrease of carriers 
at the emitter. This is shown in Fig. 1B, which shows a typical oscilloscope 
trace obtained in the experiments. The asymmetry of the arrival pulse is to 
be compared with that in Fig. 1A, which shows a trace recorded in a drift 
mobility experiment where the carriers are transported under an applied 


1Manuscript received November 29, 1960. 

Contribution from the Solid State Physics Laboratory, Research and Development Division, 
Northern Electric Company Limited, Montreal, Que. 

2Now at Physics Department, McGill University, Montreal, Canada. 


Can. J. Phys. Vol. 39 (1961) 
754 





PLATE I 


Fic. 1. Variation of collector signal with time. A shows an oscilloscope trace obtained 
during a drift mobility experiment. B shows the trace with the circuit of Fig. 3, i.e., no applied 
electric field. 
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electric field. In the experiments reported in this paper there was no applied 
electric field. The quantities measured were the time (¢,) from the center of 
the injection pulse to the maximum in the collector signal, the decay of the 
collector signal following the maximum, and the emitter-collector spacing (r). 


THEORY OF DIFFUSION FROM AN INSTANTANEOUS POINT SOURCE 


Before the details of the experimental work are presented, it is helpful to 
review theoretical results for the case of diffusion from an instantaneous point 
source firstly in an infinite and secondly in a semi-infinite medium. 


Infinite Solid 

Suppose N, holes are injected in a very short interval of time at a point 
r = 0 inside an infinite n-type semiconductor. Suppose, too, that there are 
no applied electric fields and that there is no space charge anywhere. With 
Ap, the excess concentration of holes, equal to An, the excess concentration 
of electrons, there is a single lifetime 7 for holes and electrons. The rate of 
change of Ap with time ¢ is then given by 


(1) oot = div (D' grad Ap) —“P 


where D’ is the ambipolar diffusivity (n+ )p)[(m/D,)+(p/D,)]“, Dy is the 
electron diffusivity, D, the hole diffusivity, and m and p respectively are the 
electron and hole concentrations. If div grad is expressed in polar co-ordinates 
with radial dependence only, equation (1) may be solved with the boundary 
conditions N, = f(Ap.4ar2dr as t—0, and Ap—0 as r—> © to give the 
well-known (Shockley 1950) solution 


‘ - —t/2 Bb r \ 
(2) Ap = N,(42D,t) exp) (5+ - 


v 


where D’ has been replaced by D, since 1 > p is assumed. 
This equation predicts qualitatively the approximate shape of the observed 
oscilloscope trace in Fig. 1B. At its maximum, equation (2) yields 


(3) t= ADsty?(r?—6Dytm)—, 
where ¢,, is the time at the maximum. If 7 > 3t,, equation (3) reduces to 
(4) tm = r°/6D, 
and if tr X 3tm we get 
~ r 
(5) bn = = (r/Dy)*. 
Thus in comparison with the time to maximum, if the lifetime is large, ty 
depends quadratically on the emitter—collector spacing but if the lifetime is 


small a linear variation is expected. For lifetimes between the two extremes 
equation (3) gives 


(6) 
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Semi-infinite Solid 

The theory of diffusion in a semi-infinite solid has been worked out in 
great detail by van Roosbroeck (1955), taking surface recombination velocity 
(s) into account. Using his theory the following three cases may be considered : 
(A) emitter and collector on the surface, (B) emitter inside the solid and 
collector on the surface, and (C) emitter on the surface and collector inside 
the solid. At the maximum of Ap the following result holds for the three 
cases 


(7) 7 = 4Dy ta” {r”—6Dy tm(1 +36) }™ 


where ¢ is a number between 0 and 1 which can be evaluated specifically for 
extreme values of s. These are as follows: 

(i) s & (D,/tm)?, i.e. s small, « = 0 for all three cases and hence the results 
obtained for an infinite solid apply, since equation (7) reduces to equation (3). 
(ii) s >> (D,/tm)?, i.e. s large, € = 1 for case (A) but for cases (B) and (C) 


(8) e = [1+(1/2stm)]“. 


(iii) s >> 7/2tm, ie. s still larger (since r/2tm 2 (D,/tm)*), ¢ = | for all three 
cases. Equation (7) then reduces to 


(9) + = ADyty?(r?—10D tm). 


Under condition (ii) for cases (B) and (C), equations (7) and (8) may be 
solved for s to get the following result 
r Xr-1 
(10) $= 3 1-¥r 
where X = (7?—6Dytm)/4D ptm? and Y = (r?—10Djtn)/4D,tn?. Since X > Y, 
we must have, for s to be a positive quantity, that 


(11) A ee eo. 


Using the above notation for X and Y we may rewrite equations (3) and 
(9) as simply 


(3A) T s—0 
(9A) T - 5s © 


EXPERIMENTAL METHOD 


The experiments were mostly carried out using samples cut from a 2-ohm- 
cm doped n-type single crystal of silicon. Subsidiary measurements made on 
the material gave a lifetime of about 6 microseconds, a drift mobility for the 
holes of 230 cm?/volt sec, and a donor concentration of about 10'%/cm*, Experi- 
ments were also made on a sample of 120-ohm-cm p-type Du Pont single 
crystal silicon for which the lifetime was quoted as 100 microseconds. 

The first set of experiments was carried out using a silicon sample in the 
form of a filament (approximately 0.20.21 cm) with electroplated nickel 
end contacts to which wires had been soldered. Two point contacts of tungsten 
wire 0.005 in. in diameter were placed close together on one of the faces of 
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the sample by means of micromanipulators. In the later experiments, the 
ends of the wires were reduced to a finer point by etching. 

In the second set of experiments, which formed the greater part of the 
work, measurements were carried out on wedge-shaped specimens. In these 
experiments, the thick end of the wedge was clamped between spring-loaded 
copper blocks (base connection) and the emitter and collector wires were 
placed opposite to each other on the two inclined faces of the wedge as shown 


in Fig. 2. 


ANGLE EXAGGERATED 









EMITTER COLLECTOR 


MINORITY CARRIERS 


Fic. 2. Cross section of a wedge of silicon showing the position of the emitter and collector 
points. 





EMITTER 





SILICON 


Fic. 3. Diagram of the circuit used to observe the diffusion pulse. 


The circuit used is shown in Fig. 3. The resistor R of 1000 ohms was smaller 
than the collector resistance (> 10‘ ohms), so that the signal across R and 
hence that shown on the oscilloscope was approximately proportional to the 
collector current. 

The forward pulse applied to the emitter was usually less than 1 micro- 
second in length and frequently less than 0.1 of a microsecond. 

The experiments consisted of obtaining different oscilloscope traces for 
different emitter-collector spacings and measuring the variation in the time 
to maximum and in the decay following the maximum. The distances were 
measured with a travelling microscope. 
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EXPERIMENTAL RESULTS 


Time to Maximum Measurements 

Flat Surface Experiments 

These measurements were carried out with a single filamentary sample A3 
and are referred to in the results as the flat sample experiments, since both 
emitter and collector points were on the same surface. 

Figure 4A shows a log-log plot of time to maximum against emitter-collector 
spacing for a surface treatment giving a shiny finish to the sample. This 
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Fic. 4. Log-log plots of the time to collector maximum against the distance between 
emitter and collector points, on the same surface of 2-ohm-cm, n-type silicon sample A3. 
A shows results for a surface etched to give a shiny finish. B shows results for a surface ground 
on carborundum paper and also a surface blasted with pumice powder. 


treatment consisted of first polishing the sample surface with progressively 
finer lapping powders and then etching it for 5 minutes in a 5:1 mixture of 
nitric and hydrofluoric acids. One point was obtained for a mat finish to the 
surface. The treatment in this case was an etch for 2 hours in a mixture of 
acetic, nitric, and hydrofluoric acids in the ratio 10:3:1. Another set of measure- 
ments was made on the same sample after the surface was rubbed on 400A 
carborundum paper. Yet another set of measurements was taken on the same 
sample, this time after air blasting with pumice powder to produce a pitted 
surface. The results for these two surface treatments are given in Fig. 4B. 
For comparison purposes, the three average lines through the points are 
shown redrawn on a single graph in Fig. 6B. The ground surface on the 
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average appears to give the largest time to maximum values and the pumice 
blasted surface the least; the shiny surface values are between the two. From 
the slopes, the following empirical laws are found: 


(12) tm = 4.9 10-4! (shiny surface), 
(13) tn = 4.8X 10-t116 (ground surface), 








(14) tn = 5.6X 10-‘r!-26 (blasted surface). 





Thus it may be said that the experiments give approximately t, « r'? 
regardless of the surface treatment. 

Wedge Experiments 

The wedge experiments were carried out to ensure that diffusion was 
taking place through the bulk and not predominantly along the surface of 
the specimen. The wedge form of geometry enabled different emitter—collector 
distances through the bulk to be used. Three wedges were used. Each was 
prepared by lapping or grinding a specimen of silicon waxed to an inclined 
face set in a block of steel. The apex angles were kept small so that the small 
areas on the faces near the emitter and collector might be regarded, for 
theoretical purposes, as parallel. 

Figure 5A shows the results for a wedge made from the same material as 
sample A3 used in the flat surface experiments. This wedge had an apex 
angle of 3° but later on it was reground without maintaining the original 
angle. Figure 5B shows the results for a sample A2 with an apex angle of 
6°. Figure 6A shows results for a 17° wedge of p-type material. 

The average lines drawn through the experimental points for the three 
wedges are redrawn together on Fig. 6B. The slopes for these lines are greater 
than those for the flat surface experiments. The empirical laws for the three 
lines are as follows: 


(15) E. 
(16) bn = 9.2K 10-7!“ (sample A2), 
(17) je 


Apart from the differences in the proportionality constants, the wedge 
experiments therefore give approximately t, « r'. 





























1.9 10-*7!-% (sample A3), 







7.1X% 10-4" (Du Pont p-type sample). 


Il 









Decay Measurements 

Figure 7 shows typical plots of the logarithm of the collector signal against 
time for three particular emitter—collector spacings. The decay curve in one 
case is a straight line. It is slightly convex in another case and slightly concave 
in the third. Despite the slight curvature observed in some of the cases, it 
was found that log-linear plots gave, in general, the nearest approach to 
straight lines. Other plots (log-log) for instance gave greater curvatures. The 
decays were thus described by a term exp(—t/71) and values of the so-called 
decay time 7; were obtained from the slopes of the lines. These values covered 
a wide range but were not completely random. Figure 8 shows 7; plotted 
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Fic. 5. Log—log plots of time to maximum against distance between emitter and collector 
points on wedges of 2-ohm-cm, -type silicon. A shows results for wedge sample A3 with an 
original apex angle of 3° and B shows results for wedge sample A2 with an angle of 6°. 

Fic. 6. Log-log plots of time to maximum against distance between emitter and collector 
points. A shows results for a 17° wedge sample of 120-ohm-cm, 100-microsecond p-type 
silicon. B shows all the average lines through the experimental points in the previous plots. 
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Fic. 7. Variation of the logarithm of the collector signal against time, taken from three 
typical oscilloscope traces. 


against ¢, and Fig. 9 shows it plotted against 7. While there is considerable 
scatter in the points, indicating that some important factor is not kept con- 
stant, there is a general tendency for 7; to increase with t, or r as indicated 
by the broken line in each case. There was no obvious difference between 
values of 7; for the flat surface and for the wedge experiments, although with 
such scatter present any small difference would not show up. At larger values 
of r (or tm), 71 appears to flatten off somewhat. Linear-log plots for 7; and tm 
and 7; and r appear to have the effect of straightening out the general trends. 
However, it is found that the variation in Fig. 8 may be described in the 
simplest way by a rough empirical formula of the form 


(18) t1 = 10-%m?*®. 


All the experiments were done in semidarkness to get the best sensitivity 
to resolve the maxima. The effect of steady light from a microscope lamp on 
the sample was generally to reduce the amplitude of the maximum without 
noticeably altering the time to maximum. While the decay time was slightly 
reduced under illumination, the variation with 4, (and r) was apparently still 
present. 

Regarding the effect of the collector circuitry, it is estimated from the 
capacitances (preamplifier 25 wuf) and resistances (load resistance 10* ohms) 
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Fic. 8. Decay time r; (determined from log-linear plots of collector signal versus time) 
plotted against time to maximum. 
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Fic. 9. Decay time 7; (determined from log-linear plots of collector signal versus time) 
plotted against emitter-collector spacing. 
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involved that they result in a time constant of less than 0.1 microsecond; the 
capacitance of the collector contact itself was less than 1 puf. This has, there- 
fore, a negligible effect on the decay curves. 


Other Measurements 

Effect of Injection Voltage on Decay Time 

Increasing the amplitude of the injecting pulse did not produce a clear-cut 
effect but the general tendency was for the decay time 7; to increase slightly. 


Effect of Emitter and Collector Fields on tm 

The influence of both fields was examined by measuring times to maximum 
firstly with increasing values of collector bias, secondly with increasing pulse 
generator amplitude settings, and thirdly with increasing injection pulse 
widths, keeping the two other variables constant each time. It was found that 
in all three cases the time to maximum decreased. However, extrapolation to 
zero values showed that a finite value of time to maximum would still be 
obtained for an injection pulse of very small amplitude and width with a 
very small value of collector bias. Thus we are justified in saying that the 
experimentally observed maximum is predominantly due to diffusion and not 
to drift in an electric field. 


ANALYSIS OF RESULTS 


Time to Maximum Results 

There appears to be no definite order in the relative vertical positions of 
the experimental lines in Fig. 6B. At a given value of 7, the ratio of t,, values 
might be expected to vary as D- (equation (4)) or as D-}74 (equation (5)) 
but this is not so. Wedge sample A3 is of the same material as wedge sample 
A2 but the ¢,, values are roughly a third higher. For the p-type sample with 
a diffusivity about 4 times and a lifetime about 17 times as great, the tp 
values are only about half those for sample A2. If an attempt is made to 
calculate values of tf, for the n-type material using equations (4) and (5), it 
is found taking 7 = 6 microseconds and D = 6 cm?/sec that the values are 
larger than the experimental ones by roughly a factor of 2. 

The empirical equations (12) to (17) show that the variation of ¢, with r 
is neither quadratic as in equation (4) nor linear as in equation (5). Further- 
more, detailed examination shows that the variation cannot be fitted by the 
more general relation (6), which applies not only to an infinite solid but also 
to a semi-infinite solid with s = 0 on the surface. In fact, even the theory 
with a non zero s in a semi-infinite solid is inadequate. This can be seen from 
Fig. 10 where r?/4D, tm, obtained from the experimental results on flat shiny 
sample A3 and wedge sample A2, has been plotted against f,. Equation (7) 
may be rewritten in the form 


r bee 
4D,h, ~ 7 tgte 


Thus the plot should have a slope 1/7 and an intercept on the ordinate axis 
of 3/2+.«. It is, however, apparent that the slope is not constant and the 
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Fic. 10. Variation of r?/4D,tm with time to maximum obtained from the experimental 
tm, r plots for the flat shiny surface sample A3 and the wedge sample A2, taking D, = 6 
cm?/sec. The broken curves were obtained by taking reciprocals of slopes of the continuous 
curves. 


extrapolated intercept does not lie between the values 3/2 and 5/2 as pre- 
dicted. In these considerations ¢ varies sufficiently slowly with ¢, to be taken 
as a constant. The theory with a constant lifetime is therefore clearly not 
sufficient. The variation of the reciprocals of the slopes, indicated by the 
broken lines in Fig. 10, shows a “‘lifetime’’ which apparently increases with 
tm; a tendency to flatten off somewhat is seen at the higher values. 

Another approach at analyzing the results is to use the one-surface theory 
to calculate s from equation (10). Unfortunately, such a calculation yields 
negative values. This is due to using the measured value for 7 of 6 micro- 
seconds which does not satisfy inequality (11). To get positive values of s not 
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only must a smaller value of 7 be used but also one which varies with fm. 
Figure 11 shows a plot of X-' and Y~ (see equations (3), (3A), (9), and (9A)) 
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Fic. Il. Variation of X~! = 4Dy tn?/(r?—6Dp tm) and Y! = 4Dp tm?/(r?—10Dp tm) with time 
to maximum obtained from the experimental tm, 7 plots for the flat shiny surface sample A3 
and the wedge sample A2, taking D, = 6 cm?/sec. The quantities X~! and Y~! represent 
values of lifetime calculated on the assumption of a zero and infinite surface recombination 
velocity, respectively. 


against ¢, for the same two samples as in Fig. 10. The two curves for each 
sample represent values of 7 calculated for the two extreme cases of zero 
and infinite s. To get positive finite values of s, by inequality (11), the calcu- 
lated values of lifetime must lie between the two curves. Since both pairs of 
curves show, for the most part, an increase, the lifetime must also increase 
with ¢,, if the theory is correct. This result is thus the same as that deduced 
by the previous type of analysis. However, the magnitudes and the shapes 
of the curves in Figs. 10 and 11 are not the same. 


Decay Time Results 

If the decay time 7; is taken as a measure of the lifetime, then the increase 
of 7; with ¢, in Fig. 8 is consistent with the variation already deduced from 
the analysis of the t,, 7 curves. However, the magnitudes of 7; are roughly 
3 and 5 times those of Figs. 10 and 11 respectively. From an analytical point 
of view, if empirical equation (18) is inserted in equation (5), it is found 
that ty = 3.6 10-*r'. Comparing with equation (16), namely f, = 0.9 107% 
r5, this shows a reasonable agreement in the exponent but a disagreement 
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by a factor of 4 in the proportionality constant. This limited amount of 
agreement is surprising when one considers that the substitution procedure is 
inherently wrong, since the derivation of equation (5) was based anyway on 
the assumption of a constant lifetime. 


DISCUSSION 


While the ‘‘lifetimes” calculated in the previous section differ depending on 
how they are deduced, they all show a “‘lifetime’’ which decreases as the 
emitter-collector spacing is reduced. The differences are no doubt due to the 
inapplicability of the theory. For consistency, the theory should include in 
the original differential equation (1) a lifetime which depends on distance or 
concentration. 

One possible way in which the lifetime could depend on concentration is 
as follows: It is estimated that about 10'° holes were injected during a pulse 
of 0.1 microsecond. This means that while the carriers were confined to a 
sphere of diameter 10~* cm, the local concentration was about 10!* holes/cm*. 
Direct electron-hole recombination can take place across the gap at a rate 
proportional to np/n{ (van Roosbroeck and Shockley 1954), where , is the 
intrinsic density of electrons. If Ap is much larger than the density of minority 
carriers (po) but much less than that of the majority carriers (mo) i.e. po K Ap 
< mo, then 7 is a constant equal to 2,7,/no, where 7, is the lifetime in intrinsic 
material; but if Ap is so large that it is even greater than the majority equili- 
brium density, i.e. Ap > mo, then the decay of Ap with time becomes hyperbolic 
instead of exponential and the “‘lifetime’’ is found to be equal to the time 
itself. In this case, the “‘lifetime’’, starting from the bulk value, would decrease 
progressively on approaching the emitter point. For the m-type material used, 
this decrease would start when Ap exceeds mp = 10"/cm‘, that is, if bimolecular 
recombination were the only factor. 

A lifetime dependent on excess concentration might also explain the fy, 7 
results on the 120-ohm-cm p-type silicon, which are particularly puzzling. 
This material shows the same type of dependence of t, on 7 as in the case 
of the 6-microsecond n-type silicon. Yet the p-type material has a 100-micro- 
second lifetime which is much larger than all the measured t, values. In this 
case we should have expected t, « r? instead of the observed t, « r'*, The 
magnitudes of f, are also smaller than expected for such a lifetime. These 
facts could possibly be accounted for with a much smaller effective lifetime 
governed by the excess carrier concentration. 

It should be mentioned that the assumption of a lifetime dependent on 
excess Carrier concentration has two serious difficulties. The first of these is 
that the decay following the maximum would not be exponential. Hence the 
dependence on concentration, if it exists at all, cannot be too strong. Secondly 
the lifetime at a given value of r should decrease with increased injection of 
carriers, whereas in fact no such decrease was observed. This is not understood. 

It should be mentioned here that Bath and Cutler (1958) have found that 
s in silicon depends on excess concentration but they do not believe body 
lifetime was affected. 

Another factor besides lifetime which could be affected by excess carrier 
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concentration is the diffusivity which arises in the original differential 
equation through Fick’s law. For large added hole concentrations Ap = An 
= p =n and the ambipolar diffusivity consequently reduces to D’ = 2D, 
(1+D,/D,)—. Since D, > D,, D’ is approximately 2D,. Thus D’ could decrease 
from 2D, near the point of injection to D, a long distance from it where the 
hole concentration is low. Besides this variation, an increase with distance 
may be obtained if D, and D, are reduced by electron-hole scattering at the 
high excess concentration near the emitter point. Such a variation would 
affect the tm, 7 curves and the calculated lifetimes. 

Trapping (Haynes and Hornbeck 1955) is not thought to be an important 
effect in this work. While it could explain the values of decay time greater 
than 6 microseconds, it could not explain the many observed times less than 
this value. 

The scatter of the experimental points in Figs. 4, 5, and 6 is somewhat 
larger than can be accounted for by an estimation of errors in the measure- 
ments of t, and r. The same is true for Figs. 8 and 9 but to an even greater 
extent. For these plots one naturally thinks of the capacitance and resistance 
of the collector point contact as uncontrolled parameters causing the scatter. 
However, these quantities have values which are apparently too small to 
affect the results. If indeed such factors had been important, there might 
have been a possible, though unlikely, explanation for the variation of decay 
time with 7; it was found necessary in the experiments to increase the pressure 
on the points as 7 was increased and this in turn could result in a change of 
capacitance and resistance of the collector point with distance. 

It is unfortunate that no two-surface theory is available for comparison 
with the experimental results in this paper, although the one-surface theory 
of van Roosbroeck does correspond to the flat surface experiments. However, 
the indications are that a two-surface theory alone would anyway be in- 
adequate. What is required is a theory which takes into account at the outset 
enhanced recombination and modification of the diffusivity near the emitter 
point. 
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NOTES 


STRAIN IN CRYSTALS OF YTTRIUM IRON GARNET 
ARCHIBALD W. SMITH AND GERALD W. WILLIAMS 


Yttrium iron garnet (Y;Fe;O,2) is a ferrimagnetic oxide which can be pre- 
pared in the form of small crystals by growth from a melt (Nielsen and Dear- 
born 1958). In a recent paper, Smith and Williams (1960) pointed out that 
an unusual type of strain occurs near (112) growth faces of these crystals. The 
strain axis, which is a direction of easy magnetization, is inclined to the growth 
face, and appears to be introduced during the growth of the crystals. This is 
to be contrasted with the strain which is often introduced during grinding 
of the crystals, and results in a strain axis perpendicular to the ground surface. 
As an extension of the work on (112) growth faces, the magnetic strain aniso- 
tropy has been investigated by means of ferromagnetic resonance. This method 
is convenient with yttrium iron garnet because of its unusually narrow resonance 
line width. The measurements have been carried out on an optically trans- 
parent disk, and the optical birefringence has been investigated with a polariz- 
ing microscope. The magnetic and optical anisotropy due to the strain have 
been correlated by considering the magnetoelastic and photoelastic inter- 
actions. An explanation of the magnetic birefringence is proposed in terms of 
the magnetostrictive and photoelastic effects. 


EXPERIMENTAL 

Two samples cut from the same crystal were investigated. They were 
prepared by grinding and polishing thin slices from the crystal, and cutting 
disks from these with an ultrasonic ‘‘cookie cutter’. Sample 1 had a (112) 
growth face on one side, and was 0.103 cm in diameter and 0.0023 cm thick. 
The crystal was strained, with the strain axis inclined to the growth face. 
The [111] direction was located in the plane by X-ray techniques. The domain 
structure took the form of strips running across the sample. The magnetization 
was parallel to the strain axis, and was antiparallel in adjacent domains. The 
strips were 2° from [II1] toward [110]. The sample exhibited a permanent 
birefringence with the direction of vibration of the slow wave 6° from [II1] 
toward [110]. 

The direction of minimum anisotropy energy in the plane of the disk is 
[111] for the intrinsic negative cubic anisotropy, and along the projection of 
the strain axis for the strain anisotropy. It will be seen below that the latter 
is nearly uniaxial in form. As the domain strips indicate the position of the 
composite energy minimum, it follows that the projection of the strain axis, 
which must coincide with the direction of vibration of either the fast or slow 
wave, was 6° from [111]. The strain axis itself was inclined at 60° to the 
plane. Thus the strain axis was close to the [001] direction, which occurs at 
55° above [111]. The inclination was determined by tilting the sample about 
an axis perpendicular to the domains until a position was found where the 
contrast between the domains disappeared. This occurs when the light path 
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inside the sample is parallel to the strain axis. The procedure for this experi- 
ment was the same as previously described (Smith and Williams 1960). 

The retardation due to the birefringence was measured with the magnetiza- 
tion parallel to the [111] and [110] directions. The results over the main part 
of the sample were 7.5X10~7 and 10.3X10~ cm, respectively. In a sharply 
defined region occupying about 15% of the disk, the results were 11.0X10~7 
and 14.2X10-7 cm. As the inclination of the strain axis was the same in this 
region, these results indicate that the magnitude of the strain was greater. 
This region was located along the edge of the disk in the [111] direction from 
the center. The difference between the two values for each region gives an 
average retardation of 1.510~7 cm for the magnetic birefringence, and shows 
that the direction of vibration of the fast wave is parallel to the magnetization. 

Sample 2 was cut from material about 0.025cm below a (112) growth 
face. It was ground on both sides and had a strain axis perpendicular to its 
plane. The diameter was 0.127 cm and the thickness 0.0036 cm. The domain 
structure took the form of an irregular maze pattern. The magnetization was 
parallel to the strain axis, and was antiparallel in adjacent domains. 

The resonance field measurements for the two samples are shown in Fig. 1. 


f = 9306 Mc/s 
T= 25°C 
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RESONANT FIELD 
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ANGLE BETWEEN FIELD AND [iil] DIRECTION 
Fic. 1. Resonant field plotted against angular position in plane of disk. The solid points 


were taken between 0 and 90°, and the open points between 90° and 180°. The curves are 
theoretical functions which have been fitted to the experimental points. 
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These were obtained using a reflection cavity technique similar to that 
described by Artman and Tannenwald (1955). The disks were mounted so 
that they could be tilted to align the plane of the disk with the field, and 
rotated to vary the angular position of the field in the plane. The field was 
measured with a nuclear magnetic resonance gaussmeter. 


THEORETICAL 
The resonant frequency w of a ferrite disk with the field in the plane of the 
disk is (cf. Artman 1957) 


9 


(2) = HH] 4M 30)+ 3 Esa Ea) [+44 (1-30) Eee 


where H is the external field, 6 and @ are the polar and azimuthal angles, 
respectively, of the magnetization M, and y is the gyromagnetic ratio. E, is 
the anisotropy energy, and E,ge and Eggs are the second order partial deriva- 
tives with respect to 6 and ¢. This expression is correct to the first order in 
the angular deviation of M and H. In the resonance measurements described 
here, this deviation is less than 2°, and it will be assumed that M and H are 
aligned. 

The anisotropy energy is the sum of the magnetocrystalline and magneto- 
elastic or strain contributions. The former can be adequately expressed by a 
first-order term only (Dillon 1957). The latter will be written for a uniform 
strain ¢, in terms of the two magnetostrictive constants Ajoo and Aq (cf. 
Baltzer 1957). The total energy is then 


: i 2 3 ee 
E,=K a aia; —Sd100(Cur— Crs)e aiyi—6rinCue >, LV LY; 
i>j “ i i>j 


where the a; and y; are the direction cosines of the magnetization and strain 
axis, respectively, referred to the edges of the unit cell. K is the magneto- 
crystalline anisotropy constant, and Cy, Ciz, and Cys are the stiffness con- 
stants. If the strain is along a cube edge, or if the crystal behaves isotropically, 
i.e. AqaiCsa = ¥A100(Cur— Cz), the strain contribution reduces to a uniaxial 
term —K, cos? 6, where K, = (3/2)A100(Ci1—Ci2) and 6 is the angle between 
the magnetization and the strain axis. The experimental results can be 
adequately represented by taking a uniaxial term only, with the projection 


of the strain axis along [111]. In terms of the co-ordinate system of Fig. 2, 


STRAIN 
AXIS 


DISK IN 


(iin x (2) PLANE 


Fic. 2. Co-ordinate system for disk in (112) plane. Polar angles are measured from [112], 
azimuthal angles from [111]. 





NOTES 


E, = K[sin‘ 0(4 sin’ 6+4 cos’) +4 cos* @—“ sin 6 cos*@ sin 
+sin’ 6 cos @ sin” (4 cos 0+ +/2 sin 6 cos ¢)] 
—K,(sin 7 sin @ cos ¢+cos 7 cos 6)’, 
and the partial derivatives are 
Exoo = K(sin® $ cos’ ¢—} cos’ ¢) —2K,{1—sin® 7(1+cos’ ¢)] 


and 
Engs = K(7 cos’ ¢ sin” ¢—sin* ¢—}cos* ¢)—2K, sin’ r(1—2 cos” ¢). 


Examination of the above equations shows that K, and 1, the polar angle 
of strain axis, enter only in the form 4%M(1—3n) —2K,/M and (2K,/M)sin? r. 
This means that K, cannot be determined independently of 4r%/(1—3n) and 
7 by resonance measurements alone. 

The permanent birefringence which occurs in sample 1 is due to the photo- 
elastic effect of the strain. This effect has been evaluated in terms of the 
elasto-optical coefficients for a uniform strain along the [001] direction, which 
is nearly the case in sample 1. The index of refraction for light passing through 
the strained crystal in the [112] direction is given by the index ellipse 


cos” ¥ sin” y 1 
oe 
ny ny 


where y is the angle of the plane of polarization with respect to [II1]. In an 
unstrained crystal ny = n, = n. In terms of the elasto-optical coefficients P1; 
and Py. (Nye 1957), my—n, = n*(P1:— P12)€/6. The measured retardation is 
equal to d(my—n,), where d is the thickness of the sample. 

The ratio of the magnetic and optical anisotropy is independent of the 
strain: 


3 


ie - 9(Cu— Cis) A100 | 
nyrn, n Py—P 2h" 


In order to establish the validity of the theory, this ratio will be evaluated 
and compared to the experimental result. The only quantity which has been 
measured is the index of refraction: » = 2.24. The others may be estimated 
as follows: Ayo0 = 5X10~° (cf. Baltzer 1957, Table II), Ci: —Ci2 = 10" dynes/ 
cm? (cf. A.I.P. Handbook 1957, Table 2f-1), and Py:—Pi. = 0.2 (cf. Nye 
1957, Table 17). The result is K,/(my—n,) = 2X10? dynes/cm?. 

The magnetostrictive and photoelastic effects provide a possible explanation 
of the magnetic birefringence. In contrast to the Faraday rotation, the magnetic 
birefringence in yttrium iron garnet is only weakly dependent on wavelength, 
which suggests that it has a different origin. The Faraday rotation arises 
from the effect of the internal magnetic field on electric dipole transitions 
in the Fe** ion (Clogston 1959). The photoelastic birefringence arising from 
the magnetostriction has been evaluated for an isotropic medium, with the 
magnetization perpendicular to the line of sight. Under these conditions, the 
strain is uniform and either the major or minor axis of the index ellipse is 





772 CANADIAN JOURNAL OF PHYSICS. VOL. 39, 1961 


parallel to the magnetization. In an anisotropic medium this would not 
generally be the case. The difference in the principal indices of refraction is 
ny—n, = $dn*(P1:—P 12), where \ is the isotropic magnetostriction. Using 
the same data as before, this quantity is estimated to be |n,—n,| = 810-6, 


DISCUSSION 


The resonance equation has been plotted in the form H vs. ¢ for sample 2, 
with 7 = 0 and taking K/M = —42 oersteds and g = 2.005 (y = ge/2mc 
= 2.8064 Mc/sec oersted) from Dillon (1957). The value 4rM(1—3n)—2K,/ 
M = 1645 oersteds was derived from the point at @¢ = 0. It can be seen 
from Fig. 1 that this gives a satisfactory fit to the experimental points. As 
the sample is appreciably strained, this means that the crystal is nearly 
isotropic in behavior, i.e. Ay00(Cir— Cw) = 2Ai11C44. The fact that the strain 
anisotropy overrides the cubic indicates that 2K,/M > 80 oersteds; this 
inequality cannot be calculated exactly because the first derivative of the 
cubic anisotropy is not zero in the [112] direction. A value of 2 = 0.021 was 
calculated on the assumption that the sample was an ellipsoid, although this 
is not an adequate representation. A value of 4rM > 1840 gauss is then 
found. This is high compared to most published values, which fall in the 
range 1720+40 gauss, but is compatible with the value of 1890 gauss recently 
found for a particular batch of crystals (Rossal 1960). 

For sample 1, the resonance equation has been plotted by deriving 
4rM(1—3n)—2K,/M = 1500 oersteds and (2K,/M) sins = 76.7 oersteds 
from the points at @¢ = 0 and 90°. The fit is not as good in this case, but is 
still satisfactory considering that the sample is not uniformly strained. The 
strain anisotropy is expected to be uniaxial in this case because the strain 
axis is close to [001]. Taking r = 30° and m = 0.017, the values 2K,/M = 300 
oersteds and 4rM = 1900 gauss are obtained. The difference between the 
two values of 4x may be due to an inadequate analysis of the non-uniform 
sample, or may indicate an actual difference between the growth face and 
the inner part of the crystal. 

In evaluating the ratio |K,/(my—n,)| experimentally, it will be assumed 
that the resonance experiment measures a weighted average of the strain 
anisotropy in the two regions of sample 1. The retardation calculated on the 
same basis is 1.0X10-*cm. The experimental ratio is then |K,/(my—n,)| = 
4.8X10° dynes/cm?, which is in satisfactory agreement with the predicted 
value. The corresponding value of the strain ¢ is about 2X10-*. Thus it is 
possible to make a unified analysis of the magnetic and optical properties of 
the strained crystals. For the magnetic birefringence, |my—m,| = 6.51075 
experimentally, which is about 10 times the predicted value. It is clear that 
the mechanism outlined here will make an appreciable contribution to the 
magnetic birefringence, although it may not be able to account for all of 
it. Measurements of \ and P:;— P12 would be required to make a more con- 
clusive analysis, as well as more extensive experiments on the magnetic 
birefringence itself. 
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THE K INTERNAL CONVERSION COEFFICIENT OF THE 
412-KEV TRANSITION IN Hg!** 


J. L. WoLFrson 


Hamilton, Frey, and Hultberg (1960) have recently measured ag, the K 
internal conversion coefficient of the 412-kev transition in Hg!®* by the inter- 
nal—-external method (Hultberg and Stockendal 1959a; Hultberg 1959) arriving 
at a value of (2.44+0.08)X10-*. Using the same method, de Vries et al. 
(1960) measured the value to be (2.52+0.13) X10-*, whereas employing the 
more common method of comparing the area of the K line to that of the beta 
spectrum, de Vries et al. arrived at the value (2.41+0.09) X10-*. These values 
are all lower than most determinations of this quantity (a complete listing is 
given in Nuclear Data Sheets). The most accurate determination, made by 
Wapstra et al. (1959), is (2.81+0.05) X10-*. 

Hamilton et al. point out that their value is considerably lower than the 
theoretical value of 3.0010-? (Sliv and Band 1956; Rose 1958) and regard 
it as anomalously low, since the total conversion coefficients of other £2 tran- 
sitions, in low Z elements, viz. Co®® (Frey, Hamilton, and Hultberg 1960) and 
Sc** (Hultberg, Frey, and Hamilton 1960), measured by the same method, 
have been found to be in good agreement with the theoretical values of Rose 
(1958). They conclude that the latest theoretical E2 conversion coefficients 
of high Z unhindered transitions may be in substantial error. 

It is the purpose of this communication to draw attention to hitherto 
unreported measurements made some years ago at the Chalk River Labora- 
tories of Atomic Energy of Canada Limited, which have a bearing on this 
question. In these measurements ax for the 412-kev transition was measured 
relative to a7, the total internal conversion coefficient, for the 1.33-Mev and 
1.17-Mev transitions in Ni®, whence it becomes possible to deduce the ratio 


*Issued as N.R.C. No. 6256. 
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of ax, for the 412-kev transition to each of the values of ag for the 1.33-Mev 
and 1.17-Mev transitions in Ni®. The method employed avoided comparisons 
of line areas with integrated beta spectra on the one hand, and on the other, 
achieved a comparison of the internal conversion ccefficients of an E2 tran- 
sition in a high Z nuclide (412 kev in Hg!**) with those of two £2 transitions 
in a low Z nuclide (1.33 an] 1.17 Mev in Ni®). The results obtained do not 
support the findings of Hamilton, Frey, and Hultberg. The ratios obtained 
are found to be somewhat higher than the theoretical ratios. 

The measurements made involved the determination of two quantities: the 
relative numbers of conversion electrons from the transitions in question, and 
the numbers of quanta. The former were determined by comparison of the 
areas of the corresponding internal conversion lines in the conversion electron 
spectra from Au'®’ and Co® sources, using a magnetic lens spectrometer. The 
latter were deduced from the measured disintegration rates of the same 
sources. This procedure avoids difficulties inherent in determining intensities 
of beta-ray spectra, but there is an assumption of constant spectrometer 
transmission over the energy range 328 kev to 1.32 kev. There are experi- 
mental checks which justify this assumption in that determinations of the 
intensities of conversion lines obtained with the spectrometer have been found 
to be in good agreement with identical measurements made elsewhere on 
other instruments. As an example, the intensity of the line at 9988 gauss cm 
to that at 2607 gauss cm in the conversion electron spectrum of RdTh active 
deposit, 0.099+0.003 as found with the present instrument (Elliott e¢ al. 
1954), compares with the value 0.101+0.002 as found by Krisiuk et al. (1956). 
This example is quoted as being particularly pertinent since in the latter 
measurement the question of penetration of counter slit edges by the energetic 
conversion electrons received careful study. 

The source of Au'®* was prepared by neutron irradiation of gold leaf of 
thickness 200 ug/cm*. The activity of the source was measured by direct 
comparison with a standard Au'®’ source supplied by the standards laboratory 
of Atomic Energy of Canada Limited. The comparison was made by gamma- 
ray counting using a sodium iodide detector. The activity of the standard 
source had been measured to better than 1% and the comparison was made 
with an accuracy of better than 1%. In deducing the number of 412-kev 
quanta from the disintegration rate, a value of a7 = 0.041 was assumed. The 
findings of Hamilton et al. (1960) and de Vries (1960) for ag would lead to a 
value of a7 = 0.036, with a resulting uncertainty, on this account, of about 
3% in the number of 412-kev quanta. The total uncertainty in the number 
of 412-kev quanta is estimated to be less than 2%. 

Sources of Co were prepared by evaporation from Co(NQOs)2 solution. 
Average source thickness was about 200 ug/cm*. Because the transitions 
studied were of high energy, and also because of the high specific activity 
of the Co® solution, uniformity of deposition was not of primary importance. 
Nevertheless since uniformity could not be assured by this method of source 
preparation, three separate sources were prepared and measurements made 
with all three. As for the Au!%* source the disintegration rates of the sources 
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were measured by comparison with a standard source, supplied by the Atomic 
Energy of Canada Limited standards laboratory. The number of 1.33-Mev 
and 1.17-Mev quanta from each source was deduced with an accuracy of 
about 13%. 

The conversion lines from the Au!*§ source and one of the Co® sources 
(source no. 2) are shown in Fig. 1. For the Ni® conversion lines, the L and 
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Fic. 1. Internal conversion lines of the 412-kev transition in Hg'® and of the 1.17- and 
1.33-Mev transition in Ni®, The vertical bars represent standard deviations of the plotted 
points. 


M lines are not resolved from the A lines. Using the data of Kamada et al. 
(1958), who found that for each of the transitions the ratio K/(L+M) was 
9.1+0.5, it was possible to make a correction for the Z and M line contribu- 
tions to the line areas. This leads to a value of ax = (0.900+0.005)a, for 
each of the transitions, where a7 = ax+az+tay. (It should be noted that 
using the theoretical values of Rose (1958) for the internal conversion co- 
efficients, and performing the necessary extrapolations, leads to a value 
ax = 0.888a7.) 

One further correction was required. The three Co® sources had diameters 
respectively of 2.35, 2.30, and 1.75 mm whereas the Au'®* source shown in 
Fig. 1 had a diameter of 2.30 mm. In order to determine the effect of source 
diameter on line area, Au!%’ sources of diameters 1.08, 2.30, 3.65, and 4.00 mm 


were prepared, their relative activities measured, and the area of the y 412- 

kev K lines from each of the sources measured in the spectrometer. It was 

found that for diameters close to 2 mm the ratio of line area to source activity 
992 


varied only slightly with source diameter. Relative to the 2.30 mm diameter 
sources, corrections of —1.4° and 0°% were required for the 1.75 mm and 


2.35 mm sources respectively. 
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The results obtained cannot be compared directly with the theoretical 
values of Sliv and Band (1956) since these authors have not computed values 
for Z < 33. Nor can they be compared directly with the theoretical values 
of Rose (1958) since Rose has not given values for energies >2.0 mc?. Extra- 
polation of the computed values of Rose indicates, however, that the theoretical 
values of ag for the 1.33-Mev and 1.17-Mev transitions are expected to be, 
within about 1%, 1.1210-* and 1.48X10-* respectively. The theoretical 
value of ag for the 412-kev transition in Hg! is 3.0010- (Sliv and Band 
1956). 

The ratios of ag for the 412-kev transition in Hg!®8 to each of the values 
of ag for the 1.33-Mev and 1.17-Mev transitions in Ni®, i.e. ag(Hg!®*)/ag(Ni®), 
were found to be: 


Co® sources 


Transition 
in Ni® (Mev) No. 1 No. 2 No. ¢ Means 





1.33 280 290 7 283412 
1.17 217 223 —_ 220+10 





A comparison with the theoretical values, and with the results of Hamilton, 
Frey, and Hultberg (1960), gives the following for the same ratios: 








ag(Hg!%, Hamilton et al.) 
Ni® (Mev) This work Theoretical axg(Ni®, theoretical) 


Transition in 


1.33 283+12 268 218+7 
1.17 220+ 10 203 165+5 





The results do not support the tentative conclusions of Hamilton e¢ al. 
that ag for the 412-kev transition in Hg!®8 is anomalously low. The experi- 
mentally determined ratios found in this work are actually somewhat higher 
than the theoretical ratios. 

The results do not of course lead directly to a value for ax for the 412-kev 
transition in Hg!*%’. If the theoretical values of ag for the 1.33- and 1.17-Mev 
transitions in Ni®® are combined with the results, values of (3.17+0.16) X10-? 
and (3.25+0.16) X10- are obtained for ag for the 412-kev transition. Most 
experimental values (Nuclear Data Sheets) of ag for the Ni® transitions are 
higher than the theoretical ones, leading to even larger values for ag for the 
412-kev transition. De Vries (1960) has pointed out, however, that if Hultberg 
and Stockendal’s (19596) data are combined with Colgate’s (1952) value for 
the photoelectric cross section of uranium at 1.33 Mev, then a value of 
0.98X10-* (46%) is obtained for ag for the 1.33-Mev transition. If this 
value be combined with the corresponding experimental ratio above, then a 
value of 2.78X10-? (47%) is obtained for ag for the 412-kev transition in 
Hg!®8, in agreement with the determination of Wapstra et al. (1959), of 
(2.81+0.05) X 10-. 
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Frequency Measurement of Standard Frequency Transmissions! 


Measurements are made at Ottawa, Canada, using N.R.C. caesium-beam frequency resonator 
as reference standard (with an assumed frequency of 9 192 631 770 c.p.s.). Frequency devia- 
tions from nominal are quoted in parts per 10". A negative sign indicates that the frequency 
is below nominal. 





GBR, 16 kc/s 
Date, —— 
February MSF, 63-hour 24-hour WWVB, 
1961 60 kc/s average* averaget 60 kc/s 

1 — 160 —151 —152 N.M. 
— 158 —150 —151 N.M. 
— 162 —151 —153 —155 
—148 —153 — 152 — 150 
—153 — 146 — 148 N.M. 
— 155 —151 —150 N.M. 
— 160 —150 — 149 — 156 
— 156 — 143 — 144 —151 
—155 —149 — 149 — 153 
—152 — 143 — 146 — 162 
—146 —144 —144 —147 
—146 — 142 —142 N.M. 
— 150 — 143 —141 —150 
N.M. —141 — 144 —150 
— 148 —140 —141 —153 
—152 —137 — 138 —154 
— 157 —144 — 143 — 157 
—155 —140 — 143 — 148 
— 143 —143 —141 —145 
—155 —142 — 142 — 156 
—155 —145 —145 —151 
— 158 —142 — 143 —152 
—162 —145 — 147 — 153 
— 160 —147 N.M. — 154 
—152 — 148 —150 —147 
N.M. —151 —148 N.M. 
—151 — 147 —152 —152 
—150 —142 —145 — 152 
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Average — 154 —145 — 146 — 152 


Midmonthly 

mean of WWV —150 
Note: N.M. no measurement. 
*Time of observations: 23.30 to 06.00 hours U.T. 
tTime of observations: 15 to 15 hours U.T. 
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The Conversion Electron Spectrum from a Mass-separated Source of Xe'*%* 


A mass separator of the Nobel Institute type (Alvager and Uhler 1958) has recently been 
installed and brought into operation in the Chemistry Division of this laboratory. It is designed 
primarily to produce isotopically pure radioactive materials for sources in nuclear spectroscopy 
or for radiochemical investigations. It can also provide small quantities of stable material of 
high isotopic purity for use, for example, in targets for nuclear reaction studies. 

As a test of the performance of the mass separator, we have prepared a source of Xe", 
1 mm wide by 20 mm long, for use with the Chalk River iron-free 8-ray spectrometer (Graham 
et al. 1960). An iron ampoule containing natural Xe gas was irradiated for 6 days in a neutron 
flux of 2.510" cm™ sec™! in the NRU reactor. We are indebted to Dr. R. O. Kelly for supply- 
ing the container and for arranging this irradiation. The gas was then admitted to the ion 
source of the mass separator at an accelerating potential of 30 kv. The focused Xe ions appeared 
at the collector as lines ~1.5 mm wide and were separated by ~15 mm per mass unit. Under 
optimum conditions the line width can be reduced to less than 1 mm. A 1-mm X20-mm slit 
defined the fraction of the Xe! ion beam striking a source foil of 800 wg/cm? Al. There was 
no deceleration voltage at the collector and it is known from other experiments (Davies 1961) 
that the mean depth of penetration of the ions into the collecting foil was ~8 yg/cm?. 

The conversion lines of the 80.99+0.02 kev y-transition in Cs! following B-decay from 
Xe! are shown in Fig. 1. The 8-spectrometer was adjusted for an instrumental resolution of 
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Fic. 1. Conversion electron spectrum of the 80.99-kev gamma transition in Cs!33, The points are the observed 
counts which include 8-continuum and background rates. 


~0.07% in momentum. The greater line widths observed experimentally are attributed to 
source thickness, i.e., to energy loss as the electrons emerge from the depth at which the Xe 
atoms are located. 

The ratios of the line intensities are: K:L,;:Le:L3 = 79+2:10.0:1.48+0.05:1.15+0.05. 
Comparison with the theoretical conversion coefficients of Sliv and Band (1957) indicate that 
the multipolarity of the 80.99-kev transition is 97.4% M1+(2.6+0.2)% E2. From this and 
the measured half-life of (6.00.4) X10~® sec for this transition (Graham and Bell 1953) one 
deduces that the transition rate of the E2 component is about five times the single particle 
estimate (see, for example, Wapstra et a/. 1959). If one takes into account the statistical factor 
(Moskowski 1953) of S = 0.19 for the spin sequence 5/2 7/2 (Bergstrém 1952) the 
enhancement factor is ~25. The hindrance factor for the 171 component is ~700. 
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